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INTRODUCTION 


A  computer  simulation  was  developed  for  a  complete  artillery  safing  and 
arming  (S&A)  mechanism  containing  a  straight-sided  verge  runaway  escapement,  a 
two  pass  Involute  step-up  gear  train,  and  a  spin  driven  rotor  which  must  operate 
in  a  projectile  that  experiences  spin,  precession,  and  nutation. 

Top  views  of  the  mechanism  planes  of  the  two  possible  configurations  are 
shown  in  figures  1  and  2.  The  position  of  this  mechanism  plane  with  respect  to  a 
projectile  which  experiences  the  aeroballistic  motion  is  shown  in  figure  3. 

While  the  basic  ideas  concerning  the  three  motion  regimes  of  the  runaway 
escapement  (coupled  motion,  free  motion,  and  impact)  are  identical  to  those 
developed  for  the  verge  type  runaway  escapement  (ref  1),  the  presence  of  the 
three  dimensional  projectile  kinematics  makes  it  necessary  to  consider  three 
dimensional  force  and  moment  equations  for  all  mechanism  components  in  order  to 
derive  the  mathematical  models  for  the  motion  regimes. 

The  following  briefly  outlines  the  course  of  the  derivation  of  the  mathemat¬ 
ical  models: 

1.  Kinematics  of  Aeroballistic  Bystems:  Absolute  angular  velocities  and 
accelerations  in  terms  of  component-fixed  and  projectile-fixed  coordinate  systems 
(app  A) 

2.  Angular  Momentum  and  Its  Derivatives  in  Various  Coordinate  Systems: 
Three  dimensional  moment  equations  in  various  coordinate  systems  (app  B) 

3.  Absolute  Acceleration  of  the  Geometric  Center  C  of  the  Mechanism 
Plane  (app  C) 

4.  Dynamics  of  Rotor  Driven  S&A  Mechanism  with  a  TV/o  Pass  Involute  Gear 
Train  and  a  Verge  Runaway  Escapement  Operating  in  an  Aeroballistic  Environment: 
The  derivations  of  the  equations  of  motion  for  both  entrance  and  exit  coupled 
motion,  free  motion,  and  impact  regimes  are  contained  in  the  appendix.  Expres¬ 
sions  for  all  types  of  contact  forces  are  given.  The  pivot  friction  forces  are 
treated  conservatively  (refs  1  and  2).  The  change  of  direction  of  the  friction 
forces  and  torques  in  the  gear  train  due  to  a  motion  reversal  of  the  mechanism 
are  handled  by  appropriate  sign  change  of  the  coefficient  of  friction  (app  D). 

5.  Projectile  Kinematics:  Since  at  the  present  time  actual  aeroballis¬ 
tic  data  are  not  available  for  incorporation  into  the  program,  a  set  of  appropri¬ 
ate  expressions,  which  allows  certain  simulation  runs,  has  been  provided  (app  E) . 


F  6.  Computer  Program  SAEROV:  The  listing  of  the  program  also  contains  a 

(J  sample  output  (app  F) . 

<  i 

;;  To  understand  the  derivations  in  the  appendixes,  it  is  suggested  that  refer- 

J  ences  1  through  5  be  consulted  concerning  the  work  on  gear  trains  as  well  as 

!'  -  rotor  and  constant  torque  driven  S&A  mechanisms  which  contain  pin  pallet  and 

verge  type  runaway  escapements.  For  general  background  and  kinematics,  reference  * 

' 


f 

g 


1 


ESCAPE  WHEEL 
&  PINION  3 


s! 

8 

I 

o 

1 

/ 

y 

N 

GEAR  & 

z 

o 

s 

s  _ 

f 

M  i 

_ ! 

3u  * 

KA 

oil 


/  ”  1 3  3 

/  sy 

/  o  2 

L  « 


0 


Figure  1.  Rotor  driven  S&A  device  with  verge,  configuration  no. 


Figure  3.  Mechanism  plane  in  projectile  which  experiences  aeroballistic  motion 


6  Is  recommended .  For  the  understanding  of  three  dimensional  rigid  body  dynamics 
as  well  as  the  relationships  between  the  non-or thogonal  coordinates  by  which  the 
aeroballistlc  motion  Is  expressed  and  the  orthogonal  ccmponent-f lxed  and  projec¬ 
tile-fixed  coordinate  systems  used,  references  7  and  8  are  suggested. 


DESCRIPTION  OF  COMPUTER  PROGRAM  SAEROV 


With  the  exception  of  the  inclusion  of  the  aeroballistlc  kinematics,  the 
programming  schemes  which  make  it  possible  to  distinguish  between  entrance  and 
exit  coupled  motion,  free  motion,  and  impact  run  parallel  to  those  given  in  ref¬ 
erence  1.  (This  reference  lists  the  control  details  applicable  to  the  separation 
of  entrance  and  exit  coupled  motion.  For  other  control  details,  see  references  2 
through  5.) 

The  main  program  starts  with  the  reading  in  and  writing  of  all  relevant 
physical  data.  This  is  followed  by  the  computation  of  gear  ratios,  fuze  body 
angles,  gear  train  constants,  and  earliest  and  latest  possible  values  of  the  gear 
angles  by  way  of  subroutine  ALFA,  as  well  as  the  initialization  of  the  gear 
angles.  The  simulation  begins  with  entrance-coupled  motion  at  a  starting  angle 
PHID,  which  represents  that  angle  $  of  the  escape  wheel  that  is  associated  with 
the  approximate  center  of  the  entrance  working  surface  of  the  pallet.  This  angle 
then  corresponds  to  a  cumulative  escape  wheel  angle  PHI TOT  of  0  degree. 


Coupled  Motion  (Location  1) 


Regardless  of  whether  entrance-  or  exit-coupled  motion  takes  place,  differ¬ 
ential  equation  D-513  must  be  solved.  (The  difference  between  entrance  and  exi- 
tion  motion  is  set  by  the  value  of  the  signum  function  Bj  as  used  in  the  computa¬ 
tions  of  parameters  Aj^,  Ao i  ,  A29,  A33,  A-j^,  and  A5j.)  To  this  end,  the  main 
program  calls  on  an  available  fourth-order  Runge-Kutta  routine.1  The  main  pur¬ 
pose  of  the  subroutine  FCT  is  to  present  the  second-order  differential  equation 
in  terms  of  two  first-order  ones  to  RKGS.  PHI(l)  and  PHI(2)  represent  the  angle 
<|>  and  the  angular  velocity  £,  respectively.  The  computation  of  all  parameters  of 
the  differential  equation  takes  place  by  way  of  subroutine  FCT,  which  itself 
calls  on  subroutine  KINEM  and  AFIVE.  The  latter  subroutine  calls  both  on  subrou¬ 
tine  ACCEL,  which  depends  on  subroutine  AERO,  as  well  as  on  the  sequential  sub¬ 
routines  AWON,  CWON,  ATWO,  CTWO,  A THREE,  C THREE,  AFOUR,  and  CFOUR. 

The  subroutine  KINEM  computes  current  values  of  g,  i|»  and  \p  (ref  l,  app  C)  as 
well  as  the  moment  arms  Aj ,  B| ,  C| ,  and  D|  (ref  1,  app  D) . 


*  RKGS  Routine,  IBM  System/360  Scientific  Subroutine  Package  (360A-CM-OX3) ,  Ver¬ 
sion  III. 
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Subroutine  AFIVE  computes  varioue  gear  mesh  parameters  and  instantaneous 
mesh  contact  angles,  as  well  as  the  signum  functions  sj,  82,  S3,  s^,  S3,  and  S7. 
In  addition,  the  parameters  Aj  to  Aj20  an<*  ^1  to  ^72  (aPP  D)  are  obtained  with 

the  previously  mentioned  subroutines The  gear-indexing  operation  (ref  A)  is 
performed  with  the  help  of  the  angle 

The  Instantaneous  rotor  angle  <J>^  4-  of  the  coupled  motion  differential 

equations  must  be  expressed  in  subroutine  rCT.  Recall  that  ij».r  is  the  initial 
rotor  angle;  ^  is  the  total  angle  of  rotation  of  the  escape  wheel.  Since  the 
angle  <{>  with  the  Runge-Kutta  variable  PHI(l)  varies  between  approximately  134  and 
144  degrees  during  entrance-coupled  motion  and  between  209  and  216  degrees  during 
exit-coupled  motion,  the  total  escape  wheel  angle  *  can  only  be  obtained  by 
continuously  adding  the  increments  due  to  each  cycle  of  Runge-Kutta  computations. 
Therefore 

♦r-W*4*  (1) 

where 

<(>  *  total  escape  wheel  angle  up  to  a  certain  Runge-Kutta  cycle. 

iu  (This  is  represented  by  PHITOT  in  the  program.) 

A<J>  -  increment  of  escape  wheel  during  this  Runge-Kutta  cycle. 

The  increment  A<p  is  calculated  as  the  difference  between  the  latest  value  of 
PHI(l)  and  its  previous  one  which  has  been  stored  as  PHIPR.  In  this  manner, 
equation  1  becomes 

♦T  -  PHITOT  +  PHI(l)  -  PHIPR  (2) 

Subroutines  associated  with  AFIVE  also  decide  on  the  values  of  IpR  and  I.R 
as  required  by  equations  D-134  and  D-135  as  well  as  equations  D-405  and  D-410 
(app  D).  The  associated  conditional  statements  assign  the  larger  values  for 
these  combined  moments  of  inertia  whenever  the  product  of  the  angular  velocity 
and  the  angular  acceleration  is  positive;  i.e.,  both  quantities  have  the  same 
sign. 


The  associated  subroutine  OUTP  is  responsible  for  printing  out  the  results 

s  **  * 

<J>,  and  <f>,  together  with  the  current  values  of  time,  g,  t|>,  41,  and  PHITOT. 
Further  all  coupled  motion  contact  forces  are  calculated  according  to  equations 
D-520,  J>-525,  D-527 ,  and  D-529,  and  the  maximum  values  of  these  forces  during  one 
arming  cycle  are  determined. 


^  The  program  uses  the  symbols  AA1,  etc.  throughout.  This  should  not  be  confused 
with  the  symbols  AAj^  to  AAjj  which  are  first  used  in  the  combined  exit-coupled 
motion  differential  equation  (D-278), 
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Free  Motion  (Location  5) 


The  differential  equations  of  free  motion,  as  given  by  equation  D-53U  for 
the  pallet  and  equation  D-531  for  the  combined  escape  wheel-gear  train-rotor 
system,  are  again  solved  by  the  Runge-Kutta  routine.  To  obtain  the  magnitudes  of 
the  variables  <p  and  ip,  as  well  as  their  derivatives  at  identical  times,  the  two 
independent  second-order  differential  equations  are  transformed  into  four  simul¬ 
taneous  first  order  equations.  (While  only  the  two  first-order  equations  associ¬ 
ated  with  each  of  the  two  variables  are  actually  coupled,  the  routine  treats  all 
four  as  if  they  were  coupled  and,  therefore,  produces  solutions  for  identical 
time  increments.)  These  four  expressions,  which  are  presented  in  subroutine 
FCTF ,  are  of  the  following  form: 

UX( 1)  -  X(2)  :  (-  (3) 
DX(3)  -  X(4)  :  (-  i*p)  (4) 


DX(2)  m-j~-  [rAilb(X(2))2-AuyX(2)  +  A^J  :  (-  <j>)  (5) 

DX(4)  - -J-  [-A  (X(4))2-A  X(4)-A  +  in  r  (X  sinB-X  cosfl)  ]  :  (“  *P)  (6) 

Iol}  1  32  31  116  p  cp  x  y 


The  subroutine  FCTF  also  calls  on  subroutine  AFIVE  for  the  computation  of  all 
gear-related  parameters. ^ 

The  associated  subroutine  OUTFF  computes  the  free  motion  contact  forces 
according  to  equations  D-537  and  D-53y  and  finds  their  maxima.  In  addition,  a 
continuous  count  of  P11IT0T  is  provided  by  OUTPF.  This  angle  as  well  as  time, 

<j>,  ip,  ^i,  and  the  contact  forces  are  printed  out.  The  same  routine  also  makes 
the  decision  of  whether  or  not  to  remain  in  free  motion.  The  sensing  variables  f 
and  g'  ■  (IP  are  used  for  this  purpose  (ref  1,  eqs  E-4  and  E-3). 


Impact  (Location  13) 


The  subroutine  IMPACT  uses  the  pre-impact  values  ^ ,  and  of  the  angular 

velocities  and  computes  their  post-impact  values  and  ip^  according  to  equations 
D-540  and  D-541  (app  D).  (Note  that  the  moment  of  inertia  of  the  escape  wheel  is 
now  expressed  according  to  equation  D-542  (app  D) ,  which  refers  the  rotor  as  well 
as  the  gear  train  inertias  to  the  escape  wheel.) 
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Whenever  1^  **  o,  the  simulation  stops  because  of  the  division  by  zero.  Should 
this  occur,  FCTK  prints  "IFK  EQUALS  ZERO— SIMULATION  TERMINATED." 


Reversal  of  Gear  Train  Motion  Due  to  Impact 


If  the  Impact  torque  on  the  escape  wheel  is  sufficiently  large,  the  motion 
of  the  gear  train  may  be  temporarily  reversed;  i.e.,  the  escape  wheel  angular 
velocity  $  may  become  negative.  This  would  cause  the  friction  forces  between  the 
gear  teeth  and  at  the  various  gear  pivots  to  be  reversed  in  direction.  (The 
normal  forces  between  the  gear  teeth  remain  unaffected,  and  the  normal  bearing 
forces  are  obtained  in  the  usual  manner.)  This  change  in  the  direction  of  the 
friction  forces  is  expressed  for  both  coupled  and  free  motion  by  letting  the 
coefficient  of  friction  y  of  all  gear  train  components  become  negative  (app  E  of 
ref  2).  This  is  accomplished  in  subroutine  APIVE  by  the  following  use  of  the 
signum  function  |/|$|  : 

MU  -  ABS(MU)  *  $/|$|  (7) 

(The  coefficient  of  friction  associated  with  the  escapement  interface  and  the 
pallet  pivot  is  called  y  and  is  read  into  the  programs  as  tiUl.)  Any  motion 
reversal  at  these  surfaces  is  accounted  for  by  the  signum  functions  s^  and  85, 
respectively. 


Termination  of  Computations 


Computations  are  terminated  whenever  the  geared  motion  of  the  rotor  ends. 
This  corresponds  to  <J>  •*  PHICUTD.  The  duration  of  the  subsequent  unretarded 
motion  of  the  rotor  is  assumed  to  be  negligible. 


COMPUTER  SIMULATION  Of’  EXAMPLE  MECHANISM 


The  mechanism  which  has  been  simulated  is  that  of  a  modified  S&A  device  of 
the  M/39  fuze.  It  has  configuration  no.  2  (fig.  2)  and  contains  a  newly-designed 
Involute  gear  train.  While  this  gear  train  has  the  same  gear  ratio  and  individ¬ 
ual  center  distances  as  the  original  design,  each  of  the  meshes  now  has  unity 
contact  ratio. ^  The  simulation  of  this  mechanism  was  accomplished  with  the  help 
of  computer  program  SAEKOV  (app  P) .  It  was  run  for  30,U0l)  rpm  to  obtain  maximum 
contact  forces  and  used  the  projectile  kinematics  (app  B) . 


4  Both  meshes  were  designed  with  the  help  of  computer  programs  1NV0L11  and 
GEARPAKAh2,  originally  shown  in  Progress  Report  No.  11  of  the  "Development  of 
Automated  Design  Optimization  Technique  for  Safety  and  Arming  Devices"  (Con¬ 
tract  No.  DAAK10-79-C-0251 ,  January  15,  1981).  Copies  of  this  report  may  be 
obtained  from  either  P.  R.  Tepper,  ARRADCOM  or  G.  G.  Lowen,  The  City  College  of 
New  York. 


The  following  shows  the  Input  requirements  of  the  program,  explains  the 
various  output  data,  and  discusses  the  manner  In  which  the  “ number -of-turna-to- 
arm"  la  obtained  for  a  given  spin  velocity. 


Input  Data 


The  first  portion  of  the  output  repeats  all  Input  data,  which  represent  the 
mechanism  parameters  of  the  M739  fuse.  Tnese  are  listed  both  as  computer  vari¬ 
ables  and  as  symbols,  according  to  the  various  appendixes  of  this  report  as  well 
as  reference  1. 


Escapement  Parameters 

A  -  a  ■  0.226  in.  (5.740  mm)  *  distance  between  pivots  0_  and  0o 

(fig.  2)  P 

B  -  b  -  0.168  in.  (4.267  mm)  •>  escape  wheel  radius 

C  *  c  ■  0.13138  in.  (3.337  mm)  »  pallet  radius  as  defined  by  figure  F-l 

of  appendix  F,  reference  1 

ALPHEN  *•  a  «  43.352  deg  -  entrance  working  surface  angle 
en 

ALPHEX  -  «ex  «  29.2981  deg  -  exit  working  surface  angle 

NT  *  4  -  number  of  escape  wheel  teeth  spanned  by  verge 

CONFIG  ■»  2  -  configuration  no.  2  (fuze  body  configuration  no.  2  in 
ref  1,  app  B) 

BREST  »  er  *  0  ■  coefficient  of  restitution 

LAMBDA  ■  X  *  92,93  deg  ■  angle  between  entrance  and  exit  pallet  radii 

(ref  1,  app  F,  fig.  P-1) 

N  *  22  *  number  of  escape  wheel  teeth 

For  details  of  the  above  nomenclature,  see  reference  1,  appendixes  C,  E, 

and  F. 


Mass  Parameters  of  Components 


Ml  *  m^  »  0.3165  x  10“^  lb-sec^/in.  (5.552  x  10*"^  kg)  ■  mass  of 
rotor 


b»2  m2  -  0.3275  x  10*5  lb-sec2/in.  (5.745  x  10“4  kg)  -  mass  of  gear 
and  pinion  no.  2 

M3  -  m3  ■  0.2631  x  lO"5  lb-sec2/in.  (4.615  x  10“4  kg)  -  mass  of 

escape  wheel  and  pinion  no.  3 

Me  -  m  -  0.1640  x  10-5  lb-sec2/ln.  ( 2 .87 7  x  JO"4  kg)  -  mass  of 

pallet 

1XX1  ■  1  »  0.1222  x  10“5  in.-lb-sec2  (1.383  x  10”?  leg  -  m2)  .  moment 

’F,l  of  inertia  of  rotor  with  respect  to  5^-axis  (through  cen¬ 
ter  of  mass,  see  fig.  A-3) 

1EE1  "I  «  0.1234  x  10"5  in.-lb-sec2  (1.397  x  10"/  ltg-m2)  -  moment  of 
nr,l  inertia  of  rotor  with  respect  to  r^-axis 

IZZ1  *  1  »  0.1967  x  lO”1*  in.-lb-sec2  (2.226  x  10”7  kg-m2)  -  moment  of 

??1  inertia  of  rotor  with  respect  to  pivot  axis  (*;,- axis) 

IXE1  -  I  -  -0.1012  x  10“6  in.-lb-sec2  (-1.145  x  10“8  kg-m2)  -  £  -n 
^nl  product  of  inertia  of  rotor 

1ZX1  «  I  -  -0.3656  x  10-'  in.-lb-sec2  (-4.137  x  10"y  kg-m2)  - 
*•*1  product  of  inertia  of  rotor 

IEZ1  -  I  -  -0.1770  x  10*7  in.-lb-sec2  (-2.003  x  lU“y  kg-m2)  -  n.-C, 
n^l  product  of  inertia  of  rotor 

1X2  -  lx2  -  0.2944  x  10” 7  in.-lb-sec2  (3.389  x  luJ'y  kg-m2)  -  moment  of 

inertia  of  gear  and  pinion  no.  2  (about  axis  normal  to 

pivot  axis) 

IY2  -  Iy2  -  0.2944  x  lo“7  in.-lb-sec2  (3.389  x  10~y  kg-m2)  -  moment  of 

Inertia  of  gear  and  pinion  no.  2  (about  axis  normal  to 

pivot  axis  and  perpendicular  to  x2-axis) 

IZ2  -  lz2  “  0.4026  x  10"7  in.-lb-sec2  (4.556  x  10“y  kg-m2)  -  moment  of 
Z  inertia  of  gear  and  pinion  no.  2  with  respect  to  pivot 
axis 

IXS  »  Ixs  -  0.2038  x  10“7  in.-lb-sec2  (2.307  x  lU~y  kg-m2)  “  moment  of 

inertia  of  escape  wheel  and  pinion  no.  3  (about  axis  nor¬ 

mal  to  pivot  axis) 

IYS  -  lys  -  0.2038  x  10~7  in.-lb-sec2  (2.307  x  10“y  kg-m2)  “  moment  of 

inertia  of  escape  wheel  and  pinion  no.  3  (about  axis  nor¬ 

mal  to  pivot  axis  and  perpendicular  to  xg-axis) 

IZS  “  1  -  0.2125  x  10“7  in.-lb-sec2  (2.405  x  10"y  kg-m2)  -  moment  of 

inertia  of  escape  wheel  and  pinion  no.  3  with  respect  to 
pivot  axis 


IXXP  -  I 


■  0.1721  x  10"tt  in. -lb-sec2  (1.948  x  10“  M  kg-rn2)  ■  moment 
^p  of  Inertia  of  pallet  with  respect  to  K  -axis  (through  cen¬ 
ter  of  mass,  see  fig.  A-2)  p 

IKKP  “I  *■  0.3038  x  10“a  in. -lb-sec2  (3.438  x  10“ 1°  kg-m2)  ■  moment 
nnp  of  inertia  of  pallet  with  respect  to  r^-axis 

IZZP  •  I  ■  0.1951  x  10"7  in. -lb-sec2  (2.028  x  10“*  kg-m2)  -  moment  of 
**p  inertia  of  pallet  with  respect  to  pivot  axis  (c^-axis) 

IXEP  "I-  *  0.  *  r  -n  product  of  inertia  of  pallet 

5np  'p 

IZXP  -  1  „  -  0.  ■  t  -E  product  of  inertia  of  pallet 

C£p  *P  P 

I1SZP  -I  ■  0.  ■  n  -r  product  of  inertia  of  pallet 
ntp  P  ^P 

General  Parameters 


RCl 

m 

rcl 

RCP 

m 

rcp 

KHOP 

m 

pp 

RPM 

m 

30,000 

PHI 1 CD 

m 

♦lc 

PSICCU 

m 

*c 

PHID 

m 

139  deg 

PHICUTU 

m 

14B5  deg 

0.0567  in.  (1.463  mm)  »  distance  from  pivot  of 
rotor  to  its  center  of  mass 

0  «  pallet  eccentricity 

0.0227  in.  (0.577  mm)  ■  pallet  pivot  radius 
spin  rate 

-120.134  deg  ■  rotor  angle  in  starting  position 
(fig.  2) 

0  deg  -  eccentricity  angle  of  pallet 

escape  wheel  starting  angle  of  initial  coupled 
motion 

cumulative  escape  wheel  angle  obtained  from  pro¬ 
duct  of  total  engaged  rotor  rotation  and  gear 
ratio.  The  total  rotor  rotation  for  the  M739 
fuze  is  46.41  deg,  while  the  gear  ratio  is  32. 
Thus,  PUICUTD  -  46.41  x  32  “  1485  deg. 


hU  “  u  -  0.10  -  coefficient  of  friction  of  gear  train 

(pivots  and  tooth-to-tooth  contacts)  and  escape 
wheel  pivot  (constant  for  a  computer  run) 

MU1  “  “  U.1U  «  coefficient  of  friction  of  pallet-escape 

wheel  interface  and  pallet  pivot  (constant  for  a 
computer  run) 


LU  ■  LL  ■  Ly  »  Ll  »  0.285  in.  (7.24  mu)  -  one-half  of  mean  distance 

between  bearing  plate  surfaces 

Gear  Parameters 

PSUbDl  “  P(Q  *  80  ■  diametral  pitch  of  mesh  no.  1  (rotor  and  pinion 

no .  2 ) 

PSUBD2  -  Pd2  -  100  ■  diametral  pitch  of  mesh  no.  2  (gear  no.  2  and 

escape  wheel  pinion) 

NG1  -  Nqj  -  64  “  number  of  teeth  of  rotor  (full  gear  no.  1) 

NG2  ■  Ny2  ■  36  “  number  of  teeth  of  gear  no.  2 

NP2  -  Np2  »  9  ■  number  of  teeth  of  pinion  no.  2 

NP3  ■  Npj  *  8  “  number  of  teeth  of  pinion  no.  3  (escape  wheel  pinion) 

CAPRP1  -  Rpj^  *  0.41214  in.  (10.466  mm)  -  pitch  radius  of  gear  no.  1 

(rotor) 

CAPRP2  ■  Kp2  “  0.19039  in.  (4.836  mm)  -  pitch  radius  of  gear  no.  2 
*^2  “  rp2  “  0.05796  in.  (1.472  mm)  -  pitch  radius  of  pinion  no.  2 

RP3  -  rp3  ■  0.04231  in.  (1.075  mm)  -  pitch  radius  of  pinion  no.  3 

(escape  wheel  pinion) 

THETA1  ■  3^  «  24.215  deg  -  pressure  angle  of  mesh  no.  1 

THETA2  ■  8^  "  27,326  deg  -  pressure  angle  of  mesh  no.  2 

R1  ■  -  0.250  in.  (6.350  mm)  ■  distance  of  rotor  pivot  from  spin 

axis 

R2  «  "  0.317  in.  (8.052  mm)  ■  distance  of  pivot  of  gear  and  pinion 

set  no.  2  from  spin  axis 

R3  -  *'3  ■  0.309  in.  (7.849  mm)  ■  distance  of  pivot  of  escape  wheel  from 

spin  axis 

R4  m  -  0.304  in.  (  7.722  mm)  -  distance  of  pivot  of  pallet  from  spin 

axis 

RHOl  "  Pj  "  0.03075  in.  (0.781  nun)  ■  pivot  radius  of  rotor 

KH02  “  P2  “  0.015  in.  (0.381  mm)  ■  pivot  radius  of  gear  and  pinion  no.  2 

RH03  -  ■  0.015  in.  (0.381  mm)  ■  pivot  radius  of  escape  wheel 


*.."1 


RH0F1  «  p„,  «  0*055  in.  (1*397  mm)  •  friction  thrust  radius  of  rotor 

*  (for  computation  of  friction  thrust 

radius  see  p  268  in  ref  8) 

RH0F2  ■  Pp^  ■  0.0294  in.  (0.747  mm)  -  friction  thrust  radius  of  gear  and 

pinion  no.  2 

RH0F3  ■  Pp^  •  0.0294  in.  (0.747  mm)  -  friction  thrust  radius  of  escape 

wheel  and  pinion  no .  3 

RHOF  -  p  -  0.1138  in.  (2.890  mm)  -  friction  thrust  radius  of  pallet 

r 

CAPRBl  ■  Rjji  ■  0.37588  in.  (9.547  ram)  ■  base  radius  of  gear  no.  1 

CAPRB2  -  Rjj2  11  0.16915  in*  (4.296  mm)  -  base  radius  of  gear  no.  2 

RB2  -  r^2  “  0.05286  in.  (1.343  ram)  -  base  radius  of  pinion  no.  2 

RB3  •  r^j  ■  0.03759  in.  (0.955  mm)  «  base  radius  of  escape  wheel  pinion 

CAPR01  "  Rqj  -  0.41425  in.  (10.522  mm)  -  outside  radius  of  gear  no.  1 

CAPR02  -  Rq2  “  0.19404  in.  (4.929  mm)  -  outside  radius  of  gear  no.  2 

R02  “  rQ2  -  0.07670  in.  (1.948  mm)  »  outside  radius  of  pinion  no.  2 

R03  -  rno  -  0.05580  in.  (1.417  mm)  ■  outside  radius  of  escape  wheel 

pinion 

J1  ■  Jj  -  0  *  Initialization  parameter  for  mesh  no.  1  [The  zero  value 

corresponds  to  earliest  possible  contact  of  mesh  (ref  3).J 

J2  ■  J2  ■  0  ■  initialization  parameter  for  mesh  no.  2 
Projectile  Parameters 

RX,  RY,  RZ  -coordinates  of  geometric  center  C  of  mechanism  plane  with 

respect  to  projectile  center  of  mass,  expressed  in 

projectile  fixed  coordinate  system  with  origin  at  point 
OpR  [This  system  is  parallel  to  mechanism  plane  fixed  XY 
system  (figs.  1  through  3  and  C— 1 )  . ] 

RX  -  Rjj  -  0.001  in.  (0.0254  mm) 

RY  -  Ry  «  0.001  in.  (0.0254  ran) 

RZ  *  R2  -  20.0  in.  (508  mm) 
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Projectile  Kinematics 


The  projectile  kinematics  are  programmed  in  subroutine  AERO  according  to 
the  expressions  given  in  appendix  B.  The  following  parameters  are  incorporated: 

For  equation  E-5,  KP  *  Kp  -  100 

For  equation  E-8,  THETIN  -  8  deg 
TV  •  2  deg 
KN  «  10 

DDZ  ■  X  -  -386.4  x  10  (corresponds  to  a  lu-g  decel¬ 
eration) 

Output  Data 

The  data  blocks  following  the  input  data  represent  the  results  of  various 
computations. 

Fuze  Geometry 

The  angles  BETA1D  -  to  BET  A3  D  -  83  and  GAMUA2D  -  y2  to  GAi-iMA4l)  -  y4 
are  printed  for  checking  purposes.  J  4 

Coupled  hot ion 

The  first  coupled  motion  output  refers  to  the  entrance  side  of  the 
verge.  For  each  time  T  of  the  coupled  motion,  the  following  variables  are 
computed : 

PHI  -  $  -  instantaneous  escape  wheel  angle  (deg) 

PHIDOT  ■  $  ■  escape  wheel  angular  velocity  (rad/sec) 

G  ■  g  -  pallet  -  escape  wheel  contact  position  (in.)  (ref  1,  app  C, 
eq  C-15) 

PS1D  ■  i(i  ■  pallet  angle  (deg) 

PSIDOT  -  •  pallet  angular  velocity  (rad/sec) 

PH  I  TOT  ■  •  cumulative,  escape  vheel  angle  (deg) 

F23  ■  f'23  -  normal  contact  force  of  gear  no.  2  on  pinion  no.  3  (lb) 


F12  *  Fj,  2  ■  normal  contact  force  of  gear  no.  1  on  pinion  no*  2  (lb) 

PN  -  Pn  ■  normal  contact  force  between  escape  wheel  and  pallet  (lb), 
computed  according  to  equation  D-529  in  appendix  D 

PNPSI  •  Pn  ■  normal  contact  force  between  escape  wheel  and  pallet 

(lb),  computed  according  to  equation  D-527  in  appendix  D 
(serves  for  checking) 

••  * 

DPHI2  -  4  ■  escape  wheel  angular  acceleration  (rad/sec4) ,  Runge-Kutta 
output 


Free  Motion 


The  first  free  motion  on  the  exit  side  follows  the  coupled  motion  on  the 
entrance  side  of  the  verge.  For  each  time  T  of  the  free  motion,  the  following 
variables  are  evaluated; 

PHI  -  instantaneous  escape  wheel  angle  (deg) 

PHIDOT  “  £  ■  escape  wheel  angular  velocity  (rad/sec) 

PSI  -  ■  pallet  angle  (deg) 

PSIDOT  -  ijt  »  pallet  angular  velocity  (rad/ sec) 

FF12  -  Fpil  "  norraal  contact  force  of  gear  no.  1  on  pinion  no.  2  for 
free  motion  (lb) 

FF23  ■  *f23  “  normal  contact  force  of  gear  no.  2  on  escape  wheel  pin¬ 
ion  for  free  motion  (lb) 


Impact 


The  first  exit  impact  follows  the  first  exit  free  motion.  Just  preced¬ 
ing  the  IMPACT  label,  the  program  prints  the  values  of  VP  ■  V™^  and  VS  »  Vgjjj 
which  stand  for  the  pre-impact  velocity  components,  normal  to  the  verge  face  of 
both  the  pallet  and  escape  wheel  contact  points  (ref  1,  app  D,  eq  D-13).  Subse¬ 
quent  to  the  IMPACT  label,  the  following  variables  are  evaluated: 

PHI  *  <j>  ■  Instantaneous  escape  wheel  angle  (deg),  same  as  before 
impact 

PHIDOT  •  |  »  post-impact  escape  wheel  angular  velocity  (rad/sec) 

PSI  -  \|>  »  pallet  angle  (deg),  same  as  before  impact 
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PSIDOT  ■  I  ■  post-impact  pallet  angular  velocity  (rad/sec) 

PHITOT  ■  A  ■  cumulative  escape  wheel  angle  (deg),  same  as  before 
impact 

VP  *  »  post-impact  normal  velocity  component  of  pallet  at  con- 

A  tact  point  (ref  1,  eq  D-15) 

VS  *  Vg^£  •  post-impact  normal  velocity  component  of  escape  wheel 
tooth  at  contact  point  (ref  1,  eq  D-13) 

In  the  present  program,  the  post-impact  VP  is  equal  to  VS  since  the  coefficient 
of  restitution  is  zero. 


Number  of  Turns-to-Arm  and  Maximum  Contact  Forces 


The  number  of  turns-to-arm  at  30,000  rpm  is  obtained  with  the  help  of 
that  time  Tj^j  which  corresponds  to  the  escape  wheel  angle  PHICUTD  ■  1485  deg. 
Thus,  with  l'i485  "  0.05094  sec, 

number  of  turns-to-arm  *  '"Qyffi  x  0 .05094  -  25 .47 

OU 

The  maximum  non-impact  contact  forces  for  the  total  cycle,  for  both 
coupled  and  free  motion,  are  listed  at  the  end  of  the  output. 


CONCLUSIONS 


While  it  was  not  the  purpose  of  this  investigation  to  undertake  a  parametric 
study  of  the  mechanism  for  which  the  program  was  written,  the  program  was  suffi¬ 
ciently  tested  to  confirm  that  such  a  study  is  possible.  It  may  include  varia¬ 
tions  In  masses  and  moments  of  inertia  of  all  components;  variations  in  the  loca¬ 
tions  of  the  centers  of  mass  of  the  verge  and  the  rotor;  variations  of  gear, 
escapement,  and  fuze  geometries;  as  well  as  various  friction  and  coefficient  of 
restitution  conditions.  In  addition,  the  aeroballistic  data  can  also  be  varied. 
This  makes  it  possible  to  determine  the  functioning  limits  of  the  mechanism  under 
pathological  projectile  flight  conditions. 

The  present  work  reports  only  on  a  single  test  run  using  the  M739  fuze  S&A 
data  with  a  system  coefficient  of  friction  of  0.1.  This  is  assumed  to  be  repre¬ 
sentative  of  actual  test  conditions  since  previous  simulations  of  pin  pallet 
escapements  showed  that  the  range  of  actual  experimental  results  (with  spin  only) 
may  be  reproduced  with  coefficients  of  friction  between  0.1  and  0.2,  and  a  spe¬ 
cial  lubricant  is  used  in  conjunction  with  the  M739  fuze.  This  choice  of  coeffi¬ 
cient  of  friction  is  proven  by  the  good  agreement  with  experimental  results.  A 
zero  coefficient  of  restitution  is  used  in  the  impact  model  (ref  1  and  2). 

Previous  high-speed  motion  picture  observations  of  pin  pallet  escapements 
showed  that  the  impacts  were  essentially  inelastic  and  that  a  zero  coefficient  of 
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restitution  was  justified.  Similar  observations  made  on  the  detached  lever 
escapement  of  the  M577  fuse  timer  confirmed  this. 

The  test  run  showed  that  for  a  spin  rate  of  30,000  rpm,  together  with  small 
precession  and  nutation  velocities  chosen  in  the  manner  shown  in  appendix  E,  the 
number  of  turns-to-arm  is  essentially  the  same  as  that  obtained  in  reference  6 
where  only  spin  was  considered. 
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ANGULAR  VELOCITIES  AND  ACCELERATIONS  IN  TERMS  OK  KROJECT1LE-EIXED  COORDINATES 


A  projectile  which  experiences  general  aeroballistic  motion,  i.s.  spin  about 
an  axis  through  its  center  of  mass  as  well  as  precession  and  nutation  of  this 
spin  axis  with  respect  to  its  center  of  mass,  is  shown  in  figure  A-l.  (In  the 
figure,  the  spin  axis  coincides  with  the  geometric  axis.) 

The  spin  angle,  spin  velocity,  and  spin  acceleration  are  expressed  by  the 

•  •• 

time  dependent  quantities  $  ,  $  ,  and  (The  subscript  E  stands  for  the  Euler 

angles,  which  are  Involved  in  this  derivation.)  Similarly,  the  kinematic  quanti¬ 
ties  associated  with  the  precession  are  vp $  ,  and  <j^,.  The  nutation  variables 
..  b  E  h 

are  0  ,  9  ,  and  0„  (refs  7  and  8). 

b  tt  zt 

With  spin,  precession,  and  nutation  angular  velocity  vectors,  together  with 
their  associated  angles  (fig.  A-l),  orthogonal  angular  velocity  components  in 
terms  of  the  projectile  fixed  x-y-z  system  may  be  obtained  as  follows: 

Let 

%/a  “  V  +  “yj  +  “s*  (A”1) 

where  w  ,  represents  the  angular  velocity  of  the  projectile  b  with  respect  to 
the  inerc&l  frame  a.  Then 

Ux  "  ®EC°8  ^E  +  ^E  8ln  9E  8in  ^E  <A-2) 

<*»y  -  -  3g  sin  <|>E  +  'I»E  ain  0g  cos  <1>E  (A-3) 

%  "  +  C°8  0E  (A-4) 

The  absolute  angular  acceleration  of  the  projectile,  i.e., 


U).,  “Ui+Mjf  +  wit 
b/  a  x  r  z 


( A-  3a ) 


is  obtained  by  differention  of  the  body  fixed  quantities  with  respect  to  time. 
Thus 


“x  “  0ECO8  »E  “  ^E^E  8in  ^E  +  *E  8in  ®E  8in  ^E  (A-5b) 

+  i}»E0E  cos  0^  sin  ifE  +  sin  9g  cos  ^ 

»y  ■  -  fljj  sin  <j>E  -  8^  cos  +  iE  sin  6g  cos  $g  (A-3c) 

+  tE0E  cos  0E  cos  <J>E  -  ^E  sin  3g  sin  4>g 
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FIGURE  A*1 


Xo 

PROJECTILE-FIXED  x-y-z  SYSTEM 

Note i  This  system  is  luter  described 
with  capital  letters  X-Y-Z 

NOTE:  This  system  is  later  described  with  capital  letters  X-Y-Z. 

Figure  A-l.  Projectile-fixed  x-y-z  system 
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-  *E  +  *B  co*  eg  - 


E  ,in  0E 


(A-5d) 


ANGULAR  VELOCITIES  AND  ACCELERATIONS  IN  TEKhS  OK  PROJECTILE- FIXED  COORDINATES 


Pallet-Fixed  Coordinates 


The  relationship  of  the  pallet-fixed  £  -  n  -  ?  system  with  respect  to  the 

projectile  fixed  X-Y-Z  and  x'-y'-s*  systems  %s  stJown  iSi  figure  A-2  (ref  1). 

The  £  -  n  plane  is  parallel  to  the  x'-y'  and  X-Y  planes  and  contains  the 
pallet  centPer  of  mass  Cp.  The  £p-axis  is  parallel  to  the  z  and  s'  axes. 

The  pallet  angles  tp  and  tp  are  measured  in  the  £  -  n  plane  and  are  other¬ 
wise  defined  as  in  reference  1?  before  determining  tt&  absolute  angular  velocity 
and  acceleration  of  the  pallet,  a  number  of  unit  vectors  should  be  defined. 
According  to  equations  b-2E  and  B-2V  of  ref  1: 


i'  »  -  cos  flji  -  sin  ft3j 


(A-b) 


j '  -  sin  g..i  -  cos  0..j 


(A- 7  ) 


Further,  when  expressed  in  the  primed  system,  pallet  fixed  unit  vectors  become: 


n-  *  cos  pi'  +  sin  Bj' 

S 

n  -  -  sin  gl'  +  cos  gj‘ 


(A-8) 


(A-y) 


where 


$  ■  ip  +  ip 


( A-lU) 


If  equations  A-6  and  A-7  are  substituted  into  the  above  expressions,  after 
some  trigonometric  simplifications,  the  following  expressions  are  obtained  for 
the  pallet  fixed  unit  vectors  in  terms  of  the  X-Y-Z  system: 


Hr  -  -  cos  a’l  -  Bin  a'J 


(A-ll) 


P3 


Figure  A-2,  Pallet-fixed  £  -n  -c  coordinate  system 


5WK 


^  ^ rv*  ^77  ■u’’  ''7*  *  r  v;-  vr 


r^'.T'.v.Tr'V'.^  'r^;r  ■..•T.’i"*  :*jcA*.ra 


and 


n  •  ain  o'i  -  cos  a'j 


(A-12) 


where 


«'  -  *  +  *  +  B 

c  j 


(A-13) 


Because  of  the  given  parallel  axes, 


n  ■»  k*  •  k 

S 


(A-14) 


It  the  relative  angular  velocity  of  the  pallet  P  with  respect  to  the  projec¬ 
tile  b  is  given  by 


Vb  ■  *  5C(>- 


(A-15) 


then  its  absolute  angular  velocity  Is  given  by 


*p/»  ‘  Vb +  “b/p 


(A-16) 


To  express  equation  A-16  in  pallet-fixed  terms,  it  is  necessary  to  transform 


equation  A-l  which  gave  (o^a*  According  to  equations  A-U  and  A-12: 


i  •  -  cos  ex '  n„  +  sin  a'  n 

€  n 

P  P 


(A-l/) 


and 


j  *  -  sin  a'  n_  -  cos  a'  n 

S  np 


(A- Id) 


Thus,  one  obtains  in  the  pallet  fixed  system: 


(ii 


b/a  “  "f^x  c08  tt'  +  ,0y  sin  °'  1  nF,  +  fwx  8in  °* 


o)y  cos  a*  ]  nn 


+  n 
2 


(A-iy) 


K 

g 

f. 


\ 

k. 


c 


tv 


u 


i 

C.'1 
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Finally,  the  absolute  angular  velocity  of  the  pallet  io  .  becomes  with  equations 
A-15  and  A-16:  p/a 


b)  ,  ■  <o_  n_  +  w  n  +  w.  n_ 

p/a  t  £  r\n  X  X  X 
P  P  P  P  P  P 


(A-20) 


where 


w_  »  -  (uj  cos  o’  +  ui  sin  a1) 

£p  *  y 


w  ■>  (o  sin  ct'  -  u»  cos  a' 

np  x  y 


w  -  0)  +  \}» 

CP  2 


(A-21) 
( A-22) 
(A-23) 


The  absolute  angular  acceleration  io  ,  of  the  pallet  is  obtained  by  the 
differentiation  with  respect  to  time  of  The  measure  numbers  of  equation  A-20, 
i  •  e  •  • 


u>  ,  r_  +  u)  n  n„ 

p/a  *■  K  n  n  C  C 
P  P  p  P  P  P 


(A-24) 


where 


w-  -  -  fw  cos  a'  -  id  t}»  sin  a*  +  a>  sin  «'  +  u»  i  cos  a*  ) 
5vx  x  y  y  J 


(A-25) 


•  •  •  .  •  • 

u>  *01  sin  a'  +  w  cos  <v  -  u>  cos  o'  +  u  tic  sin  a' 
Tip  x  x  y  y 

*  • 

0)  *  CO  +  tl) 

s  z 


( A-26) 
( A-27) 


Rotor-Fixed  Coordinates 


The  relationship  between  the  rotor-fixed  C,.  -  n,  -  C.  system  and  the  projec¬ 
tile  fixed  X-Y-Z  system  is  shown  in  figure  A-3  (see  also  ref  1),  The  £.  -  n, 
plane  is  parallel  to  the  X-Y  plane  and  contains  the  center  of  mass  Cj  of  the 
rotor  (referred  to  as  link  1  below) .  The  ^-axis  connects  the  point  Oj  on  the 
rotor  pivot  centerline  and  point  . 


•anti  t 

x-axis. 


Tha  rotor  angles  $  and  ^  ara  measured  In  tha  5  -  n,  plana.  rapra- 

tha  Initial  posltioft  of  thi  5^-axla,  l.a.  tha  ClQ  axis/ with  r aspect  to  the 


The  unit  vectors  associated  with  tha  rotor-fixed  system  ara  given  by: 


n^  -  cos  ($lc  +  <^>1  +  sin  (<frlc  +  <^>3 


(A-2B) 


-  -  »in  <$lc  +  ♦j)!  +  cos  (*lc  +  ♦1)J 


(A-29) 


n  »  Ic 

Because  of  the  gear  train,  it  is  best  to  let  (see  ref  1,  eq  B-123): 
*1  "  N31  *T 


(A-30) 


(A-31) 


Y  "  ♦ic  +  N31  S 


(A-32) 


The  absolute  angular  velocity  u  .  of  the  projectile  is  expressed  in  terms  of  the 
rotor  fixed  coordinates  with  the't&lp  of  equations  A-l  and  A-28  to  A-30: 


“b/a  -  “a1  +  V  +  “*k 


-  o>  ( cos  "Y  n-  -  sin  Y  n  ) 

x  _ni 

+  a)  (sin  Y  n-  +  cos  y  n  )  +  u>  n 
y  ^  nl  2  ;1 


“b/a  “  K  C08  Y  +  Uy  8in  YJ  \ 


+  [-  w  sin  y  +  uj  cos  y1  n 
1  x  y  J  n. 


+  u>  n 


(A-33) 


To  obtain  the  toj:al  angular  velocity  of  the  rotor,  its  relative  angular 

velocity  must  be  added  vectorially  to  equation  A-33: 
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f 5T^  'M 'V  If ?»,7k: *."<  ’VJ*  t"¥  mn  vj  v 


S(E) 


Rotor-fixed  ^1-n1-C1  coordinate  system 
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5l/»  "  V.  +  N3l*  \ 
Then,  with  equation  A- 3 3 : 


(A-34) 


where 


W,  /.  ■  wr  nr  +  »n  n.  +  w  n 
u*  h  *1  "l  "1  '1  'l 


(A-35) 


wr  ■  (0  cos  y  +  w„  sin  y 
%,  x  y 


w  *  -  to  sin  y  +  id  cos  y 
n.  x  y 


V  “  “a  +  N31  ♦ 


(A-36) 

(A-37) 

(A-38) 


To  obtain  the  absolute  angular  acceleration  id.,  of  the  rotor,  differentiate  the 
measure  numbers  of  equation  A-35  with  respect  to  time.  Therefore 


where 


1/a 


uL  n,  +  w  n_  +  u>  n 

f*i  f*i  ni  ni  h  h 


W, 


ID 


W.. 


U)  cos  Y  -  «*  N_,  $  sin  Y  +  w  aln  y  +  u>  N,.  d>  cos  y 
x  x  31  y  y  31 

-  w  sin  y  -  <11  N„.  $  cos  y  +  10  cos  y  -  m  N  $  sin  y 
x  x  ji  y  y  j! 

+  N31  ♦ 


( A-39) 

(A-40) 
( A-4 1 ) 
( A-42) 


Two  Ways  of  Obtaining  Expressions  for  Absolute  Angular  Velocities  and  Accelera¬ 
tions  of  Components,  Such  as  the  Pallet,  in  Terms  of  Projectile-Fixed  Coordinates 


Let  it  be  required  to  express  the  absolute  angular  velocity  and  the  absolute 
angular  acceleration  of  the  pallet,  as  given  by  equations  A-20  and  A-24,  in  terms 
of  the  projectile  fixed  x’-y’-z’  system  (fig.  A-2).  This  is  accomplished  by 
substitution  of  equations  A-8,  A-9,  and  A-14  for  the  unit  vectors  of  the  above 
expressions . 
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v;’*.- rr  tTTTT.TT?  ..*"T'v -.t  v  *  >'  *-rrr*-. x-.-r-,-*-.  »-;< 


s  ■ns’rs 


This  leads  to 


Wp/a  -  “(<t>x  COS  e3  +  wy  sin  P3)i'  +  (wx  sin  B3  -  wy  cos  0^) J* 

x  y 1  s' 


+  (u)^  +  $i)k' 


(A-43) 


and 


“p/a  "  l"“x  co8  ®3  ”  “y  8in  +  ♦  (wx  8in  P3  “  “y  008  Bj)]*' 


x'y's 


+  [«x  sin  $3  -  ui  cos  Bj  +  I  (ux  cos  B3  +  u  sin  ($3)]j' 


+  [it»z  +  t)»]k' 


(A-44) 


The  same  results  may  be  obtained  if  the  vector  i|i  is  interpreted  aa  a  vari¬ 
able  vector  in  the  primed  system,  which  is  attached  to  the  projectile,  first  let 
the  absolute  angular  velocity  of  the  projectile  be  expressed  in  terms,  of  the  x'- 
y'-r'  system.  This  is  accomplished  by  substituting  equations  A-6  and  A-7  for  the 
unit  vectors  I  and  3  in  equation  A-l.  (The  unit  vector  E  is  also  E'  in  the  pre¬ 
sent  systems.) 


Then 


-  -(u>x  cos  B3  +  w  sin  33)i* 


b/a  , 
x  y  z 


+  (wx  sin  f?3  -  a)  cos  P3)J'  +  u^k' 


y  '  3 z  ( A-45) 

The  absolute  angular  velocity  of  the  pallet  becomes,  according  to  equation  A-16 

a  .  m  $  k 1  +  w,  ,  (A-46a) 

"/a,'y'*'  b/ax','*' 

or 


u_/_  •  -(w  cos  3,  +  at  sin  B..)i'  +  (w  sin  -  a>  cos  P..)j"' 

x  s  y  j  x  j  y  j 


P  t  •  vi  9  ■ 
x  ;  2 


+  (u>  +  ^)k* 


(A-46b) 


T: 


t 


r 

i» 


w 

r. 


This  is  identical  to  equation  A-43. 

To  obtain  the  absolute  angular  acceleration  of  the  pallet,  the  following 
expression  in  which  $  E'  is  treated  as  a  variable  vector  in  the  x'-y'-z1  system 
must  be  evaluated.  Therefore, 
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(A-47) 
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APPENDIX  B 

ANGULAR  MOMENTUM  AND  ITS  DERIVATIVES  IN 


VARIOUS  COORDINATE  SYSTEMS 
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Angular  Momentum 


M 

l;| 

1% 

r 


:v 

is 


<  , 

i 


The  angular  momentum  vector  Hq,  with  respect  to  a  point  0.,  has  the  general 
expression: 


m 

H„  - 

fl  U>  -  I  U  -  I  111  ll 

0 

1  XX  X 

xy  y 

zx  z- 

+ 

[-1  Ui 

+  I  U) 

-  I  u  ]j 

1  xy  x 

yy  y 

yz  zJJ 

[-1  m 

-  i  u 

+  I  m  lk 

L  XX  X 

ya  y 

zz  ZJ 

The  above  holds  for  all  types  of  body-fixed  and  space-fixed  coordinate  systems. 
If  principal  axes  are  involved,  the  products  of  inertia  vanish.  Note  that  the 
angular  velocity  component  must  be  absolute. 

Derivative  of  Body-Fixed  Angular  Momentum  Vector:  Torque  Equation 


Body  b  in  general  motion  la  shown  in  figure  B-l.  It  contains  the  body  fixed 
X-Y-Z  system  and  its  angular  momentum  may  be  expressed  with  the  help  of  equation 
B-l. 


The  time  derivative  of  the  angular  momentum  with  respect  to  the  inertial  X0~ 
Y0-Z0  system  is  obtained  from: 


where 


**0/X  Y  2  dt  ^O^XYZ  +  “  X 


U-2) 


o  o  o 


dt  ^O^XYZ 


■  derivative  of  the  measure  numbers  in  equation  B-l 


and 

S  x  H0  -  (u»xi  +  uiyj  +  u>zk)  x  HQ 

It  is  to  be  recalled  at  this  point,  that  the  absolute  angular  acceleration  of 
body  b  is  given  by: 

«  ■  «  1  +  u J  +  J  k  (B-3) 

X  y  Z 

Both  3  and  (!)  are  now  expressed  In  terms  of  the  body  fixed  coordinates. 


'•i’V'CTs.T 
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Upon  performing  all  oparationa  of  equation  B-2,  the  torque  equation  with 
respect  to  point  0  results: 


%  "  “o/X  Y  Z  “ 

0  0  0 


flu»+wu)fl“I  1+1  fu  to  -  (I)  1 
1  XX  x  y  z'  zz  yy;  xy  z  x  yJ 

“1  fo)  +  (j  b  ]  •  I  fa)  -  a)  ^lli 
zx  ^  z  x  yJ  yz  K  y  z 

+  [l  id  +  u  in  (i  -1  )  +  I  (ut  0)  -  0)  1 

1  yy  y  :  xx  zzJ  yz  K  x  y  zJ 

-1  (a)  +  w  w  )  -  I  f  0)^  -  0)  ^llj 
xyv  x  y  s-1  zx  v  z  x  ;JJ 


+  fl  (a  +  a)  a)  fl  “I  1  +  I  faio)  —  ail 
L  zz  z  x  y  *■  yy  xxJ  zx  1  y  z  x' 

-I  (a)  +  a)a)l-I  fai^-o)  ‘hlk 

yz  v  y  x  z ;  xy  ^  x  y  ' J 


(B-4) 


When  Ixy  -  Iyz  -  Izx  ■  0,  the  above  expression  becomes  the  well  known  Euler 
e  equation. 


torque  equat 


Derivative  of  Angular  Momentum  Vector  Which  is  Described  in  Terms  of  the  Body- 
t’ixed  System  of  a  Carrier:  Vector-Torque  Equation 

The  carrier  body  b  which  has  general  rotational  motion  is  shown  in  figure  B- 
2.  Its  absolute  angular  velocity  and  acceleration  are  given  in  terms  of  the 
indicated  body-fixed  system,  i.e. 


“b/a  "  “b/a*Z  +  “b/ayj  +  “b/az* 


(B-3) 


Ub/a  "  Ub/ax*  +  ^b/ay^  +  ^b/az^ 


(r-6) 


respectively. 

The  symmetrical  body  c  rotates  about  an  axis  parallel  to  the  Z-axis  with 
respect  to  body  b  with  the  relative  angular  velocity 


Sc/b  ■  “c/b(t)C 


(B-7) 
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If  th«  absolute  angular  velocity  of  the  body  is  expressed  in  the  X-Y-Z  system  and 
fixed  to  the  carrier  (projectile),  the  following  is  obtained: 


“c/a  "  “c/b  +  %/a 


(B-8) 


(For  comparison  see  equation  A-46.) 


To  obtain  the  absolute  angular  velocity  i  ,  in  terms  of  the  projectile 
fixed  system,  Interpret  as  a  variable  vector  in  the  X-Y-Z  system.  Then, 


where 


®X  Y  Z 
o  o  o 


“c/b  +  Vc  x  “c/b  +  V. 


(B-9) 


to  .  -  do  .  k,  the  relative  angular  acceleration  of  component  C  with 

c  *  respect  to  projectile  b 

-  given  by  equation  B-6 

Because  body  c  is  symmetrical,  its  products  of  inertia  with  respect  to  Its 
center  of  mass  Cc  are  zero.  This  symmetry  also  makes  it  possible  to  express  its 
angular  momentum  with  respect  to  point  Cc  in  terms  of  the  body-fixed  system  of 
the  carrier  b.  (Regardless  of  the  angle  of  body  c  with  respect  to  body  b,  the 
moments  of  Inertia  Ixx  and  lyy,  expressed  in  terms  of  body  b,  remain  invariant.) 
The  angular  momentum  vector,  with  respect  to  point  Cc,  appropriately  reduced, 
hecomes  according  to  equations  B-l  and  B-8: 


H„  -  I 


;%/ax*  +  *yy%/ay^  +  *zz^%/az  +  “c/b  ^ 


(B-10) 


The  vector  H-,  must  be  interpreted  as  a  variable  vector  in  the  carrier-fixed 
c 

coordinate  system.  Its  time  derivative  with  respect  to  the  inertial  system  is 
therefore  obtained  by  the  following  operations: 


'Cc/X  Y  z  -  *T  <*=>.  +  V.  *  acc 

OOO 


(B-l 1 ) 
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When  applied  co  equation  B-10,  the  following  is  obtained: 


‘xx  "b/ax1  +  lyy 


w 


(“k/..  + 


+  (“b/ax*  +  “b/ay^"  +  “b/az“)  *  lIxx“b/»x1'  +  *yy  “b/ayJ 
+  ^zz (“b/az  f  “c/b^( 

(B-12) 


The  above  becomes  the  torque  equation  with  respect  to  point  Cc: 

H_  “  h„  ■  [l  .  +  I  w.  ,  ftiJ.  ,  +<*>,.)-!  u),  2  w.  .  ji 

Cc/X  y  z  Cc  1  xx  b'ax  zz  v  Waz  c/h;  yy  *>/ay  VazJ 
o  o  o 

+  ^yyub/ay  +  ^xx^b/ ax“b/ 22  “  *zz“b/ax  (“b/az  +  “c/b^ 

+  rzz  (“b/az  +  “c/b)E 


APPENDIX  C 

ABSOLUTE  ACCELERATION  OF  GEOMETRIC  CENTER  C 


OF  THE  S&A  PLANE 
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The  relationship  between  the  center  of  mass  CpR  of  the  projectile  end  the 
center  of  the  plena,  where  the  S&A  mechanism  la  loceted,  la  shown  In  figure  C- 
1.  Note  that  the  origin  of  the  projectile  fixed  X-Y-Z  ay  a  tern  Ilea  on  the  geo¬ 
metric  axle  of  the  projectile*  This  axis  la  aaeumed  to  be  parallel  to  the  spin 
axis  of  the  projectile.  The  canter  of  nasa  of  the  projectile,  about  which  all 
rotation  takes  place,  Ilea  In  the  name  plana  as  the  origin  0pR  of  the  body  fixed 
system*  The  position  vector  of  point  C  with  reapect  to  the  center  of  mass  la 
given  by: 


R-Rl+Rj+Rk 
x  y  a 


(C-l) 


It  la  assumed  that  the  deceleration  of  the  center  of  mass  due  to  drag  is  only  In 
the  Z-dtrection  and  that  it  la  given  by 


ACpR/ground  Z  k 


(C-2) 


The  absolute  acceleration  ^c/Rroun<i  of  point  C,  may  then  be  obtained  from: 


Wound  ’  *0,,/  ground  +  «(•  **>+««* 


(C-3) 


where- w  and  i  are  obtained  from  equations  A-l  and  A-5,  respectively.  (For  clar¬ 
ity  they  were  designated  as  and  u^a  in  appendix  A.) 

When  the  operations  of  equation  C-3  are  carried  out  and  equation  C-2  is 
substituted,  the  following  is  obtained: 


where 


Ao/  -Gi+Gj+Gk 

C/ground  x  yJ  z 


(C-4) 


(yy  +  - 

KRx  +  “xRx  >  “y  ' 


(<D  2  +  w  2  )  R  +  (u>  R 
v  y  z  J  x  '  y  z 


f id  +0)  )  R  +  (w  R 

v  x  z  1  y  v  z  x 


MR  ) 
z  y ' 


»R  ) 
x  z  * 


(C-5) 


(C-6) 


(uR  +  m  R  ]  0)  - 
v  x  x  yy  z 


fu>  +  w  2  )  R  +  f i  R  -  uR  )  +  Z 
'  x  y  '  z  v  x  y  y  x ' 


(C-7) 
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APPENDIX  D 

DYNAMICS  OF  ROTOR-DRIVEN  S&A  MECHANISM  WITH  A  TWO-PASS  INVOLUTE  GEAR  TRAIN 
AND  A  VERGE  RUNAWAY  ESCAPEMENT  OPERATING  IN  AN  AEROBALLISTIC  ENVIRONMENT 
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Geometry  of  fuse  Body  Configuration* 


Tha  too  poaslblo  fuse  body  configuration*  of  rafaranc*  i  ara  alao  accoramo- 
datad  in  th*  praaant  simulation*  Tharafora,  all  vorlc  concarnlng  fus*  body  anglaa 
in  this  aarliar  report  ara  also  applicable  hare. 


Dynamics  of  Pallet  and  Escape  Wheal  In  Coupled  Motion 


Absout*  Acceleration  of  Pallet  Pivot  0, 


Tha  position  of  th*  pallet  pivot  0P  with  respect  to  the  geometric  center 
C  of  tha  mechanism  plana  Is  shown  In  figure  D-l .  In  addition,  tha  relationship 
of  tha  projectile  fixed  x'-y'-s'  system  to  tha  projectile  fixed  X-Y-Z  system  is 
indicated. 

Th*  absolute  acceleration  of  point  0p  is  given  by: 


A0  /ground  "  ^  /C  +  AC/ground 
P  P 

where,  ^j/ground  iB  *iveQ  by  equation  C-4  of  appendix  C  and 


CD-I) 


A0  /c  "  “  x  Cu  x  <^)  +  «  x 


(D-2) 


In  the  above,  u>  and  w  are  obtained  from  equations  A-l  and  A-5,  respectively,  and 


n4 


(D-3) 


where,  in  the  primed  coordinate  system 


n4  -  cos  Ypi'  +  sin  Ypj ' 


(D-4) 


After  transformation  into  the  X-Y-Z  system  with  the  help  of  equations  A-6  and  A-7 
and  some  trigonometric  rearrangement,  the  following  is  obtained: 


n4  -  -  cos  (y*  +  03)i  -  sin  (y‘  +  03)j 


’*  <“•'  Km  «r*  ■*-, 


(D-3) 


r»  rr.  viiTvwr  TiTr^Tr^r.  ”■%  i'.r. 


rsvw.'^rw.'TOTme 


Equation  D-3  My  now  bo  written  as: 


&4  *  +<?yJ  "  "*4  cot  (Yi  +  “*4  Bin  (Yp  +  *3® 


(D-6) 


With  tho  above  ,  equation  D-2  la  now  evaluated: 


where 


*0  /C  "  Hxl  +  V  +  *ik 
P 

a  o  j 

H  *  ffl/  u  ti)  —  j?.  U  +  w  1  -  <R.  u  1 

X  l^4y  x  y  »‘4X  V  y  z  >  °%y  ZJ 


<D-/) 


(D-8) 


[tft,  b)  M  *•  4t.  fb)^  +  w^)+(R,  0)  ] 

1  4x  xy  4y  v  x  t  ‘  4x  z J 


(D-9) 


[  (<?,  u>  +  ft.  Ill  ]  III  +  (fl,  u  -  ft .  u  )  1 

1  v  4x  x  4y  y ;  a  v  4y  x  4x  y  ' J 


(D-lU) 


The  acceleration  A^  /groun^  i8  evaluated  according  to  equation  D-l  with  the  help 
of  equations  C-4  and  l>-7,  i.e.t 


\  /ground  *  +  +  (°y  +  V5  +  K  +  “J1 

P 


(0-11) 


for  later  computational  convenience,  the  above  expression  is  transformed  into  the 
x'-y’-z'  system: 


where 


A..  <  •  Ki1  +Kj'  Hk1 

Op/ ground  x  yJ  z 


Kx  “  "  (Gx  +  Hx)  008  03  ■  ^Gy  +  rty)  8itl  S3 


(0-12) 


(0-13) 


K.  -  (g  +  H  )  sin  &  -  (g  +  H  )  cos  S. 

y  xx'  3  s  y  y  3 


(0-14) 


Kz  -  Uz  +  Hs 


(0-15) 


Acceleration  of  Pallet  Center  of  Mass  Cp  with  Respect  to  the  Pallet  Pivot  0p 


When  the  relative  acceleration  of  the  pallet  center  of  mass  with  respect 
to  the  pallet  pivot  is  formulated  in  terms  of  the  pallet-fixed  K  -  n  -  C  coor¬ 
dinate  system,  the  following  is  obtained  (fig.  D-2):  P  P  P 


51 


where 


V°p  ■  V.  *  <v.  *  'cp^  +  "p/a  *  rcp”Cp 


(U-16) 


w_/a  “  “r  nr  +  n«  +  Wr  n». 

P/«  5_  n_  n  C  5_ 

P  P  P  P  P  P 

mm 

“p/a  •  \  \  +  %  \  +  \  \ 

P  P  P  P  P  P 

r  »  0  -  C 
cp  p  p 

Appropriate  substitution  and  evaluation  of  equation  D-16  furnishes: 

/0  “  rc  It"  (“x  8itl  *'  ”  “y  cos  “')2  -  («z  +  *)*]  n 

P  P  P  ^p 


+  [sin  o'  cos  o'  (to  2  -  u i  2)  +  uu  (coe2a*  -  sin2a' ) 

/  x  x  y  7 

+  +  ♦] 


+  [-  sin  o'  (a>x  +  2aiy  tp  +  u)yWz)  +  cos  a'  (ai 


-  2(0.  *  -  uyoj]  ^  } 


(U-17) 


The  above  expression  is  now  transformed  into  the  x'-y’-z'  system  (again 
tor  later  computational  convenience)  with  the  help  of  equations  A-b,  A-y,  and  A- 
14: 


where 


\  10  •  \  1‘  *  Ty  5'  +  T*  E' 

P  P 


2  2 

T  “  r  [-W  sin  6-  sin  o'  -  w  cos  0,4  cos  o' 
x  cp  x  j  y  j 


(U-lba) 


+  oj i  to  sin  (o'  +  0  )  -  (u)  +  ?»)2  cos  0  -  (u  +  (p)  sin  0]  (D-lbb) 
x  y  j  z  z 
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* 

% 


r  2  2 

T  -  r  I  -  u>  cos  (L  sin  a*  +  u>  sin  jL  cos  a1 

v  or*  *■  v  ^  v  ^ 


y  cp  1  x 


+  <o  (o  cos  (a1  +  3*)  —  (to  -h  ijj)^  sin  3  +  (i  +  ^)  cos  3]  (L)-18c) 


T  *  r  [-  sin  o'  (io  +  2  w  t|i  +  to  io  ) 
z  cp  1  '  x  y  r  y  zJ 


+  cos  a'fu  -  2  to  i  ■  «  u  )] 

\  y  X  X  Z  ^  J 


(D-lttd) 


Absolute  Acceleration  of  Pallet  Center  of  Mass  C. 


The  total  acceleration  of  the  pallet  center  of  mass  is  given  by: 


AC  /ground  "  ACp/0p  +  ^Op/ground 


(u-iy) 


Substitution  of  equations  D-12  and  D-L8a  into  the  above  yields  the  following 
expression: 


ACp/ ground 

2  2 

{rcp  [-  u>x  sin  sin  o'  -  toy  cos  83  cos  o'  +  8ln  (<*'  +  B3) 

-  (to  +  cos  B  -  (u  +  <J»)  sin  3]  +  K  } i ' 

Z  Z  X 

2  2 

+  {r  [-  <0  cos  3t  sin  o'  +  to  sin  3-,  cos  o'  -  10  to  cos  (o'  +  3^) 

cp  x  j  y  j  x  y  j 

-  (to  +  t|t)^  sin  3  +  (to  +  \|»)  cos  e]  +  kv}j’ 
z  z  y 

+  {r  f“  (to  +  to  to  )  sin  o'  +  (to  -  to  to  )  cos  o' 

lcpl'-x  y  z '  '■y  x  z ' 


-  2  t|>  («ox  cos  o’  +  u  sin  o’  )]  +  Kz  }R’ 


(U-20) 


Signum  functions 


Before  deriving  the  equations  of  motion  of  the  pallet  and  the  escape 
wheel,  it  is  necessary  to  introduce  two  signum  functions  which  will  be  used  to 
determine  the  directions  of  the  friction  forces  at  the  pallet-escape  wheel  inter¬ 
face  and  at  the  pallet  pivots,  respectively  (ref  l). 
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The  relationship  between  the  contact  point  S  on  the  escape  wheel  and  the 
coincident  point  T  on  the  pallet  is  shown  in  figure  D-3a.  The  signum  function  84 
makes  use  of  the  relative  velocity  Vg/.^;  i.e. 


VS/T 


(D-21) 


The  signum  function  s^,  which  is  associated  with  pallet  rotation,  is 
defined  by: 


8 


5 


(0-22) 


Pallet  and  Escape  Wheel  in  Entrance  Coupled  Motion 


A  free  body  diagram  of  the  pallet  with  the  normal  force  Pn  and  the  fric¬ 
tion  force  uP  acting  on  its  entrance  working  surface  is  shown  in  figure  l)-3a. 
The  normal  an<f  friction  forces  acting  on  the  upper  and  lower  pivots  of  the  pallet 
are  shown  in  figure  D-3b.  Por  verge  nomenclature  see  reference  1. 

Force  Equations  for  Pallet.  The  force  equations  for  the  pallet  in 
entrance  coupled  motion  are  obtained  from  Newton's  law  according  to 


E  1  ■  mp  Aq  /groun(i 
P 


(D-23) 


where  the  acceleration  A,,  ,  .  of  the  pallet  center  of  mass  is  given  by  equa- 

Cp/ground 

tion  i)-2t).  The  sum  of  the  forces  is  obtained  with  the  help  of  the  figures  men¬ 
tioned.  (For  escapement  forces,  see  equation  E-43  of  reference  1.)  Equation  D- 
23  becomes : 

en  “n  +  ul*4Pn"t  +  ‘V  E'  "  ‘'xu'1' 


'  ly'u  V  *  3'  +  VlVu  r 


+  Fx,l  i'  +  Fy,L  j'  +  U^F^  J’  -  PiVy'L  11 


r  2  2 

■  mp  Mr  (-  u>x  sin  sin  a'  -  uiy  cos  B3  cos  a’ 


+  to  u>  sin  (a*  +B,)  -  (w  +  t|>)  cos  B  -  (u>  +  <|>)  sin  B]  +  K  }i' 
x  y  j  z  z  x 
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•oM...^i  Vf  r^.iNArdt  i! 


2  2 

+  {r  [-  u>  cos  0„  sin  a*  +  w  sin  0,  cos  a'  -  u>  w 
cp  x  <j  y  j  x  y 

-  (U)  +  $)2  sin  0  +  (w  +  ij/)  cos  g]  +  K  }j' 

£  2  y 

+  {r  f-fu  +  w  u  )  sin  a'  +  (u»  -  w  «  )  cos  o' 

1  cp1  v  x  y  z  y  x  * 

-  2  $  (u  cos  a*  +  ui  sin  a' )  ]  +  K„}k1 

p  x  y  j 


cos  (a*  +  p^) 


(D-24) 


where 


F’xiu  and  Fyiu  are  the  normal  force  components  acting  on  the  upper  pivot 
*VL  and  tfy'h  act  on  the  lower  Pivot 

f'zi  represents  a  thrustforce  exerted  on  the  pivot  shaft 


Note  that,  as  In  reference  1,  the  force  and  moment  equations  of  the  pallet  are 
given  in  the  x'-y'-z'  system  for  computational  convenience. 

The  unit  vectors  (L  and  nn  are  expressed  according  to  equations  B-5 
and  B-6  of  reference  1  in  the  primed  system  as  follows: 


nt  -  cos  (ij>  +  <x)i'  +  sin  (i|»  +  o)j'  (D-25) 

n  ■  -  3in  (i|»  +  a)i'  +  cos  (i|>  +  «)J'  (U-2b) 


The  angle  a  is  associated  with  the  pallet  and  is  different  for  entrance  and  exit 
contact . 


Substitution  of  equations  D-25  and  D-26  into  equation  D-24  fur¬ 
nishes  the  following  component  expressions: 

x1  -  component  of  force  equation 


-  i*  sin  (vp  +  o)  +  !is,  F  cos  (tji  +  a)  -  F  +  p.s,  F  , 
n  1  4  n  x’u  1  5  y'u 


+  F  -  p,se  F 
x'L  1  5  y'L 


m 


r  r  2 

fr  -  a)  sin  g..  sin  a' 
1  cpL  x  j 


P  1  cp 


2  .  2 
-  to  cos  0,  cos  a'  +  ui  a)  sin  (o'  +  $_)  -  (w  +  <f>)  cos  0 
y  3  x  y  3  z 


-  (iz  4-i|i)  sin  g]  +  KjJ 


(D-27 ) 
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My  ■  sum  of  external  moments  about  point  Up.  It  Is  assumed 

p  that  Up  lies  In  the  plane  of  the  center  of  mass  of  the 

verge  luormal  to  the  verge  pivot  axis).  It  Is  also 

assumed  that  the  forces  f»  and  p.s.P  lie  In  this  plane. 

n  K1  4  n 

A0  /ground  "  absolute  acceleration  of  point  Op  according  to  equation  D- 

p  8  12 

mm 

U  ,  .  ,  the  rate  of  change  of  the  angular  momentum  of  the  verge 
0  X 1  V 1  z  * 

p  J  with  respect  to  point  Op.  This  expression  is  obtained  by 

adapting  equation  b-4  to  the  parameters  of  the  pallet  and 
transforming  the  result  into  the  x'y'z'  system. 


Determination  of  M 


The^jaoments  due  to  the  verge  contact  force  ?n  and  the  associated  fric¬ 
tion  force  UjS^i^  are  taken  from  equations  E-48  of  reference  1.  The  moments  due 
to  the  pivot1  forces,  both  normal  and  frictional,  are  obtained  with  the  help  of 


r‘  **  component  of  force  equation 


Vn  coe  <♦  +  «)  +  wle4  tin  <*  +  a)  -  *y,u  -  *x,u 

+  Fy'L  +  M18S  Fx'L  “  “x2  C08  P3  *in  0,1  +  wy2  #in  c0# 

-  m  w  coe  (o'  +  B-)  -  (w  +  i  )  sin  3  +  (w  +  ip)  coe  $] 

X  jr  •  •  P  *  \ 


+  K  } 

yl 


s'  -  component  of  force  equation 


(0-28) 


“  “p  ^cpt*’  <«x  +  wy“*^  8in  +  “  uxV  c08  °' 

2  (u  cos  a'  +  u  sin  a*)]  +  K  } 
x  y  z 


(d-29) 


Moment  Equation  for  the  Pallet.  The  moment  equation  for  the  pallet  must 


be  written  with  respect  to  the  accelerated  pivot  point  Op: 


where 


h.  -  -  A, 


Op/ ground  p  cp 


x  m  r  (cos  ffl*  +  sin  8J')  +  A.  ,  ,  . 

n  /in  J  M  v  •  ti  •  ' 


0  x'y’z' 
P 


(D-3U) 
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the  figure  D-3b,  The  symbols  p  end  p.  stand  for  the  pellet  pivot  redius  end  the 
pellet  thrust  friction  redius,  respectively."  Therefore, 


\  “  Di  Pn  k'  “  wlBACi  *n  *'  “  Wf  V  k’ 
P 


+  (Lu  *  +  pp  V)  X  (-  Fy,uj*  +  Wls5  Fy,u  i') 
+  Ch,  fe'  +  pp  i')  x  (-  fx,u  i-  -  Wls5  *x,u  I*) 

+  (-4.E'  -Opi'JMVL  i‘  -  Vs  Vi.1'') 

+  (-  &'  «*  p  I')  X  (fx.l  I'  +  v1b5  fx,l  j') 


(l>-31) 


The  above  becomes: 


*\>  "  lLu  V»  +  Lu  “1*5  Vu  +  Lt  'y'L  +  LL  Vs 

P 

+  K  Vi  Ky'“  -  Lu  Vu  +  \  “1*5  fy'L  *  LL  V  JP 
+  K  (Di  -  V4C1>  -  U1  V 5V  -  Vl*5  Vu 
-  °pVs  Vu  -  Vi’s  Vl  ‘  Vl*5  VjE' 


(i>-32) 


Determination  of  -  A  .  .  x  m  r  (cos  0  i*  +  sin  0  J') 

Op/ground  p  cp  _  J 


With  the  help  of  equation  D-12,  for  the  above  cross-product  the  following  is 
obtained: 


-  (kx  i'  +  Kv  j'  +  K„  k’)  x  mprcp  (cos  0  i'  +  sin  0  j') 


-  m  r  K  sin  0  i'  -  m  r  K  cos  0  j’  -  m  r  (K  sin  0  -  K  cos  0)k' 
p  cp  z  p  cp  z  p  cp  x  y 

(D-33) 


M  keference  8  for  determination  of  thrust  friction  radius,  p  268. 
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Determination  of  lL 

**p*  y  * 

As  Meted  earlier,  equation  b-4  must  first  be  adapted  to  the  pallet  fixed 
coordinate  system  with  pallet  related  nomenclature.  This  leads  tot 

Ml 

\  •  i1^  \*  vc  (I«s "  H1  +  1{rt  4") 

+  i5 ^  («e+  «*,«„)  +  iw  (%  -  »t  )F^ 

+  [ln^  “i  +  Vc  (IC^  '  *5^, )+  InS  (*|V  “cJ 
-  ^  (“c+  Vt)-1«  K2  ■  “C2)P^, 


+  [I 


*s 


V 


Vi  (tn^,  '  ^  (\V 

Vc1'  (“c2  •  “i  ^Cp 


(D-34) 


The  angular  velocities  and  accelerations  of  the  pallet  are  now 
expressed  according  to  equations  A-21  to  j4-23_and  equations  A-25  to  A-27 ,  respec¬ 
tively.  Subsequently,  the  unit  vectors  n  ,  n  ,  and  n  are  substituted  according 
to  equations  A-8,  A-9,  and  A- 14. 


V 


These  operations  result  In  the  following  component  expressions  for 


x'y' z*  * 


where 


H0  x’  *  A1  +  A2  *  +  A3  *  +  \  m 


”0py*  "  A5  +  A6  *  +  A7  $  +  A8  ♦ 


"  A0  +  A1q  ^ 


*0  z'  9 
P 


(D-35) 
(D-36) 
(D-37 ) 


A^  -  cos  B  {-  (<» ^  cos  o'  +  uy  sin  a'  ) 

+  fl  -  I  1  a)  f  to  sin  a'  -  tu  cos  a'  ) 
nT\»  z  x  y 

-  I,  [to  fu  cos  o'  +  u»  sin  o'  1  +  (u  sin  o'  -  u>.  cos  a'] 
*•  z  '  x  y  x  y 
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+  I  .  [(i^  eo.  o'  +  ay  »tn  a'  )(uy  sin  o'  -  ay  co.  a*  )  -  ay  ] 
p 


-  In^  [(‘^  «in  tf  -  ay  co.  o'  f  -  ay2  ]} 


-  tin  B  {I  [uy  «i«  ^  608  «’  ] 

p 


-  (l„  -  I„)  «y  co.  a’  +  uy  .in  a*  ) 

p  v 

-  I  [(ay  cot  a’  +  ay  Bin  a’  )(ay  sin  tf  -  ay  co.  o'  )  +  ay  ] 

nS> 

+  I  [(ay  cos  a'  +  ay  .in  a’  )  “  ay  («y  .In  o'  -  ay  co.  a'  )] 
p 

(D-38) 


-  1^  K2  “  (“x  C08  +  <*y  8in  «'  f  ]} 

A2  “  K  8in  “  “y  C08  a'  +  *«p  “  608  0+  2  x5iy  sin  ^ 


-  (ay  cos  o'  +  ay  .in  a’  )  [2  I  ^co.  $4-  (i  ^  -  I  ^  +  1  ^ )  8ln  •] 


+  2  oj 

(i  co.  8  + 

IrF 

2 

ns> 

?5P 

A.  "  I 

cos  8  +  I 

sin  8 

3  nCp 

55» 

A  -  I 

.in  8  “  I  -c 

cos  8 

4  n?p 

!5p 

(D-39) 


(D-40) 

(D-41) 


sin  B  {  -  I ^  (ay  cob  a'  +  ay  sin  o'  ) 


+  fl  -  I  ]  a)  (a>  .in  o'  -  as 
1  nnJ  z  y 


cos 


o'  ) 


+  X  [—  a>z  ( cos  o'  +  ay  tin  at  )  -  (ay  sin  o'  -  ay  cos  a  )] 

P 


-I  ["  (ay  CO.  o'  +  ay  .in  o’  )(ay  .in  at  -  ay  co.  a’  )  +  uz  ] 
C> 
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r! V,V,V*V.  <  •%'sV.rV*%  ’■***•" 


-  In^t(\  tin  rf  -  ^  cot  o'  )2  -  wR2  ]} 


+  cot  8  {  I  tin  o'  -  cot  o'  ) 


-  -  I5^]  ^  («  cot  o’  +  sin  o’  ) 


+  I  C0B  0,1  +  “V  8in  0,1  )(“ic  ain  “  °y  008  )”  ] 

p  y 


-  I_  f-(i  cos  o'  +  u  sin  o'  )  +  tu  (w  sin  o'  -  u  cos  o'  )1 
5n  *■  v  x  y  ^  ss  v  x  y 


”  Ic5p  K2  “  K  C0B  *'  +  "V  8in  °'  f  H 


(D-42) 


A,  *  fw  sin  o'  -  a)  cos  o'  )[(l,.,  +  I  -  I  )  sin  0  -  2  I_  cos  0] 

6  *  y  1  «p  %  J 

+  fw  cos  o'  +  0)  sin  o'  Iffl  -  I...  +  I  „  )  cos  0>  2  I.  sin  0] 

*  y  %  ^  V  ^  1 


♦2*V  V 

A?  -  I^„  sin  0  -  I  rr  cos  0 


‘S  "  'S 


*J  ‘  -  &;£  8l"  6+  ^  81)9  ») 


CD-43) 

(D-44) 

(D-45) 


Ag  -  I?c  wz  -  [I  -  1^  ]((Vos  o’  +  u  sin  o'  )(uxsin  o'  -  iuycos  o'  )] 


+  I  -  [m  (id  sin  o'  -  u>  sin  o'  )  +  (u  cos  o'  +  o>  sin  o'  )] 
s  s»_  z  x  y  x  y 


-  I  _  [(to .  sin  o'  -  w  cos  o'  )  -  u>  ( to  cos  o'  +  w  sin  a'  1] 
nCp  1  v  x  y  J  z  K  x  y  J  J 


-  I_  f  f to  cos  o'  +  to  sin  o'  l2  -  (o»  sin  o'  -  oj  cos  o’  )2  1  (D-46) 

Cn  lv  x  y  '  '  x  y  7  J 


A10  “  XCC 


(D-47 ) 
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simplification  of  Force  and  Moment  Equations.  In  order  to  be  able  to 
solve  for' the  upper  end  lower  pivot  forces,  both  the  force  and  moment  component 
equations  are  now  rewritten  in  an  appropriate  simplified  form. 

x* -Component  of  the  Force  Aquation 
liquation  D-27  becomes! 


-Px.u  +  Au  Py«u  +  ?xtL  “  AH  Fytt 

■  *u  +  *13  i+  v  *“  +  *15  *  +  Pn  *16 


A12  -  mp  rcp  [-  u>x2  sin  Bj  sin  tf  -  uy  cos  ftj  cos  o’ 


+  u^ui  sin  (o'  +  ftj  )  -  Uj2  cos  B  ■“  w,,  ®ln  B]  +  ®p  Kx 


A13  *  “  2  “*  “p  rcp  coa  6 


A,,  ■  -  m  r  cos  B 
14  p  cp 


Aie  *  -  m  r  sin  B 
1“  o  cp 

A16  *  "  [yis4  008  (’l,+  °0  “  8in  ^  + 

y' -Component  of  the  Force  Equation 
Equation  IV- 28  becomes: 

-Au  Fx,u  -  Fy,u  +  An  Fx,l  +  Fy.L 

-  a17  +  a18  4>+  a19  t2  +  a20  i+  a21  Pn 
2  2 

A17  *  mp  rcp  1“  “x  C08  h  8in  *  "V  8in  H 

”  °k(JV  008  (“'  +  %  )  ”  sin  B  +  ujj  cos  p]  +  mp  Ky 

Air  -  -  2  mp  rcp  u>8  sin  B 


cos  a' 


(D-48) 


(D-49) 


(D-50) 

(0-51) 

(D-52) 

(D-53) 

(D-54) 


(D-55) 


(D-56) 

(0-57) 
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Aio  "  "  b«  r„«.  *in  8  CD-58 > 

**  p  cp 

•  %  rcp  coa  »  (D-59) 

A21  *  “  [co*  (♦+  °)  +  sin  (♦+  ««)]  (D-60) 

^-Component  of  the  Force  Equation 

Equation  T>— 29  is  rewritten  to  read: 

f«>  ‘  *22  +  *23  *  (D-61) 


The  tilde  is  now  used  to  indicate  the  conservative  nature  of  this  force,  when  the 
terms  A22  and  A23  are  tade  absolute. 

Thus 

A22  *  l“p  rCp["  (“x  +  uy“  )  8ln  «'  +  (“  -  )  cos  1  +  %  KJ  <D~62> 

and 

A23  -  I-  2  mp  rcp  (<^  cos  o’  +  sin  o'  )|  (D-63) 

The  absolute  values  in  the  above  expressions  will  be  useful  later  (eq  D-123). 
x'  -  Component  of  Moment  Equation 

The  x'-component  of  equation  D-30  is  obtained  with  the  help  of  the  x'- 
components  of  equations  D-32  and  D-33,  as  well  as  equation  D-35.  Therefore, 

T'u  All  Fx'u  +  Fy'u  +  k.  Ali  fx'L  +  Lh  Fy'L 

-  m  rcp  Kz  sin  $  +  Aj  +  }  A2  +  ?  ^  +  J  A4  (D-64) 

y' -Component  of  the  Moment  Kquation 

The  y’-component  of  equation  D-30  is  obtained  with  the  help  of  the  y'- 
components  of  equations  D-32  and  D-33,  as  well  as  equation  D-36: 

Ki  Fx'u  +  Lu  A11  Fy'u  ”  FxfL  +  W.  A11  Fy'L 

"  "  mp  rcp  Kz  C08  e+  A5  +  A6  $  +  A7  ?  +  Ag  ip  (D-65) 
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s' -Components  of  the  tenant  Equation 

The  i' -component  of  equation  D-30  is  composed  of  the  z' -components  of  equa¬ 
tions  D-32  and  D-33,  as  wall  as  equation  l>-37 : 


P  (uj  -  p  a.  C' )  -  p  A  P  ,  -  p  A, ,  f  .  -  p  A. .  P  . 

n  v  1  14  1'  f  11  a'  p  11  y'u  p  11  x'u 


-  p  A. ,  V  , ,  -  p  A, ,  P  , , 
p  11  y'L  p  11  x'L 


-  «p  rcp  (Kx  «i«  ®  *  *y  coa 


+  Ay  +  A10  4» 


(D-6b) 


Solution  for  Pallet  Pivot  if orcea.  The  force  Px'u,  f’y'u*  f’x'L*  and  *'y'L 
are  obtained  from  tne  simultaneous  solution  of  equations  D-4tt,  D-55,  l)-64,  and V 

65.  The  force  PE,  given  by  aquation  D-61.  These  five  forces  are  eventually 

substituted  into  equation  D-66,  and  the  resulting  expression  is  solved  for  the 

contact  force  P„. 

n* 

The  simultaneous  set  of  equations  becomes 


where 
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l 

~A11 

-1 
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LuAll 
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U.Ail 

~Lu 

WAn 
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f'x'u 

“pi 

1 

Py'u 

m 

tfp2 

ll 

f’x'L 

UP3 

LLAu 
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Py'L 
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Bp4 
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1P1 


-  ♦  +  Ai/.  +  A,t  ♦ 


12 


13 


14 


15 


P  A, 
n  16 


(l)-b/ ) 


(D-bd) 


8p2  “  a17  +  Aiu  ^  +  Ai9  ^  +  A20  *  +  Fn  A21 


(D-6y) 


Bp3  ■  mp  rcp  \  8ln  8  +  A1  +  A2  *  +  A3  ^  +  A4  ♦  (D-70) 

*p4  “  “  mp  rcp  Kz  cos  B  +  A5  +  A6  *  +  A?  $2  +  Ay  *  CD-71) 


Cramer's  rule  will  now  be  used  to  determine  the  four  pivot  forces 

*x'u>  f’y'u.  f'x'Lt  and  Py'L*  To  this  end,  the  coefficient  determinant  D  must  be 
found  first. 
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Evaluation  of  the  Coefficient  Determinant  D 


The  coefficient  determinant  of  equation  l>-67  la  given  by: 


-l 

All 

1 

-An 

“All 

-1 

All 

1 

LUA11 

Lu 

lLa11 

Ll 

-Lu 

LuAll 

“lL 

llaii 

(D-72) 


Evaluation  of  the  above  furnishes: 

2u2 


»-[(>•„  +  \)  U  +  \S)Y 

Since,  according  to  equation  D-49 


(D-73) 


A  m  u  s 
11  yl  5 


(D-74) 


and  is  always  equal  to  unity  (eq  0-22),  the  coefficient  determinant  becomes 


2  1 12 


(D-75) 


Evaluation  of  Pivot  Force  F  . 

xu 


The  pivot  force  Fx<u  i8  obtained  from 
D, 


*  ...  “ 


Vu 


x‘  u  D 


(D~76) 


where 


px'u 


“pi 

A11 

1 

-An 

“p2 

-1 

Ali 

1 

»p3 

Lu 

LlAh 

Li 

v 

LuAii 

-ll 

llah 

(U-77) 
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Evaluation  of  0V  furnishes: 

*x'u 

°F ,  -  (Lu +  h.)  (1  +  An2)  t-  \  v  -  *11  H  v  +  *11  flp3  -  VJ  (|W8) 

X  u 

After  substitution  of 

Au2  -  u2  (D-7!/) 

the  following  is  obtained 


B,  ,  ■  ('•u  +  h.)!1  +  \  »pl  -  Au  \  Bp2  ♦  Au  »p3  -  bpJ  OMtW 


x'u 


Subsequently,  equations  D-4y  and  D-68  to  D-71  are  substituted  into  the  above  and 
the  coefficients  of  similar  terms  are  collected  and  made  absolute.  The  latter  is 
done  to  get  conservative  pivot  and  friction  forces.  This  leads  to: 


where 


x'u 


p  cp  z  '  1  5 
C2  I,ll  A13  +  Vs  K  ‘ 

^3  “  I"  ^14  +  ^1*5  (A3  ”  Aiy  )  “  A7 

C4  ■  I"  h.  A15  +  V5  'A4  H,  A20^  *  \ 

C5  “  I*  LL  Alb  ”  V5  LL  A2ll 


+  C2 

|  +  c3  *2  +  c4  J  +  c5  fn] 

(D-81) 

,*17) 

-  *3 

+  cos 

e)l 

(D-82) 

A12 

(D-83) 

1  Aiy) 

-  A7l 

(D-S4) 

A20  ^ 

-  Ay| 

(D-85) 

(D-86) 

finally,  substitution  of  equations  1W5  and  D-81  into  equation  D-76  gives  the  now 
tllded  pivot  force 


x'u 


(Lu  +  I^JCl  +  V\) 


(Cl  +  C2  *  +  C3  f-  +  C4  ’i  +  C5  pj  CD-8/) 
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Evaluation  of  Pivot  force  7  , 

y'u 

The  pivot  force  Pyiu  is  again  obtained  with  Cramer's  rule,  i.e., 
DP 

Xjl 


'y'u  D 


where 


(D-88a) 


Dk 
“  y'  u 


-1 

BP1 

1 

“An 

~Au 

Bp2 

An 

l 

LuAl  1 

bp3 

LLa11 

ll 

_Lu 

Bp4 

-!*L 

lLa11 

(D-88b) 
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Evaluation  of  Dw  furnishes: 
'y'u 


UE  ,  -  (Lu +  +  An  )  Ki  *1  v  -  h.  V  +  *»P3  +  AU  v] 

y  u 


(D-89) 


and  again,  with  A]^  • 


3 


UVu  ■  (Lu  +  \.>(l  +  "l)  (“l*5  \  “pi  -  h.  V  +  “p3  +  "l*S  V> 


(D-yu) 


M 

H 


Appropriate  substitution  into  equation  D-88a  and  proceeding  in  a  manner  parallel 

to  that  followed  in  the  determination  of  7  ,  ,  the  following  is  obtained  for  7  . 

x'u  y'u 


where 


<v  J 

,,u"Kr+-hJ(r+ut) 


IT  (C6  +  S  {  +  %  *2  +  S  *  +  U10  '.I  (^,U 


C6  ■  lAi  "  h,  A17  +  w185  (lL  a12  +  *5)  +  mp  rcp  Kz  (sin  0  "  ^1*5  008  p)l 

(U-V2) 

c7  "  Ia2  “  lL  A18  +  ul95  +  LL  a1  3 )  I  CU-93) 

C8  "  lA3  ”  A19  +  ul®5  (A7  +  LL  A14 )  i  (U-94) 
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“  1^4  "  A20  +  wl®5  A15  A8^ 

c10  "  l^is5  LL  A16  “  \  A21 I 


(D-95) 

(D-96) 


Evaluation  of  Pivot  Force  Fxil 


’•j 


The  pivot  force  Fx»l  is  obtained  from 
D. 


Fx'L  - 


Fx'L 


where 


<Vl 


-1 

A11 

Bpl 

"Ali 

“Ml 

-1 

Bp2 

1 

luMi 

Lu 

Mp3 

ll 

~Lu 

V*ll 

Bp4 

llMi 

(D-97) 


(D-98) 


Evaluation  of  D_  furnishes: 
fx'L 

»F  ,  ■  (L„  +  't-Hl  +  An2)  (>■„  »pl  +  \  A1I  V  +  A11  “p3  -  V> 

x’L 


(D-99) 


and  again,  with  *  s^u^ 


(Lu  +  Lj(l  +  uj)  (Lu  B  +  V5  Lu  Bp2  +  V,  BpJ  -  Bp4) 


Fx'L 


(D-100) 


Proceeding  as  before  to  obtain  ?  ,  ,  the  following  is  found: 

X  Li 


3.1 

ji 


Fx'L 


(Lu+  h.)(l  +  »l) 


[Ca  +  C12  *  +  C13  J,2  +  C14  i  +  C,3  F„1  <»-101) 


where 


C11  *  K  A12  ~  ^  +  A17  +  Al)  +  mp  rcp  KzK85  8ln  B  + 


c12  ‘  !L0  a13  ~  A6  +  Vs  ^Lu  A18  +  A2^ 


(D-103) 
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Evaluation  of  Pivot  Force  F 


y'L 


The  pivot  force  FytL  is 

obtained 

from: 

df 

F  ,  - 

ry'L  D 

where 

-1 

A11 

1 

BP1 

D-  - 

-AU 

-1 

A11 

Bp2 

Fy'L 

LUA11 

Lu 

LLA11 

®p3 

“Lu 

LuAn 

“lL 

Bp4 

Evaluation  of  D„ 

*  y  L 


furnishes : 


(D-104) 

(D-105) 

(D-106) 


(D-107) 


(D-108) 


df  ,  ‘  (H,  +  H.)t>  *  *n2)[-  K,  *11  v  +  S,  V  +  V  +  *11  BP?1 

y'L  tu-iu*#; 


and  again,  with  Ajj  “  S5,J1 


°F  "  <L„  +  h.K1  +  »!)  I"  “1*5  S.  "pi  +  Lu  Bp2  +  V  +  V5  “pJ 

y'L  QD-llu; 


F  Is  found  from  equation  D~107  in  a  manner  shown  earlier: 

y'L 

•2 


p  .  — - - - - [C.,  +  Cl7  ♦  +  C. jj  r  +  C  g  t  +  C  P  ] 

y  (L„+  Ll)(1  +  »{)  16  U  ‘8  19 


(D-lll) 


°16  -  !>•„  *17  +  *1  +  V5(*5  -  \  *12>  +  “p  rcp  Kplsin  6  -  “195 


where 


(D-113) 


C17  "  \\  Aifi  +  A2  +  gl85  ^A6  ”  Lu  A13  )l 

C18  *  K  Ai9  +  S  +  Vs  -  Lu  Au)|  (D-114) 

C19  "  K  a20  +  A4  +  wl®5  (A8  “  Lu  Al5)l 

C20  "  lLu  A21  ”  W1S5  Lu  A16 I  (D-116) 


Substitution  of  Conservative  (Tilded)  Pivot  Forces  Tnto  the  z1 -Moment 

Equation.  Rather  than  substitute  the  forces  Fxiu,  Fy»u,  FX»L,  FytT,  and  Fzi 

Into  the  z' -moment  equation  D-66,  the  associated  t lined,  conservative  expres 

slons,  as  given  by  equations  D-6 1 ,  0-87,  D-^91,  D-101^  and  0-11 1  are  used.  To 

simplify  matters,  let  the  sum  of  F  ,  ,  F  ,  ,  F  ,  and  F  ,,  be  first  determined: 

r  x’u’  y'u*  x'L*  y'L 


F  ,  +  F  ,  +F,T+FfT 

x'u  y'u  x'L  y’L 


where 


•2 


A24  +  A25  ♦  +  A26  *  +  A27  ♦  +  ^8  Pn 


*24 


VC6+  C11  +  C16 

^  (1  +  »?r 

c2  +  c7  +  C12  +  C1? 

'25  “  ;  t, 

4  (1  +  wj 


A26 


C3  +  C8  +  Ci3  +  C18 

\  (l  +  pj) 

C4  +  C9  +  Cl4  +  c19 


27  lt  (1  +  u2) 


A28  “ 


C5  *  C10  *  C15  +  C20 

4  (1  +  U2) 


and 


(0-117) 

(D-118) 

(D-119) 

(D-120) 

(D-121) 

(D-122a) 


Lt  -  Ou  +  Ll  (0-122b) 

Substitution  of  the  above,  as  well  as  equation  D-61  into  equation  D-66,  and  let¬ 
ting  Ajj  ■  u^s-  according  to  equation  D-49  leads  to  the  following  provisional  z'- 
moment  expression: 
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Pn  (°1  "  W1#4C1^  “  pf  V5  (A22  *  A23  ^ 

“  pp  wlg5  1^4  *  ^  ^5  *  ^  A26  *  ^  ^7  +  Pn  A28^ 

•t 

“  ~  «p  tcp  (**  »ln  8  -  Ky  co»  8)  +  Ag  +  Aiq  <|> 


CD— 1 23) 


To  make  sure  that  all  friction  momenta  act  In  a  direction  opposite  to  the  instan¬ 
taneous  rotation  of  the  pallet,  the  eigne  of  thoee  friction  terma  which  depend  on 
i|i,  |i,  or  i|>  have  been  left  undecided  for  the  moment.  They  will  be  reaolved  below. 

Before  these  decisions  are  made,  let  equation  0-123  be  rewritten  as 

follows: 


Pn  “  C[  V4  “  PP  Ml*5  A28^  "  Ml85  ^Pf  A22  +  Pp  A24^ 

*  yis5  [Pf  A23  >  pp  ^5^  *  *  pp  wl®5  A26  ^  *  pp  yl®5  A27  ^ 
•• 

"  Al0  ♦  +  \  ‘  “p  rcp  8in  $  “  Ky  cos  ^ 


(D-124) 


With  83  positive  for  positive  rotation  of  the  verge  and  vice  versa 
and  with  all  other  parameters  positive  at  all  tines,  the  following  moment  compo¬ 
nents  of  equation  D-124  must  have  negative  signs  during  positive  rotation: 


'  Pn  %  V5  *28  (D-‘25) 

-  V,  *22  +  %  *24]  CD-‘26) 

"  pp  V5  *26  *2  (D-127) 

The  sign  of  the  terra  containing  $  must  be  negative  for  a  positive  $  and  vice 
versa.  Therefore,  the  sign  of  $  can  be  used  to  control  the  sign  of  this  term, 
and  the  signum  op«  *ator  S5  may  be  omitted.  This  term  becomes: 

-  M1  [P£  a23  +  Pp  ^5]  *  (D-i2K) 

The  choice  of  sign  for  the  term  containing  the  pallet  angular  acceleration  is 
discussed  in  detail  in  appendix  F  of  reference  2.  This  work  leads  to  the  compu¬ 
tational  rules  of  equations  D-134  and  D-135  below.  These  rules  deal  with  the 
sign  in  the  effective  moment  of  inertia  IpR.  (Note  that  the  signum  function 
S3  has  been  omitted  in  these  expressions.) 
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With  the  above  considerations,  equation  D-124  becomes: 


J 


J 

i 


where 


A 


Pn  A29  "  A30  “  Si  *  “  A32  ♦ 

"  IPR  ^  +  N  ”  mp  rcp  ^Kx  ,ln  0  ~  Ky  co#  0) 

(D— 129) 

a29  *  D|  "  ci  hjS^  -  pp  PjSj  A28 

(D-130) 

A30  "  Vs  (Pf  Aj2  +  pp  A24) 

(D-131) 

A31  *  »*!  (P*  A23  +  Pp  ^5^ 

(D-132) 

A32  "  V5  pp  A26 

(D-133) 

•  •• 

IpR  -  1^  +  A^jj ,  when  41  and  ij>  have  identical  signs 

(D-134) 

.  •• 

IpR  "  *££  "  ^33*  when  ♦  and  ♦  have  opposite  signs* 

(D-135) 

333  "  Pp  A27 

(D-136) 

Equation  D-129  is  now  rewritten  to  find  an  expression  for  the  contact  force  Pn: 


P 

n 


IPR  ♦  +  *>  +  *30  +  *31  ♦  *  *32  ~  “p  rcp  (K,  «  -  Ky  °°»  ») 

(D-137) 


'29 


The  above  expression  is  now  changed  to  reflect  the  escape  wheel  angular  velocity 

and  angular  acceleration  $  and  respectively,  so  that  it  may  later  be  equated 
to  ah  expression  for  the  escape  wheel.  Equations  B-19  and  B-26  of  appendix  B  of 
reference  1  show  the  following  relationships: 

ij/  »  q>  U  (D— 138) 


*  Ip  -  A333  oust  not  become  negative, 
zero. 


If  this  occurs  Ipg  must  be  set  equal  to 
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V  •  • 

A.  JO. 


'i^^r*.*  ■  *.  r  i%jn.  *  .  -  *y  m  r.;.  v  %  i  #7 ,'  ’rv%"  *.’*■„ ^tt  i  v  •W7"£W*"^V5*;  rcr.’rv*;  T^r^..  iwvr.'**~v  .t,; 


*•  .  #j 

♦  -  u  4  ♦  v  r 

U  and  V  may  ba  obtained  from  rafaranca  1.  Thla  leads  tot 
'n  *  iq;  t1™  »  ♦  +  (*32  »2  +  *FR  V>  *’ 


(D-139) 


+  *31  u  ♦  +  A,  +  -  mp  rcp  (K*  win  8  -  Ky  cos  $)] 


(0-140) 


Fores  Equations  for  Escape  Wheal  and  Pinion  No.  3.  How  the  contact 
forces  P_  end  P53,  together  with  their  associated  friction  forces,  act  on  the 
escape  wneel  ana  pinion  no.  3  combination  is  shown  in  figure  D-4a.  A  separate 
free  body  diagram  of  the  pivot  shaft  of  the  escape  wheel  is  shown  in  figure  D- 
4b.  All  these  forces  are  now  defined  in  the  projectile  (fuze)  fixed  X-Y-Z  sys¬ 
tem.  This  makes  it  necessary  to  transform  the  unit  vectors  nt  and  nn  from  the 
x'-y*  to  the  X-Y  system  (eqs  D~25  and  D-26  as  well  as  eqs  B-79  to  B-82  of  refer¬ 
ence  1). 

Since 


and 


i'  -  -  cos  0-  i  -  sin  8«»  J 


j '  -  sin  $3  I  -  cos  03  J 


(D-141) 


(D-142) 


the  previous  unit  vectors  become: 


n  ■  -  cos  (i|i  +  a+03)I  -  sin  (^  +  a  +  03 )  J 


n  ■  sin  (t|>+a+B,  )I-  cos  (\ji  +  a  +  B.>)  j 


(D-143) 

(D-144) 


The  force  equations  for  the  escape  wheel  in  coupled  motion  are 
generally  obtained  from  Newton's  law: 


E  F  -  m3  /ground 


(D-145) 


where 


E  F 


■  sum  of  pivot  forces  as  well  as  contact  forces  Pn  and  F23 
and  their  associated  friction  forces 


m3 


mass  of  escape  wheel  and  pinion  no.  3 


view  free  body  diagram  of  escape  wheel  and  pinion  no.  3  in  entrance  coupled  notion 


Escape  wheel  rotates  ccw 


TE,  *  ^*£3**,  (Thrust  friction 
3  3  torque) 


A  ,  .  ■  acceleration  of  escape  wheel  center  of  mass,  which  lies  on 

g/grou  axis  of  rotation,  with  respect  to  ground.  Therefore 


*0  /ground  ^0  /C  +  ^C/ ground 


C  D— 146 ) 


In  the  above,  \j/gr0und*  t*'°  acceleration  of  the  fuze  geometric 

center  C  with  respect  to  the  ground  is  given  in  terms  of  the  X-Y-Z  system  by 
equation  C-4  of  appendix  C.  The  acceleration  of  the  escape  wheel  center  of  mass 


with  respect  to  the  above  point  C,  i.e.,  A  .  ,  becomes: 

08/c 


°8/C  X 


x  (ux  x  n^ )  +  w  x 


where 


n,  -  cos  y-  i  +  sin  Yo  j 


Substitution  of 


ui  -  ux  i  +  wy  j  +  wz  k 

according  to  equation  A-l  and 

<o*u>  i  +  i  j+w  ic 
x  y  z 

according  to  equation  A-5,  with 
3x  "  “3  “*  T3 


3,  '  fl3  8ln  V3 


results  in: 


Art  J  i  +  J  j+J  k 
0  /C  x  y  J  z 

8 


where 


7  1  , 

J  *  (jj  w  -  fw  +  oj  )  <R~  -  a)  <R- 

v  ^  ^  v  %  '  3x  z  3y 


x  x  y  3y 


7  7  • 

J  ■  (I)  b)  <?„  -  f U)  '  +  U)  )  (R-  +  U) 

y  x  y  3x  x  z  ;  3y  z  3x 
Jz  »  (wx<R3x  +  “y  ^3y )  l0z  4'  “x  fl3y  "  “y  ^3x 


(0-147) 


(D-148) 


(D-149) 


(D-150) 


(0-151 ) 


(0-152) 


(D-153) 


(D-154) 


(D-155) 


(0-156) 
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Equation  I>-1 33  la  now  subatitutad  together  with  equation  03  Into  equation  D-146: 


where 


7L  ,  .  -  N  I  +  N  J  +  N  t 

0  /ground  x  y  J  t 

s 


Ny  ■  Civ  *4“ 


Ny  m  Gy  +  Jy 

N_  -  G,  +  J. 


(D-157) 

(D-158) 

(D-159) 

(D-160) 


The  vectorial  force  equation  Is  now  obtained  with  the  help  of  fig- 
ures  D-4a  and  D-5b,  and  equation  D-145: 


“Pnnn  "  ul®4  Pn  nt  +  F23  n23  +  u  ®2  P23  ^23 
+  Fz3  £  “  Fx3u  1  ’  Fy3u  j~  "  U  Fx3u  +  11  Fy3u  1 
+  Fx3L  *  +  Fy3L  J  +  y  Fx3L  ^  "  u  Fy3L  * 


(Nx  I  +  N  j  +  Nz  k)  »3 


(D-161) 


nt  and  nn  are  now  substituted  according  to  equations  D-143  and  D-144.  The  unit 
vectors  ^3  and  nN23  are  taken  from  reference  1,  l.e., 


‘23 


sin  (b_  +  0- )  1  -  cos  (p,  +  02)  J 


nN23  "  COfi  +  02  )  *  +  ®in  ^2  +  02^  ^ 

This  leads  to  the  following  force  component  equations: 

-P  sin  (if;  +  o  +  8^ )  +  PjS^  cos  (\[i  +  a  +  0^) 

+  F23  sin  (B2  +  02)  +  u  s2  F23  cos  (8,  +  0, )  -  F 


+  g  F  ,  +  F_T“ViF,,  “N  n„ 

y3u  x3L  y3L  x  3 


x3u 


(D-162) 

(D-163) 


(D-164) 
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&| 


11 


Pn  cos  (tj»  +  a  +  03)  +  u^  Pn  «in  (♦  +  a  +  &3) 

-  f23  co.  (b2  +  e2)  +  w  .2  f23  (e2  +  e2)  -  Fy3u 

“  “  Fx3u  +  Fy3L  Fx3L  -  Ny  ”3 
F.3  “  N.  m3 


CD-165) 


(D-166) 


Moment  Equation,  for  the  R.cape  Wheal  and  Pinion  No. _ 3^  Since  the 

escape  wheel  and  pinion  no.  3  represents  a  symmetrical  body,  ita  moment  equation 
may  be  expressed  in  terms  of  the  projectile-fixed  X-Y-Z  system  according  to  equa¬ 
tion  B-13  of  appendix  B.  Adaptation  of  this  expression  to  the  escape  wheel  sys¬ 
tem  furnishes: 


M  •>  f  I  U>  +  I  tii  f  «i)  +  A 1  -  I  u,  U.J  i 
0  L  xs  x  as  y  '  z  ys  y  zJ 

s 

+  iV  "y  +  'x.  “x  “z  -  ‘.i  “x  K  +  ♦>]  5 

*  I,.  +  ♦)  E 


(D-167 ) 


The  moment  sum  Rq  about  point  08  is  now  found  with  the  help  of  the 
free  body  diagrams  of  figures  D-4a  and  l)-4b.  Reference  l  is  also  useful. 


>\>  ■  -PJA1  -  B1  V*)S  +  rb3  F23  5  '  w7.  ( d2 -  *2)  F23  F  * 

s 

+  (Lu  k  +  p3  j)  x  (-Fy3u  J  +  v  Fy3u  i) 

+  ^Lu  ^  +  p3  x  Fx3u  *  ”  v  Fx3u  ^ 

+  (-Ll  k  -  p3  j)  x  (Fy3L  J  “  V  Fy3L  I) 

+Ki-  P3  i)  *  (Fx3L  1  +  M  Fx3L  5) 


M  Pf3  Fz3  k 


(D-168) 


The  term  -u  pf3  Fz3  represents  the  thrust  friction  moment  due  to  force  Fz3  (eq  D- 
166).  The  term  p^3  stands  for  the  thrust  friction  radius  of  the  escape  wheel 
pivot.  Equation  D-168  becomes: 
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-  K  F,3. 

t  +  Lu 

“  f«3«  +  L: 

L  Fy3L 

+  ll 

**  Fx3lJ  1 

+  [Lu  U 

Fy3u“ 

Lu  Px3u  + 

Ll  p  1 

Py3L  " 

LL  Fx3lJ 

J 

+  [-  p 

L  n 

(a;- 

>1  »,%)  + 

rb3  F: 

23  “  V 

1  S2  (d2  - 

*2 

)  F 

J  23 

-  u  Pf3 

?z3  ~ 

p3  “  V  ' 

*  P3  W 

F  , 
x3u 

-P3  PF 

y3L 

p3  11  Fx3l1 

lc 

(D-169) 

Substitution  of 

equation 

D-169 

into 

equation 

D- 

167  leads 

to  the 

following  moment  component  expressions: 


L  iiF,  +  L 
u  x3u  u 

F  ,  + 

y3u 

h.  »■ 

F  + 

x3L 

\  Fy3L 

*  t  Oi  + 

1 

to  (to 

+ ♦) 

-  I 

(0  (0 

XS  X 

zs 

y  v  z 

ys 

y  z 

-L  F  _  + 
u  x3u 

L  | 
u 

"  Fy3u  • 

■  LL  Fx3L  + 

LL  W  Fy3L 

“I  10  + 

I 

(0  (0  “  I 

to  f(0 

+  ♦) 

ys  y 

XS 

X  z 

ZS 

X  v  z 

Pn  (A'l  “  B1  VJ  +  F23  ^rb3  V  ®2  ^d2  “  ®2^ 


-UP 


-  P 


f  3  z3  H  3  1  x3u  y3u  x3L  y3L 


I  (to  +  a) 
Z8  K  Z  J 


(D-170) 


(D-171 ) 


(D-172) 


Simplification  of  Force  and  Moment  Equations.  To  solve  for  the  pivot 
forces  Fx3u,  Fy3u,  Fx3l,  and  Fy3L»  the  X  and  Y  components  of  the  force  and  moment 
equations  must  be  rewritten  in  an  appropriate  form. 


X-Component  of  Force  Equation 


Equation  D-164  becomes 


where 


Fx3u  '  “  Fy3u  "  Fx3L  +  v  Fy3L 
"  p„  *33  +  F23  *34  +  *35 


A33  -  cos  (i|»  +  a  +  33)  -  sin  (<|»  +  a  +  83) 


(D-173) 


(D-174) 


A34  ■  sin  (32  +  63)  +  w  s2  cos  (02  +  02^ 


(  D—  i  7  5 ) 


V35  “  ~  \  m3 


(D-176) 


Y-Component  of  Force  Equation 
Equation  D-165  becomes: 


where 


U  Fx3u  +  Fy3u  "  P  Fx3L  “  Py3L  “  Pn  A36  +  F23  A37  +  A38 


A36  -  cos  (i|»  +  a  +  83)  +  y^4  sin  (t|»  +  a  +  83) 


(D-177 ) 


(D-178) 


A37  *  p  s2  sin  (b2  +  02 )  -  cos  (b2  +  63) 


(D-179) 


A38  “  "  Ny  m3 


(D-180) 


Z-Component  of  Force  Equation 

Equation  D-166  cannot  be  further  simplified. 


X~Component  of  Moment  Equation 


Equation  D-170  becomes 

pL  F.  +  L  F  ,  +  p  L  F  +  L  F 

u  x3u  u  y3u  L  x3L  L  y3L 


where 


-  a39  +  a4Q  * 


A39  “  *xs  ux  +  uy  uz  (**s  ”  *ys) 


CD-181) 


(D-182) 
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(D-ia3) 


^40  "  “y 


Y-Component  of  Moment  Equation 

Equation  D— 171  becomes: 

-L  F,  +  L  |iF,  ■  L  F  „  +  L  |i  F  ,, 

u  x3u  u  y3u  l  x3L  l  y3L 

“  A41  +  A42  * 

where 

A/,  ti)  +  w  uj  ( I  —  X  ) 

41  Ays  y  x  z  ^  xs  Z8' 


A/.,  *•  -  I  to 

42  *28  "x 


Z-Component  of  Moment  Equation 

For  present  purposes  equation  D-172  remains  as  it  is. 


(D-184) 

(D-185) 

(D-186) 


Solution  of  Escape  Wheel  Pivot  Forces.  To  derive  expressions  for  the 
escape  wheel  pivot  Forces^  equations  D-173 ,  D-177,  D-181,  and  D-184  must  be 
solved  simultaneously.  To  obtain  the  same  general  form  as  in  equation  D-67, 
equations  D-173  and  D-177  are  multiplied  by  (-1).  The  resulting  form  may  then 
use  the  solution  to  equation  D-67.  Note  that  Ajj  in  equation  D-67  is  now 
replaced  by  u.  Then, 


- 

- 

•  - 

-1 

U 

1 

-y 

Fx3u 

Bsl 

-11 

-1 

11 

1 

Fy3u 

Bs2 

u  L 
v  u 

Lu 

M. 

fx3L 

bs3 

-Lu 

»*  Lu 

•*  K 

py3L 

Bs4 

where  now: 


Bsl  “  ”  lpn  a33  +  f23  a34  +  A35 i 
BS2  ■  “  lpn  a36  +  F23  a37  +  A38J 
®s3  "  a39  +  a40  ♦ 

®s4  “  A41  +  A42  $ 


(D-187) 


(D-188) 

(D-189) 

(D-190) 

(D-191) 


Evaluation  of  the  Coefficient  Determinant  D 


The  solution  for  the  coefficient  determinant  D  of  equation  D-187  is 
identical  to  equation  D-72.  With  Ajj  now  being  equal  to  u,  the  following 
parallel  to  equation  D-75  is  obtained: 

0  -  [(Lu  +  LL)(l  +  w2)]2  (n-192) 

Evaluation  of  Pivot  Force  F  „ 

x3u 


The  pivot  force  Fx3u  is  obtained  from: 
°F 

f  -  x3u 
x3u  D 


vD-193) 


where 


%3u 


BS1 

U 

1 

-M 

bs2 

-1 

P 

1 

Bs3 

Lu 

11  ll 

IjL 

bs4 

11  Lu 

-Ll 

mLl 

(D-194) 


If  u  is  substituted  for  Aij  in  equation  D-77,  the  identical  form  as  above  is 
obtained  and  the  solution  of  equation  D-80  can  be  adapted, 


D  -  (1  +  U2KLu  +  B  -  „  B  +  »  I  -  B  ]  (0-195' 
x3u 

Now  equations  D-188  to  D-191  are  substituted  into  the  above  expression  and  the 
coefficients  of  similar  terms  are  collected.  In  order  to  get  consei native  pivot 
and  pivot  friction  forces,  the  latter  terms  are  made  absolute.  Final Ly,  the 
tllded  force  ?  ^  is  obtained  from  the  appropriate  change  of  equation  D-193: 


where 


x3u 


1 


x3u  ’  °  '(■*„  +  H)U  9-  „2)  [C2‘  +  °22  !>n  +  °23  +  024 


C21  "  III  A35  '  A4l  +  w  (LL  A38  +  A39^ 


(D-197) 


r.-l ...  H  -B 


84 


C22  " 

IH.  (A33  +  u  A36^ 

(D-198) 

C23  " 

1 LI,  (A34  +  w  A37  )  1 

(D-199) 

C24  “ 

|h  a4Q  ”  a42« 

(D-200) 

Evaluation  of  Pivot  Force  F  » 

y3u 


The  pivot  force  Fy3u  is  obtained  from: 


D. 


y3u 


_^3u 


where 


-1 

Bsl 

1 

“P 

-v 

Bs2 

P 

1 

11  Lu 

Bs3 

k. 

-I-,, 

B 

s4 

“lL 

11  Ll 

(D-201) 


(D202) 


Since  the  form  of  the  above,  is  the  name  as  that  of  the  determinant  of  equation  D- 
89,  equation  D-90,  which  represents  the  solution  of  the  latter,  may  be  adapted  as 
follows: 


°F  ,  -  (Lu  +  h,)(l  +  U2)[U  Lj,  V  -  h,  Bb2  +  Bb3  +  u  Bb4)  <D-202a) 

y3u 


Again,  substitute  the  B^  terms  of  equations  D-188  to  D-191,  collect  cVie  coeffi¬ 
cients  of  similar  terms,  and  make  the  result  absolute.  The  tilded  pivot  force 

F  _  then  becomes  parallel  to  equation  D-201: 
y3u 


y3u 


y3u  _ 


(Lu  +  \)[l  +  /) 


[C 


25  +  C26 


Pn  +  C27  F23  + 


C28  ♦) 


(D-203) 


1  **  *  * 


i 


* 


where 

C25  "  I1*!,  Ajg  +  Ajg  +  p  (A^j  -  A^5 )  |  (D-204) 

c26  ■  l\  (*56  '  B  *3j)l  (D-205) 
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JS, 

& 

1 


;N 


c>:- 


i 

!!•’* 


c27  “  Ih,  (^37  "  W  A34 )  I 

C28  "  lA40  +  u  a42  I 
Evaluation  of  Pivot  Force  F 


x3L 


The  pivot  force  Fx3l  is  obtained  from: 
D* 

F 


F 


x3L 

x3L  “  0 


where 


rx3L 


-1 

W 

Bsl 

-w 

-u 

-l 

bs2 

1 

»*  Lu 

^u 

®s3 

lL 

"Lu 

“  Lu 

bs4 

M  \ 

(D-206) 

(D-207) 


(D-208) 


(D-209) 


H.V 


* 


Since  the  form  of  the  above  is  the  same  as  that  of  equation  D-98,  equation  D-100 
may  be  adapted,  i.e., 


°F  ,  ‘  K.  +  H.K1  +  “  )tLa  Bsl  +  “  Lu  B»2  +  “  Bs3  -  B9«)  (D-210) 

x3L 


Again,  the  Bsj_  terms  are  substituted  according  to  equations  D-188  to  D-191  and 
the  requisite  work  obtains  the  tilded  determinant  D-,  »  Then, 

x3L 


'  x3L 


x3L  D 


- - 1 - -  [C  +  C,_  l  +  C-.  F-.  i  C»,  $] 

(Lu  +  LL3(l  +  M2)  29  30  31  23  32  (D-211) 


where 


C29  ly  (A39  “  Lu  A38  ^  "  Lu  A35  “  A41 1 
C30  “  lLu  (^3  +  M  A36^ 


CD-212) 

CD-213) 
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JT-a* 


C31  "  iKi  (a34  +  M  a37^I 
C32  “  ly  A/»0  ”  a42  I 


(D-214) 


(D-215) 


Evaluation  of  Pivot  Force  F 


y3L 


The  pivot  force  Fy3L  t8  obtained  from 


J Qk 


y3L  D 


where 


•y3L 


-1 

M 

1 

Bsl 

“U 

-1 

V 

Bs2 

l‘  Lu 

Lu 

B83 

"Lu 

“Lu 

~lL 

bs4 

(D-216) 


(D-217) 


Since  the  form  of  the  above  is  the  same  as  that  of  equation  D-108,  equation  D-110 
may  be  adapted  to  the  present  situation,  therefore, 

DF  .  -  (LU  +  h.)(l  +  “2)  I-"  *•„  B.!  +  L„  8.2  +  8.3  +  “  8s4}  CD-218) 

y3L 

The  Bs^  terms  are  now  substituted  according  to  equations  D-188  to  D-191,  terms 
are  collected  and  the  tilded  pivot  force  is  defined: 


ilk 


1 


y3L  °  (L  +  L. 1(1  +  „*) 


23  tC33  +  C34  Pn  +  C35  F23  +  C36  ^ 
u  L''-1  ’  H  >  (D-219) 


where 


c33  "  I*1  (A41  +  Lu  A35 )  +  A39  "  Lu  A38l 
c34  "  K  (w  a33  “  a36 ) I 
C35  “  lLu  a34  “  A37 )  I 
C36  -  !A40  u  A42^ 


(D-220) 
CD-221 ) 
(D-222) 


}’  '* 
i:/- 


(D-223) 


rpz^n^rpr^TKSjKrr'triTja^.*'-..^ . 


Substitution  of  Conservative  (Tilded)  Pivot  Forces  Into  the  t-Coraponent 
of  the  Moment  Equation*  Again,  let  tine  sum  of  the  tilded  pivot  forces  be  first 
determined.  Subsequently,  this  expression  is  substituted  into  the  moment  equa¬ 
tion  D-172.  Then, 


Fx3u  +  Fy3u 

■  a43  +  a44 


+  Fx3L  +  F 


Pn  +  A45  F 


y3L 

23  +  A46 


<fr 


where 


LT  -  Ly  + 


A43 


'44 


A45 


A46 


C21 

+  C25  + 

C29  + 

C33 

Uj  (1  + 

'  u  J 

C22 

+  C26  + 

C30  + 

S4 

4  (l 

+  u2) 

C23 

+  C27  + 

°3l  + 

C35 

W  (i  + 

2, 

U  ) 

C24 

+  c28  + 

C32  + 

C36 

lt  (1  + 

2  > 

W  J 

(D-224) 

(D-225) 

(D-226) 

(D-227 ) 

(D-228) 

(D-229) 


Equation  D-224  is  now  substituted  into  equation  D-172.  Further,  Fz3  of  equation 
D-166  is  made  conservative,  i.e., 


Equation  D-172  then  becomes: 

-  (a;  -  »,  VJ  +  f23  I'm  -  » »2  ("2  '  *2)1 

-M  Pf3  A47  -  u  P3  [A43  +  A44  Pn  +  A45  F23  +  A46  $] 

“  lz»  “a  +  1*8  ♦ 


(D-230) 


fD-231) 
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The  above  expression  must  now  be  solved  for  Pn.  Before  this  Is  possible, 
consider  the  sign  of  the  friction  aoment  component: 


-w  P3  A46  I 


(D-232) 


Since  e  reversal  of  gear  train  motion  after  Impact  will  again  be  expressed  by 
letting  u  become  negative,  as  described  originally  in  appendix  E  of  reference  4, 
equation  D-232  is  modified  to  read: 


‘  “  °3  *46 


Ji 


(D-233) 


In  this  way,  the  sign  of  p  alone  decides  the  direction  of  this  friction  torque 
component.  Pn  is  then  obtained  from  equation  D-172: 


P„  I-  *1  +  “1  “l\  ■  11  P3  *44) 

+  F23  trb3  M  *2  (d2  ■  *2^  “  “  P3  *45 1 

* 

l2 

~  p  p3  A46  |-j  “  y  tpf3  A47  +  p3  A43^ 

X..  ^  +  ^zs  uz 

(D-234) 

Then 

•• 

Izs  ♦  +  A48  *  +  F23  A49  +  A50 

Pn  A 

(D-235) 

*51 

where 

p  p3  A46 

A48  *  ... 

Ul 

(D-236) 

A49  "  v  [82  (d 2  a3)  +  p3  a45 ]  ”  rb3 

(D-237) 

A50  "  Izs  K  +  u  fpf3  A47  +  p3  A43  J 

(D-238) 

A5i  *  Bi  V4  “  Ai  -  v  p3  a44 

(D-239) 

t;; 

C 


v 
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Combined  Entrance  Coupled  Motion  Differential  Equation.  Equations  D-140 
and  D-235  are  now  set  equal  to  each  other”  This  furnishes  the  following  combined 
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coupled  motion  differential  equation  for  the  escapement  under  entrance  condi¬ 
tions: 

1*51  l?R  0  -  *29  '..]  ♦  +  t*51  {*32  1,2  +  *9*  V)  -  *29  *48 1  »* 

'  +  A51  A31  U  *  “  F23  A29  A49  +  ^9  A50  "  A5l  +  A30^ 


+  A,,  ra  r  (K  sin  0  -  K  cos  0) 
51  p  cp  v  x  y  ' 


(D-240) 


i/ 
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Pallet  and  Escapement  Wheel  in  Exit  Coupled  Motion 


Pallet  Equations.  The  free  body  diagram  of  the  pallet  for  exit  coupled 
motion  ia  given  by  figures  D-5a  and  D-5b.  Now 


Vn  m  ~  pn  *n 


(0-241) 


This  sign  change  will  be  reflected  both  in  the  force  and  in  the  moment  expres¬ 
sions.  The  following  shows  the  relevant  changes  in  equations  D-23  and  0-140: 

Changes  in  Force  Equations  of  Pallet 

Equation  D-24  is  modified  to  accommodate  equation  0-241.  The  associated 
friction  forces  have  their  directions  determined  by  the  signura  functions  s^  and 
s5,  as  before.  Therefore,  equation  D-24  is  changed  in  its  first  term  only: 


-P  n  +  p, s,  P  n  + 
n  n  1  4  n  t 


(D-242) 


With  the  unit  vectors  of  equations  D-25  and  D-26,  the  x' -force  equation  D-27  is 
modified  to: 

P  sin  ( ip  +  a)  +  p  s.  P  cos  (<|>  +  a)  -  F  ,  -  p,  8,  F  , 

n  1  4  n  x'u  1  5  y  u 


+  Fx’L  + 


(0-243) 


The  terms  in  the  y' -expression  D-28  are  changed  as  follows: 


-  P  cos  ( \p  +  a)  +  p.s.  P  sin  (v|>  +  a)  -  F  ,  -  p,  s.  F  , 

n  14n  y'u  1  5  x'u 


r 

ir. 

r. 


t  \ 

P 

u. 


+  F  + 
y'L 


(D-244) 


The  expression  for  Pz*  remains  as  given  by  equation  0-29. 

Changes  in  Moment  Equation  of  Pallet 

The  form  of  P_,  according  to  equation  D-241,  also  reflects  itself  in  the 
expression  for  TL  (eq  0-61).  Therefore,  for  the  exit  case 
UP 
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Pallet  rotates  cv.  T  ,  *»  ^.a-  p-F  .  (Thrust  friction 

The  actual  directions  of  the  j  Z  torque) 


Figure  D-5b.  Pallet  in  exit  coupled  motion.  Normal  forces,  friction  forces,  and 
thrust  friction  torque  acting  on  pallet  pivots. 


."■  rrrrr.r^-r^.?-?  —^.rvA;;?;..’.  „.,-  x::  n- 


M, 


‘0  •  -  »!  pn  k’  -  **1*4  °1  pn  k'  •  • 
P 


(D-245) 


Simplification  of  Force  and  Moment  Equations . 


x* -Force  Component 


Due  to  the  change  shown  in  equation  D-243,  the  parameter  AJ6  in  equation  D- 
48  must  be  changed  to  become: 


AA^  -  -[p^  cos  (*}*  +  <*)  +  8in  +  <*)] 


(D-246) 


y1 -Force  Component 


Similarly,  because  of  the  change  in  equation  D-244 ,  the  parameter  In 

equation  D~55  must  be  changed  to: 


AA 


2i  "  ”  tui84  8in  ^  +  "  co8  (1»  +  a)  ] 


(D-247) 


z* -Force  Component 


The  z'-force  component  remains  as  that  given  by  equation  D-61,  as  stated 
earlier. 


x1-  and  y* -Moment  Component  Equations 


The  x'-  and  y’ -moment  component  equations  remain  unchanged  from  equations  D- 
64  and  D-65,  respectively,  since  they  do  not  contain  Pn. 


z* -Moment  Component  Equation 


Because  of  the  changes  shown  in  equation  D-245,  the  z* -moment  component 
expression  D-66  must  now  be  modified  to  read: 


-  P  (DJ  +  u  s  C!  )  -  p  A  F  ,  -  p  A  F  -  p  An  F  , 
n  v  1  14  1 '  f  11  z'  p  11  y'u  p  11  x'u 


—  p  A.  F.  -p  A..F. 
p  11  y'L  p  11  x'L 


■  -  w„  r„„  (K*  8in  0  -  K  cos  0)  +  A9  +  Aio  + 


p  cp  v  x 


(D-248) 


Solution  of  Pallet  Pivot  Forces.  The  solution  for  the  pivot  forces 
Fx»„,  Fy»u,  Fxil,  and  Fyi^  is  identical'  to  that  shown  for  the  entrance  coupled 
motiun,  keeping  in  mind  that  now  the  parameters  AA.^  and  AA21  are  used  Instead  of 
and  A2i*  Equation  D-117  must  subsequently  be  cnanged  to: 


’  V  XT.  '-T.Sr.W'.p 
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Is 


H 


Fx’u  +  Fy’u  +  Fx’L  +  Fy'L  " 


♦  2 


A24  +  ♦  *25  +  ♦ *26  +  ♦  ^7  +  Pn  ^ 


28 


(D-249) 


A2/t  to  A27  remain  the  same;  ao  does  AA2g  as  long  as  it  Is  realized  that  it  con¬ 
tains  AA^g  and  AA2^« 


Substitution  of  Pallet  Pivot  Forces  into  z'-Moment  Component  Equation: 
Determination  of  Pn.  Because  of  the  changes  in  equation  D-248,  and  using  the 
same  reasoning  as  employed  for  equations  D-123  to  D-128,  equation  D-128  becomes 
for  exit  coupled  motion: 


Pn  M29  "  A30  "  A31  *  “  A32  * 


•2 


■  IPR  *  +  Aq  -  mp  rcp  (**  8ln  8  -  Ky  cos  3) 


(D-250) 


where 


M29  ' 


-  [o'  +  c;  Vi  ♦  P  aa28] 


(D-251) 


Finally,  parallel  to  equation  D-137,  the  contact  force  Pn  becomes: 


•  2 


XPR  *  +  A9  +  A30  +  A31  *  +  A32  *  "  wp  rcP  K  8in  0  ~  Ky  co8  8) 


AA 


29 


(D-252) 


, 

If  this  expression  is  rewritten  in  terras  of  the  escape  wheel  variables  and  $ 

th^  following  equation  which  is  similar  to  equation  D-140  is  obtained: 


p„  ‘  53!^  (lro  u  ♦  +  (*”  1,2  +  I»«  v)  ♦ 


32 


PR 


+  A31  U  ♦  +  S  +  ^0  '  mp  rcp  ^Kx  8in  8  "  Ky  cos  8)J 


(D-253) 


Escape  Wheel  Equations.  The  free  body  diagram  of  the  escape  wheel  in 
exit  coupled  motion  is  given  in  figures  D-6a  and  D-6b.  The  change  in  the  contact 
force  Pn  must  again  be  accounted  for. 
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Changes  in  Force  Equations  of  Escape  Wheel 

Equation  D-161  must  be  modified  to  read: 


F»  \  -  V*  Pn  "t  +  F23  °21  + . 


(B-254) 


With  the  unit  vectors  of  equations  D-143  and  D-144,  the  X-force  component  equa- 
tion  (changed  from  equation  D-164)  becomes: 

P^  sin  (t|»  +  a+  ftj  )  +  Pn  cos  (t|i  +  a  +  ) 

+  F23  .in  (a,  +  e2 )  +  «  .2  r23  CO.  («J  +  % )  -  Fx3u 


+  U  F  +  F  -  11  F  „  “N  m, 
y3u  x3L  y3L  x  3 

The  Y-force  component  is  changed  from  equation  D-165  to  read: 

-  P^  cos  (i|>+  0+  ^)+  pis4  Pn  sin  (i|i  a  +  03) 

-  F23  “*  («!  +  ^  )  +  “  “2  P23  (%  +  <*>  )  -  Py3u 

-  "  Fx3u  +  Fy3T.  +  “  Fx3l  *  Ny  "3 

The  Z-force  component  remains  as  In  equation  D-166,  i.e., 

Fz3  *  N3  m3 


(D-255) 


(D-256) 


(D-257) 


Changes  In  Moment  Equations  of  Escape  Wheel 

The  moment  equation  D-167  for  the  escape  wheel  and  pinion  no.  3  must  also 
reflect  the  change  in  P  .  The  left-hand  side  of  the  above  expression,  as  given 
by  equation  D-168  must  be  modified,  because  now  the  cross  product 


A,  n.  x  P  n  -  P  A!  k 
It  n  n  n  1 


(D-258) 


This  results  in  the  following  change  to  equation  D-168 


P„  (*i  +  Bi  V*>  S  +  rb3  F23  1  * 


(D-259) 


The  right-hand  side  of  equation  D-167  remains  unchanged.  The  resulting  X  and  Y 
moment  component  expressions,  i.e.,  equations  D-170  and  D-171,  respectively,  are 
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not  influenced  by  the  above  change.  The  Z-moment  component  expression  D-172  must 
now  read: 


P 

n 


(A1 


+  B! 


V< 


)  +  F 


23 


w82  (d2  “  a2^“ 


(D-260) 


Simplification  of  New  Force  and  Moment  Equations  of  Escape  Wheel. 


X-Force  Component 

Due  to  the  change  in  equation  D-255,  the  parameter  A33  in  equation  D-173 
must  be  changed  to 


AA33  ■  PjS^  cos  ( i|»  +  a  +  B3  )  +  sin  (tj)  +  a  + 

Y-Force  Component 

Similarly,  because  of  the  change  in  equation  D-256,  the  parameter  A-j^  in 
equation  D-177  now  becomes: 


AA36  -  sin  (\I)  +  a+  63  )  -  cos  (<|>  +  a  +  ^  )  (D-262) 

Z-Force  Component 

The  Z-force  component  remains  presently  as  given  by  equation  D-257. 


As  stated  earlier,  the  X-  and  Y-components  of  the  moment  expres¬ 
sions  for  the  escape  wheel  need  not  be  changed.  They  are  used  in  their  final 
form  as  given  by  equations  D-181  and  D-L84,  respectively.  Therefore,  the  X- 
component  of  the  moment  equation  is  given  by: 


V  L  F  +  L  F  +  p  L.  F  +  L_  F 
u  x3u  u  y3u  L  x3L  L  y3L 


A39  +  A40  ♦ 


(D-263) 


The  Y-component  of  the  moment  equation  is 


-  L  F  ,  +  L  pF  ,  -  LF0I+L  p  F  „ 

u  x3u  u  y3u  L,  x3L  1.  y3L 


-  a41  +  a42  * 


(D-2f>4) 
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Solution  of  Escape  Wheal  Pivot  Forces  for  Exit  Coupled  Motion.  Since 
only  the  parameters  AA33  and  AAjg  differ  in  the  set  of  simultaneous  equations  D- 
173,  D— 177 ,  D-263,  and  t>-264  from  those  used  in  the  solution  for  the  pivot  forces 
in  entrance  coupled  motion,  the  latter  is  adapted  to  the  present  situation. 
Then,  according  to  equation  D-196 


x3u 


(Lu  +  Ljfl  +  p  ) 


2  \  fC21  +  CC22  Pn  +  C23  F23  +  C24  ^  (D-265a) 


where,  now 


CC22 

■  is.  (A 

and , 

as  hefore 

C21  “ 

Is.  A35 

C23  - 

IS  1*3* 

C24  - 

1  p  A40 

According  to  equation  D 

<v 

F  1 
y3u 

(l  +  ■ 

A  u 

where 

now. 

CC26 

-  |ll  (», 

and , 

as  before 

C25  * 

IS,  A38  ' 

C27  “ 

IS  lA37 

F'28  ' 

lA40  +  “ 

36 


(D-265b) 


37 


r;  tC25  +  CC26  S  +  C27  F23  +  °28  ♦)  (°-266> 


33 


'34 


(D-267 ) 
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According  to  equation  U-212: 


F 


x3L 


(Lu  +  1^(1  +  y2 ) 


[C29  +  CC™  Pn  +  C 


30 


J31  23 


+  c32  l] 


(D-268) 


where  now, 


cc30  "  I  Tjn  (^33  +  u  M36^l 
and  as  before 


(D-269) 


C29  “ 

'u  ^A39  “  Lu  A38  ^  "  Lu  A35  “  A41 1 

C31  " 

lLu  ^A34  +  u  A37^ 

C32  = 

> 

■£> 

o 

1 

Ni 

According  to  equation  0-219: 

where 

F 

y3L 

now, 

(  W,  2  a  tC33  +  ct34  Pn  +  C35  P23  +  C36  ^ 

+  vl1  +  v  j 

(0-270) 

and , 

CC34 
as  before 

-  |Lu  („  AA13  -  M,6)| 

(0-271) 

C33  " 

|p  (A/4l  +  Lu  a35)  +  A39  _  \  A38 1 

C35  " 

•Lu  ^  A34  "  A37^ 

C36  " 

!A40  +  w  A42 l 

Substitution  of  Conservative  (Tilded)  Pivot  Forces  into  Z-Moment  Expres- 

sion . 

The  sum 

of  the  tilded  pivot  forces  is  identical  in  form  to 

equation  D- 

224.  Therefore,  with  equations  0-265,  D-266,  D-268,  and  D-270,  the  following  is 
obtained: 
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(D-272) 


Fx3u  +  Fy3u  +  Px3L  +  Py3L 


-  a43  +  aa44  Pn  +  a45  f23  +  a46  4» 


where  now, 


M44  + 


and,  as  before 


cc22  +  cc26  +  cc30  +  cc34 


Uj,  (l  +  u2) 


V43  " 


C21  *  C25  *  C29  +  C33 

~L"(i  +  u2) 


A  .  C23  +  C27  +  C3i  +  C35 

A45  “ - >1 - 


W  (1  +  W  ) 


A, 


46 


C24  +  C28  *  C32  *  C36 
Ur"_(l  +  M  ) 


(D-273) 


Substitution  of  equations  D-257  and  D-272  into  equation  D-260  furnishes  the 
complete  Z-component  of  the  escape  wheel  moment  expression  for  exit  coupled 
motion: 


Pn  (*1  +  SVJ  +  F23  (rM  '  u  *2  •  *2)1  '  u  pf 3  *47 

-  u  o3  [*43  +  aa44  pn  +  a45  r23  +  a46  ♦]  -  i„  \  »  (0.274) 


Using  the  same  reasoning  as  given  in  connection  with  equations  D-232  and  D-233, 
equation  D-274  now  solved  for  p  Therefore, 


Pn  lAl  +  B1  m184  "  W  P3  M44^ 


+  F23  lrb3  ‘  P  S2  (d2  "  a2}  "  W  P3  A45 


“  w  p3  A46  -jij*  ~  u  ipf3  A47 


p3  A43  -1 


zs 


4>  +  I 


zs 


(D-275) 
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and,  similar  to  equation  TV-235 

?2 


P  - 

n 


lzs  ♦  +  A48  41  +  F23  A49  +  A50 


AA 


51 


where  now 


M51  "  A1  +  B1  uls4  M  p3  M44 


while  as  before. 


U  p,  A, 


A48  * 


p3  A46 

♦I 


A49  “  y  [«2  (<*2  -  a2  )  +  P3  A45  ]  "  rb3 


A50  "  1zb  K  +  U  ^pf3  A47  +  p3  A43  J 


(D-276) 


(D— 277  ) 


Combined  Exit  Coupled  Motion  Differential  Equation.  Equations  D-253  and 
D-276  are  now  set  equal  to  each  other  in  order  to  obtain  the  combined  coupled 
motion  differential  equation  of  the  escapement  under  exit  conditions: 

[M51  *  +  Kl  ^32  U?'  +  V  V  ) 

“  ^29  A48  ^  ^  +  M51  A3i  U  $ 

"  F23  M29  A49  +  M29  So  "  M51  ^9  +  A30  ^ 

+  AA  m  r  (K  sin  8  -  K  cos  8)  (D-278) 

51  p  cp  x  y 

The  above  expression  has  the  same  form  as  equation  D-240,  and  the  difference 
between  entrance  and  exit  coupled  depends  on  the  value  of  the  signum  function  S7 
which  is  introduced  in  the  next  section. 


Common  Expressions  for  Entrance-  and  Exit-Coupled  Motion 


It  is  possible  to  obtain  common  expressions  for  all  the  A's  and  C's 
(i.e.,  AA's  and  CC’s)  associated  with  entrance-and  exit-coupled  motion,  if  the 
signum  functions  ay  are  introduced,  where 

87  *  positive  for  entrance-coupled  motion 
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"P" 


"TT" 


By  m  negative  for  exit-coupled  motion 
With  the  above,  equations  D-54  and  D-246  are  satisfied,  if 

A16  "  "  [yls4  cos  ^  +  “  ®7  8itl  +  «)] 

Equations  D-60  and  D-247  are  satisfied,  if 

A21  “  “  [87  cos  +  <0  +  ein  (ty  +  «)] 
Equations  D-130  and  D-251  are  satisfied,  if 

A29  “  s7  Di  “  C1  wi®4  “  pp  wls5  A28 
Equations  D-174  and  D-261  are  satisfied,  if 

A33  *  pi84  cos  ('|i  +  a  +  83)  -  87  sin  (ty  +  a  +  83) 
Equations  D-178  and  D-262  are  satisfied,  if 

A36  "  pls4  8*n  ('ll  +  ct  +  83 )  +  By  cos  (i|»  +  a  +  83 ) 
Finally,  equations  D-239  and  D-277  are  satisfied,  if 

A51  ■  BJ  ul84  -  s?  -  u  P3  a44 
Dynamics  of  Rotor 


(D-279a) 


(D-279b) 


(D-280) 


(D-281) 


(D-282) 


(D-283) 


Before  the  force  and  moment  equations  for  the  rotor  can  be  considered,  it  is 
first  necessary  to  obtain  expressions  for  the  absolute  accelerations  of  the  rotor 
pivot  Oj  and  the  rotor  center  of  mass  .  A  top  view  of  the  rotor  in  the  mechan¬ 
ism  plane  is  shown  in  figure  D-7  (also  rig.  A-3). 
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Absolute  Acceleration  of  Rotor  Pivot  0, 


where 


while 


The  absolute  acceleration  of  the  rotor  pivot  0}  is  given  by 
^0^  ground  "  A0  ^  /C  +  AC/ ground 


(D-284) 


^C/ground  “  given  by  equation  C-4  of  appendix  C  in  the  projectile  fixed 
8  X-Y  system, 


VQ  ■  u  x  (to  x  +  w  x 


(D-285) 


After  substituting 

flj  -«j  I  (D-286) 

and  equations  A-l  and  A-5  for  ui  and  u»  respectively,  the  following  is  obtained: 


Vc-L*1  + v  +  L«k 


(D-287) 


where 


“,2>  «1 

(D-288) 

(D-289) 

u  )  <R. 
y;  i 

(D-290) 

Together  with  equations  D-187  and  C-4,  the  following  is  obtained  for  equation  D- 
284: 


where 


Art  ,  -0i  +  0j+0k 

0^ /ground  x  y  z 


°x  "  Gx  +  W 
Oy  *  ®y  + 


(D-291) 


(D-292) 

(D-293) 
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(D— 294 ) 


Absolute  Acceleration  of  the  Rotor  Center  of  Mass  Cj 

To  determine  the  absolute  acceleration  of  the  rotor  center  of  mass  in 
the  X-Y-Z  system,  it  is  first  necessary  to  find  Ac  ^  ,  the  acceleration  of  the 

rotor  center  of  mass  with  respect  to  the  rotor  pivot  0^ ,  in  the  E  -  a  -  r 
system  (fig.  B-7).  Subsequently,  this  expression  is  transformed  into1  the1X-Y-Z 
system  and  added  to  the  absolute  acceleration  of  point  0j  as  given  by  equation  D- 
291.  Therefore, 


^/ground  "  +  ^/ground 


The  term  A_  ,  is  obtained  from: 

Cl/01 

fyo,  -  \  *  *  Jcl  )  +  \  *  5 


(D-295) 


cl 


where 


r  ,■  r  ,  n. 
cl  cl  C, 


(D-296) 


(D-297 ) 


The  terms  and  are  taken  from  equations  A-35  and  A-39,  respectively. 
When  all  operations  are  performed,  equation  D-296  becomes: 


V, ' ' 

rcl 

[(02  + 
nl 

2  1  - 

uc  J  nr 
n  t 

• 

+ 

rcl 

l\ 

0)  +  U) 

"l  V 

+ 

rcl 

l\ 

0),  -  0)  1 
h  V 

r\ 

(D-298) 

With  the  help  of  equations  A-28  and  A-31b  substitute  for  the  above  body-fixed 
unit  vectors;  i.e., 


n_  ■  cos  i  1  +  sin  y  j 

n 


n  ■  -  sin  y  i  +  cos  y  J 


(D-299) 

(D-300) 
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(D-301) 


where 


n_  •  it 
?1 


Y  "  he  +  N31  h 


(D-302) 


This  results  In 


Cj/Oj  ■  '  rcl  +  “^)  =°s  If  +  (“?1  +  “Cl)  sl”  f)  1 


2  2  # 

-  *cl  [(u  +  0>c  )  sin  y  -  (w^  «  +  )  cos  y]  j 


+  r  .  fu-  a).  -  i  lk 

cl  I  et  5!  n,' 


(D-303) 


Finally,  equations  A-36  to  A-38  and  A-40  to  A-42  are  used  to  express  the  angular 
quantities: 

_  2  9  2 

Aq  /0  -  -  rcl  {[a»y  cos  Y  -  wx  uy  sin  Y  +  (u^  +  N^)  cos  y 


•  **  ** 

+  (wz  +  N31  <}> )  sin  y]  i 

2  *2 
+  fa)  sin  y  -  u>  u)  cos  y  +  (u  +  N_,  i)  sin  y 
1  x  x  y  v  z  31  T' 

•  • 

-  (ms  +  N31  ♦)  cos  y]  j 

-  [(a)  cos  y  +  u  sin  y)(a)  +  2  l  N».  )  +  a )  sin  y 
x  y  z  o  i  x 


-  U)  COB  y]  k} 


(D-304) 


The  total  acceleration  A„  ,  .  then  becomes  according  to  equation  D- 

295  with  equation  D-291:  Cj/ground 
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/ground  ‘  '  ‘  rcl  l“y2  coa  ^  ‘  “y  8ln  1  +  +  N31  *)2  c°8  1 

+  (to  +  N31  $)  sin  y]  +  Ox}  I 

2  t  2 

+  {-  rcl  [«*>x  sin  y  -  «xw  cos  y  +  (uz  +  N31  ♦)  sin  y 

-  (»x  +  N31  4>)  cos  y]  +  Oy}  J 
+  {rcl  [(ux  cos  y  +  toy  sin  y)(to2  +  2  $) 

(D-305) 

+  to  sin  y  -  to  cos  yj  +  0  }  k 
x  y  z 

Force  Equations  for  the  Rotor 


A  top  view  of  the  rotor  together  with  the  mechanism  plane  is  shown  in 
figure  D-8a.  It  indicates  all  required  geometry,  the  base  circle  radius  Rbl,  and 
the  contact  force  F  as  well  as  the  associated  friction  force  F^j  (for  details, 
see  refs  1,  2,  and  4).  A  free  body  diagram  of  the  rotor  pivot  with  all  normal 
and  friction  pivot  forces  is  shown  in  figure  D-8b. 

The  force  equation  for  the  rotor  is  given  by: 


£  F  -  ml  ^/ground 

where  ^n/groun<j  *8  f?iven  by  equation  D-305.  Therefore 


“^12  n12  +  ^  8i  nNi9  +  k  -  Fvlii  i  -  F„,„  j 


'1  1  12  N12  T  rzl 


xlu 


ylu 


w  Fxlu  J  +  y  Fylu  1  +  FxlL  1  +  FylL  J  +  w  FxlL  3 


y  FylL  *  *  ml  ^Cj/ground 

According  to  equation  E— 127  and  E-128  of  reference  1 


12 


-  sin  (B1  -  0j_ )  i  +  cos  ($1  -  Qj)  j 


"N12  “  ”co8  (®i  “  9i )  ^  '  sin  (Sj  -  )  j 


(D-306) 


(D-307 ) 


(D-308) 

(D-309) 
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Rotor  rotates  ccw 


Rotor  rotation  is  ccw 


(Thrust  friction 
torque ) 


(L  and  L  represent  one- 
ha'ff  the  distance  between 
the  pivot  bearings  in  upper 
and  lower  plates.) 


Figure  D-8b. 


Rotor  and  gear  no.  1.  Normal,  forces,  friction  forces,  and  thrust 
friction  torque  acting  on  rotor  pivots. 


Tha  above  unit  vactors  ara  now  substituted  Into  aquation  D-294.  Subsequently, 
the  component  expressions  of  this  equation  are  written  with  the  help  of  equation 
D-305i 

X-Component  of  Rotor  Force  Equation 


*12  *ln  <*l  *  8[)  -  “  *1  *12  “•  (S1  -  #l)  '  *,i„  +  »  *ylu 

+  **11.  "  "  *,1L  ’  *1  I'  rcl  t“y2  c°*  * 

^2 


-  wv  uv  sin  Y  +  («„  +  N31  cos  Y 


x  y 


+  (»,  +  H3l  ♦)  sin  y  }  +  Oj 


(D-310) 


Y-Component  of  Rotor  Force  Equation 


-  f12  coa  (Bj  “  0j)  '  V  f12  8in 

+  *,11  +  "  **1L  '  "l  !-  rcl  t“*2  8l"  * 

-  u>x  w  cos  y  +  (,<*>z  +  N31  <t»}  sin  y 

“  (w,  +  N31  4>)  cos  y  }  +  0y] 


-9  1  -  F 
1  J  ylu 


U  F 


xlu 


CD-311) 


Z-Component  of  Rotor  Force  Equation 


Fzl  “  ml  ircl  [K  cos  y  +  “y  8in  y)(»b  +  2N31  4) 


+  (o  sin  y  -  «.  cos  y]  +  0  } 

X  y  Z 


(D-312) 


Factors  Entering  Into  Moment  Equations  for  the  Rotor 


The  moment  equation  for  the  rotor  must  be  written  with  respect  to  the 
accelerated  pivot  point  Oj«  (This  is  similar  to  the  manner  in  which  the  pallet 
moment  expression  D-30  was  written  with  respect  to  point  Op.)  Therefore, 

**>1  ■  “  ^/ground  x  ml  rcl  (c08  Y  1  +  sin  y  J)  +  A0l  (D-313) 
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where 


/ground 


sum  of  external  moments  about  point  Oj .  it  is  assumed 
that  0}  lies  in  the  plane  of  the  rotor  center  of  mass,  and 
that  Fj2  and  M|  Fjj  a^80  *n  t^1-8  P^ane« 

Absolute  acceleration  of  point  Oj  (eq  0-291) 


■  Time  rate  of  change  of  angular  momentum  of  rotor  with 
U1  respect  to  point  0^.  It  is  obtained  by  adapting  equation 

B-4  of  appendix  B  to  the  -  n  -  r  system.  The  appro¬ 
priate  angular  velocity  and  acceleration  components  are 
given  by  equations  A-35  and  A-39,  respectively.  The 
transformation  into  the  X-Y-Z  system  is  accomplished  with 
the  help  of  the  unit  vector  expressions  of  equations  D-299 
to  D-301. 


Determination  of  Hq^ 


The  moments  due  to  the  gear  contact  forces  F?.  and  p  s.  F?.  are  identical  to 
those  of  equation  B-129  of  reference  1.  The  moments  due  to  the  various  bearing 
forces  may  be  adapted  from  equation  D-32.  (Since  the  rotor  always  has  counter¬ 
clockwise  rotation,  let  8j  *  +  1  in  equation  D-32,  change  p^  to  p,  and  adjust  the 
subscripts  from  the  primed  to  the  unprimed  coordinate  system.) 


-fL  F,  +  pL  F,  +  L,  F  +  |i  L,  F  ,,  1  I 
Oj  L  u  ylu  M  u  xlu  L  ylL  L  xlL J 

+  ["  \  Fylu  -  Lu  F*lu  +  11  \  FyU  ‘  \  FxlL  1  * 

+  [-  xM  f12  +  “Vi  Fi2  -  u  »fl  Fs 


-  p,  p  F 


ylu 


-  ”l  "  F*lu  -  P1  “  F 


ylL 


-  p,  «  lIlt]  I 


(D-314) 


Determination  of  Right  Hand  Side  of  Equation  D-313 

With  the  help  of  equation  D-291,  the  following  is  obtained  for  the  right 
hand  side  of  equation  D-313 


-  (°x  *  +  °y  J  +  °2  *  )  x  rai  rci  (  c09  Y  1  +  sin  Y  J  ) 

■  m.  r  ,  [0  sin  y  i  -  0  cos  y  j  -  (o  sin  y  -  0  cos  y)  h]  (D-315) 
x  cx  z  z  x  y 
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Determination  of  Tina  Rate  of  Change  of  Angular  Momentum  with  Respect  to  Rotor 
Pivot  Oj 

To  obtain  an  expression  for  fl  in  the  r  -  n  -  r  rotor-fixed  system, 

”j  111 

equation  B-4  is  first  adapted  from  the  X-Y-Z  system*  Subsequently,  the  angular 
velocity  and  acceleration  of  equations  A-35  and  A-39  are  substituted  as  follows: 


U)_  »  U>  COS  Y  +  0)  Sin  Y 

x  y 


(ii  *  -  o>  sin  y  +  oi  cos  y 
Hj  x  y 

\  '  +  N3I  1 


(D-316) 

(D-317) 

(D-318) 


«r  ■  w  cos  y  “  u  L  |  sin  y  +  w  sin  y  +  w  N_.  $  cos  y  (D-319) 

*»  2  x  x  j  i  y  y  j  l 

ui  ■  -  w  sin  y  ~  oi  N*,  $  cos  y  +  w  cos  y  “  oi  N,,  J  sin  y  (D-320) 

x  x  31  y  y  31  r 


oi  ■  oi  +  N«,  <b 
^ 


CD-32 1 ) 


Finally,  the  body-fixed  unit  vectors  n^  ,  n^  ,  and  n^  are  given  in  the  X-Y-Z 
system  according  to  equations  D-299  to  D-3fcl .  Ws  furnishes 

Hqi  -  {i^  («x  cos  Y  -  N31  |  sin  y  +  sin  Y  +  o>y  N31  $  cos  y) 


+  (-  uix  sin  Y  +  o)y  cos  y)(«£  +  N31  *){l^  -  1^) 


+  I5n1  +  N31  ♦  cos  Y  +  sin  y)  -  (-  sin  y 


-  uix  N31  p  cos  y  +  ov  cos  y  -  u)v  p  sin  y)] 


y  "31 


-  I  [(01  cos  y  +  o)  sin  y)("  w  8in  y  +  oi  cos  y) 
tCj  x  y  x  y 

+  (o»  +  N31  ♦)]  -  I  ((-  to  sin  Y  +  0)  COS  Y) 


**1 


-  (u>  +  N31  i)2]}  (cos  Y  1  +  Sin  Y  j) 


z  '  "31 
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1 »  '  A*  *  S .  ‘  11  '  >» 


^  .  i  'n 


i 

i 

r. 

£ 

i  • 

F: 


k- . 

■ 


+  {lnr)  (-  ux  sin  y  -  wx  H31  |  cos  y  +  wy  cos  y  -  wy  N31  $sin  y) 


+  (ii»x  cos  y  +  wy  sin  y)(«,  +  N31  -  I?5  ) 

+  Inc  [(wx  cos  y  +  <oy  sin  y)(“  «x  sin  y  +  oiy  cos  y)  -  (<*>z  +  N31  ♦  )] 

-  I^n  [(wx  cos  y  -  ujx  N31  0  sin  Y  +  ay  sin  y  +  uy  N31  $  cos  yj 

+  (-  u»x  sin  Y  +  wy  cos  y)(wz  +  «31  $)]  -  [(«# 

•  2  2  „  * 

+  N_.  $)  -  (a>  cos  y  +  u  sin  y)  ] }  (-  sin  y  i  +  cos  y  j) 

Ji  x  y 

+  I1,;  +  »31  <j>)  +  (ux  cos  y  +  aiy  sin  y)(-  aix  sin  y  +  u)y  cos  y) 

x  (ln^  '  *«,>  +  Icel  (1‘  “x  8ln  Y  +  “y  cos  y)K  +  “n  ♦) 

-  (ix  cos  Y  -  U)x  N31  $  sin  y  +  <oy  sin  y  +  uiy  N31  l  cos  y)] 

“  [(“  «x  sin  Y  -  wx  N31  J  cos  y  +  u)y  cos  y  -  uy  N31  l  sin  y) 


+  (ojx  cos  y  +  0)y  sin  y)(u>z  +  N31  ♦  )]  -  I?  [(wx  cos  y 

2  9 

+  U)  sin  y)  -  (-  w  sin  y  +  ui  cos  yf  ])  k 
y  x  y 


(D-322) 


The  components  A  ,  ft  ,  and  ft  must  now  be  determined  from  equation  D-322. 
This  leads  to:  ^lx  ly  lz 


H0,  "  A52  +  A53  ♦  +  A54  *  +  A55  * 

lx 


r .  , 


"■  v  A 


■■■'4 

■:7 
..  -1 

■'  .j 

■ .  i 
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«  cos 


Y  { I ^  COS  Y  +  <^  sin 


+  (I«,  '  ^  “•  \  *in  Y+  “y  c<“  y) 

+  I  [u)z  (ujx  COS  Y  +  W  Sin  y)  +  («>x  sin  Y  “  <*>  cos  y]] 

5nl  Z  X  y 

-  t5  [(i»x  cos  Y  +  uy  sin  Y)(-  \  8in  Y  +  <*>y  cos  y)  +  <*>2  ] 

2  2 

-  I  [(-  <i)x  sin  Y  +  wy  cos  y)  “  \  ]} 

-  sin  Y  {l^  ("  ^  sin  y  +  wy  cos  y) 

+  fl  -  I  )  id  (u  cos  y  +  os  sin  y) 

J  z  1  x  y 

+  I  [(u^  cos  Y  +  <*y  Sin  y)(-  ux  sin  y  +  osy  cos  y)  -  \  ] 
.  [(u)  cos  Y  +  <*>  sin  y)  +  t*>  (~  ol  sin  Y  +  <»L  cos  y)] 


-  Ir_  [u2  -  («  COS  Y  +  uv  sin  y)“  ]} 
1  y 


(0-324) 


-  N3i  H-  “x  8in  Y  +  Wy  C08  +  "  Inn1 3  cos  Y 

+  21  ^  sin  y]  +  [wx  cos  Y  +  <*>y  sin  Y][(lnT^  "  1  ^  +  1  ^  )  sln  Y 


+  21  ^  cos  y]+  2  u»z  [i^  cos  Y+  sin  y] 


(0-325) 


V‘  ".'.V 
» ,  » i 


N3i  ^ne-  co0  Y  +  ^  8in  y] 


N_.  f-  I  cos  Y  +  I  sin  y] 


(0-326) 


(D-327) 


Further , 


V  -  A56  +  A57  *  +  A58  ♦  +  A59  ♦ 

iy 

A56  -  ein  Y  [lpF  (<»>*  r-os  Y  +  <*>v  ain  y) 


+  fl  -I  If-  w  ein  y  +  «*>  cos  y)  u> 

1  nnj  *■  x  y  '  z 

+  Ip  [wz  (wx  COS  Y  +  Hi  sin  y)  “  (-  sin  y  + 


-  I_  [fui  cos  y  +  w  sin  y)(“  o>  sin  y  +  o>  cos 

x  y  x  y 

2  2 

-  I  ff-  a)  sin  y  +  w  cos  y)  -  a)  ]  1 

nc1  u  x  y  ;  z  JJ 

+  cos  y  1 1  f-o)  ain  y  +  o>  cos  y) 

1  nn.  v  X  y  ; 

+  (1K  - )(»x  «*  i +  »y sln  v) », 


+  I  „  ffoi  cos  y  +  u>  sin  y)(“  <“  sin  Y  +  <*>  cos 
1  v  x  '  y  '  'v  x  1  y 


-I.  fa)  cos  Y  +  a>  8in  y  +  (o  (-  u)  sin  y  +  u 
Cnx  1  x  y  z  x  y 

2 


-I.  fa)  ^  -  fa)  cos  Y  +  w  sin  y)  If 
SCj  1  z  v  x  y  J  J  J 


A57  "  N^»’  {[-  uiv  sin  y  +  <»„  cos  y] 


31  11  x  y 


x  [ ( I _ _  +  I  -I  1  sin  Y  -  2  cos  y] 

WL  HHj  j 

+  fa)  cos  Y  +  a)  sin  y1 

i  x  y  J 


x  2  I,  sin  y  +  (l-„  -  I  -  I  1  cos  y] 

Cn  v  f  r  r  r  nn  '  J 


Cn 


5^  nnl 


+  2  a)  fl  sin  Y  -  I 

<*  nr 


z  1  n? 


55, 


cos 


y]\ 


(D-328) 

ay  cos  y)] 

y)  +  wj 


Y)  ~  ig] 

cos  y)] 
(D-329) 


(D-330) 
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(D-331) 


*58  '  "31  [%?,  ,ln  V"  “*  1,1 
*59  ■  '  "31  t1?;,  ,l”  'nc,  “*  yl 


(D-332) 


Finally, 


where 


H.  ■  Atri  +  A,.  $ 
0.  60  61  T 
lz 


1^l  \ + t1™,-  cos  ’'+  “v sln  v) 


x  f-  u  aln  y  +  oj  cos  yl 
v  x  y  * 


(D-333) 


+  1  _  f  f-  (i)  sin  y  +  (d  cos  y!  w  -  u>  cos  y  -  u  sin  yl 
5£  1  *■  x  y  '  z  x  y  J 

+  I?T^  li  8ln  Y  “  “y  cos  Y  “  w,  (u»x  cos  Y  +  (^  sin  y  )  ] 

-  T._  [f  u>  cos  y  +  w  sin  Yf  ~  f-  w  sin  y  +  w  cos  yf  1  (D”334) 
£n.  1  *  x  y  *  *  x  y  J  1 


A,,  -  N  I 

61  31  55, 


(D-335) 


Simplification  of  Force  Equations  for  the  Rotor 


X-Component  of  the  Force  Equation 

Equation  D-310  is  now  rewritten  in  the  following  manner: 


-F  +  uF  +  F  -  uF 
xlu  ylu  xlL  ylL 

"  A62  +  A63  ♦  +  A64  f  +  A65  ♦  +  A66  F12 


(D-336) 
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where 


2 

A62  "  m!  rci  ["  <*y  coe  y  +  (i ^  sin  y 

-  ti)  ^  cos  y  -  ai  sin  y]  +  nt  0 

Z  Z  A  X 

(D-337) 

A63  -  -  2  mj  rcl  o)z  N31  cos  y 

(D-338) 

A64  "  “  mi  rci  N3i  cos  Y 

(D-339) 

A65  "  "  "j  rcl  n31  8in  V 

(D-340) 

A66  “  W  8j  cos  (6,  -  61  )  -  sin  (S1  -  0t  ) 

(D-341) 

Y-Component  of  the  Force  Equation 

Equation  D-311  becomes: 

-  F,lu  -  11  F,lu  +  Fyll  +  “  FxlL 
'  *67  +  A68  ♦  +  *69  +2  +  *70  **  *71  F12 

(D-342) 

where 

2 

A67  “  rai  rcl  ["  “x  8in  Y  +  \  wy  cos  Y 

2  •  , 

-  o)  sin  y  +  w  cos  y  +  m,  0 
z  z  1  1  y 

(D-343) 

A68  “  "  2  mi  rcl  N3i  “z  8in  Y 

(D-344) 

A69  "  -  ml  rcl  N31  8in  Y 

(D-345) 

A70  "  -  rai  rci  N31  COS  y 

(D-346) 

A7l  -  cos  (8j  -  01 )  +  y  Sj  sin  (fy  -  Oj  ) 

(D-347) 

Z-Componant  of  Force  Equation 


Equation  D-312  Is  rewritten  In  Its  tilded  form  directly: 


f,l  ‘  *72  +  *73  * 

where 

a72  *  \mi  rei  [“2  (wx  coa  ay  sin  y) 

+  (ii  sin  y  -  w  coa  yl  +  m,  0  I 
x  y  J  1  z ' 

a73  I2  mi  rcl  N31  K  coa  y  +  ay  sin  y]|  (D-350) 

Simplification  of  Moment  Equations  for  the  Rotor 


(D-348) 


(D-349) 


The  components  of  the  rotor  moment  equations  are  now  written  according 
to  equation  D-313. 

X-Component  of  Moment  Equation 


With  the  help  of  equations  D-314,  D-315,  and  D-323,  the  following  is 

obtained: 


y  L  F  +LF  +  u  L  F  +LF 
u  xlu  u  ylu  M  L  xlL  L  ylL 

"  ral  rcl  °z  8in  T  +  A52  +  A53  ♦  +  A54  ^  +  A55  ♦ 


(D-351 ) 


Y-Component  of  Moment  Equation 

Again  with  the  help  of  equations  D-314  and  D-315,  i.e.,  its  y-factors,  as 
well  as  equation  D-328,  the  following  is  found: 


-L  F  +  u  L  F  +  \i  L  F  ,T  -  L  F  lt 
u  xlu  u  ylu  L  ylL  L  xlL 


"  ml  rcl  °z  COS  A^  +  A<?  *  +  A<«  <>  +  A*o  ♦ 


56 


57 


58 


59 


(D-352) 


Z-Component  of  Moment  Equation 

Again,  using  the  Z-componenta  of  equations  D-314  and  D-315,  together  with 
equation  D-333,  obtained  for  the  Z-coraponent  of  the  moment  expression 
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”  Rbl  F12  +  w  81  al  F12  11  pfl  Fz  v  P1  ^Fxlu  +  Fylu  +  FxlL  +  *ylL  ^ 

i* 

-  -  *1  *ci  [Ox  Bin  Y  -  Oy  cos  y]  +  A60  +  A61  *  (D-353) 


Solution  of  Rotor  Pivot  Forces 

To  obtain  the  rotor  pivot  forces,  equations  D-336,  D-342,  D-351,  and  D- 
352  must  be  solved  simultaneously.  Therefore, 


-F  +  u  F  .  +  F,t  “  |i  P  „  *  B.. 

xlu  ylu  xlL  yll  11 


-  “  Fxlu  '  Fylu  +  "  F«ll  +  Vt.  ■  B12 


“  lu  Fxlu  +  Lu  F,lu  +  “  S.  Fxl,  +  lL  FyU.  “  B13 
'  L„  Fxlu  +  “  Lu  Fylu  '  H  FxlL  +  "  \  V  '  Bl* 


where 


+  A63 

♦  + 

A64 

;2  +  a65 

♦  +  A66 

F12 

+  A68 

l  + 

A69 

^  +  A70 

4>  +  A7i 

F12 

• 

•2  .  . 

cl  °z 

sin 

Y  + 

A52  +  A53 

<f>  +  a5A 

+  A, 

*14  "  -  mx  rci  °z  cos  y  +  A56  +  A57  +  +  A58  *  +  A59 


(D-354) 

(D-355) 

(D-356) 

(D-357) 

(D-358) 

(D-359) 

(D-360) 

(D-361) 


Since  equations  D-354  to  D-357  together  have  the  same  general  form  as  equation  D- 
67  for  the  pallet,  the  forms  of  the  pallet  pivot  force  solutions  for  the  rotor 
pivot  forces  may  be  used.  It  must  be  kept  in  mind  that  for  the  rotor  the  factor 
p  must  be  substituted  for  A^.  Then,  accordinp  to  equation  D-73: 


D1  “  [(Lu  +  LL^  +  ^  ^ 


(D-362) 


Parallel  to  equation  D-80,  the  determinant  Dp^  becomes: 

D,  -  (L  +  1^)(1  +  }  )[-  \  »u  -  »  \  »12  +  “  B!3  '  BU  1  <D-363) 


After  appropriate  substitution  of  equations  D-358  to  D-361,  parallel  to 
equations  D-81  to  D-B7,  the  following  is  obtained  for  the  conservative  rotor 
pivot  force: 


TJ~  ■  7 — 77~~TT  fC37  +  C38  ^  +  °39  *2  *  °40  ..  .. 

1  4  +  “  )  CO-364) 


♦  +  °A, 


where 


C37  “  1“  h.  A62  +  »*  (A52  “  h  A67 )  “  A56 


+ 

fcl 

0  [ 
z  L 

u  sin 

Y  + 

CO  8  y]| 

(D-365) 

C38  * 

1“  ll 

A63 

+  u 

(A53  “ 

A68 )  ‘  A57  1 

(D-366) 

C39  “ 

1-  ll 

A64 

+  u 

(A54  “ 

A69  ^  ‘  A58 1 

(D-367) 

C40  ” 

I'  h. 

A65 

+  u 

(A55  " 

h 

A70  )  "  A59  1 

(D-368) 

C41  “ 

1-  h- 

^A66 

,  +  u 

A7l)l 

(D-369) 

Parallel  to  equation  D-89,  the  determinant  Dp  becomes: 

rylu 

°F  "  (Lu  +  +  LL  B11  “  h  B12  +  ®13  +  W  B14  ^  (D"370) 

ylu 

After  appropriate  substitution  of  equations  D-358  to  D-361,  parallel  to 
equations  D-91  to  D-96,  it  is  found  that: 


*2  + 


'45 


<f>  +  C 


46 


CD-37  1) 


where 


C42  “  1“  ll  a67  +  u  (A56  +  H  a62  )  +  A52 
+  ml  rcl  °z  f8in  y  ”  w  cos  y]| 


(D-372) 
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C43  "  I"  \  a68  +  11  (lL  a63  +  a57  )  +  A53  ^ 

c44  "  I-  4,  a69  +  u  (ll  a64  +  a58)  +  a54| 

C45  "  I-  \  a70  +  »*  (LL  A65  +  A59)  +  A55  S 

c46  "  IH,  (**  a66  ”  a7  1 )  i 


(D-373) 

(D-374) 

(D-375) 

(D-376) 


Parallel  to  equation  D-99,  the  determinant  Dp  becomes: 

xlL 


D 


F  "  (Lu  +  Ll)(1  +  U  )  lLu  Bll  +  W  Lu  B12  +  u  *13  "  *14  J  (D"377) 
xlL  “ 


Again,  equations  D-358  to  D-361  are  substituted  into  the  above.  Then 
proceed  parallel  to  equations  D-101  to  D-106.  Finally, 


Dt 


xlL 


1 


xlTj  -  -D—  -  ,  f,  'TT  tC47  +  C48  +  +  C49  ♦  +  C50  +  +  C51  F1?J 

1  +  «  )  (D-378) 


c47  "  lLu  A62  +  u  (Lu  A67  +  A52  )  “  A56  +  ml  rcl  °z  ^  sin  y  +  cos 

(D-379) 


C48  “  K  A63  +  M  (Lu  A68  +  A53^  ”  A57 I 

Afc&  +  W  (LII  A6Q  +  As^)  "  ASft  I 


J49  -  l“u 


C50  “  lLu  A65  +  u  (Lu  A70  +  A55^  A 


59 


C51  -  lLu  Kb  +  M  A 


71 


(D-380) 

(D-381) 

(D-382) 

(D-383) 


Parallel  to  equation  D-109,  the  determinant  D-.  becomes: 

FylL 


ylL 


(Lu  +  +  y2)  I-  u  Lu  Bll  +  Lu  *12  +  *13  +  w  *14}  (D'384) 
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After  substitution  of  equations  D-358  to  D-361,  proceed  parallel  to 
equation  0-1 11: 


where 


D. 


'  ylL  D 


ilk 


1 


1  4  (1  +  /) 


[Cs,  +  Cs,  i+  C_.  I2  +  C„  ;+  c,,  F.J 


56  12 

(D-385) 


C52  "  !LU  A67  +  y  (A56  ~  Lu  A62^  +  A52  +  ml  rcl  °z  tsln  Y  -  U  cos  yjj 

(D-386) 


C53  *  lLu  A68  +  W  ^A57  "  Lu  A63 )  +  A53 1 


CD-387) 


C54  *  lLu  A69  +  u  (A58  “  Lu  A64  )  +  A54  ^ 


(D-388) 


lLi.  A7 n  +  v  (A59  “  Lu  A65)  +  A55  I 


J55  "  l\,  70 


C56  “  K  CA7i  “  v  aaa)I 


66 


(D-389) 

(D-390) 


Substitution  of  Tilded  Pivot  Forces  Into  Z-Coraponent  of  Moment  Equation 

The  sum  of  the  pivot  forces  in  equation  D-353  is  replaced  by  the  sum  of 
the  tilded  pivot  forces,  as  given  by  equations  D-364,  D-371,  0-378,  and  D-385. 
Then 


where 


p  J.  p  +p  J-p  s>p  +  F  +  F  +p 

xlu  ylu  xlL  ylL  xlu  ylu  xlL  ylL 
"  ^4  +  *75  ♦  +  *76  *2  +  A77  *  +  *78  F12 


74 


7  5 


76 


C37  +  C42  +  C47  +  C52 
Lt  (l  +  u2) 

C3B  +  C43  *  C48  *  C53 

lt  ( 1  +  )7) 

C39  +  C44  +  C49  +  C54 

~  (V  +  u2) 


(D-391) 


(D-392) 


(D-393) 


(D-394) 
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C40  +  C45  *  C50  +  C55 
L_  (TT7T 


C41  +  C46  +  C51  +  C56 


LT  (i  +  m  ) 


CD-395) 


(D-396) 


The  above  is  now  substituted,  together  with  the  thrust  friction  accord¬ 
ing  to  equation  D-348,  into  the  moment  expression  D-353: 


Rbl  F12  +  11  81  al  F12  “  w  pfl  tA72  *  A73  ^ 

,  .  .2  , 

w  P1  H4  *  A75  ^  *  A 76  ^  *  A77  ^  +  A78  F12-‘ 

•  ml  rcl  l°x  8ln  Y  "  °y  cos  +  A60  +  A61  + 


(D-397) 


The  above  is  rearranged  to: 


F12  tRbt  U  8i  ai  +  u  pi  ]  +  p  [pfl  A72  +  p1  A74] 

i  *  *2 

*  p  l Pfl  ^73  +  Pj  a75  j  4*  *  U  Pf  ^7(j  *!*  *  U  PL  A77  + 

-  A6q  +  Aftl  i  -  m1  rcl  [Ox  sin  y  -  0  cos  y] 


(D-398) 


How  consider  the  signs  of  the  various  friction  moments,  recalling  that  a 
reversal  in  the  gear  train  motion  will  cause  a  change  in  the  sign  of  u  in  the 
program.  The  following  moment  components  must  have  negative  signs  during  posi¬ 
tive  rotation  (note  also  that  is  positive): 


l!  '  u  F12  °1  *78 


(D-399) 


since  and  p^  are  positive,  and  A7g  is  a  sum  of  absolute  values. 
2:  -  P  [pfi  A72  +  Pi  A7^] 


(D-400) 


since  pj^  and  pi  are  positive,  while  A72  and  A7^  are  both  absolute  values. 


3: 


“  »*  Pi  ^6  *2 


(D-401) 


since  A75  is  also  a  sum  of  absolute  values. 

The  sign  of  the  term  containing  |  must  be  decided  by  the  sign  of  this 
angular  velocity  only.  Therefore,  the  coefficient  of  friction  must  not  change 
sign  on  motion  reversal,  and  the  expression  takes  the  form: 


-  | u|  [pfl  A73  +  Pj  A75]  $  (D-402) 


The  choice  of  signs  in  the  coefficient  of  the  angular  acceleration  <|>  is 
discussed  in  detail  in  appendix  F  of  reference  2.  This  leads  to  the  computa¬ 
tional  rules  of  equations  D-409  and  D-410  below. 


With  the  above  considerations,  equation  D-398  becomes: 


"  A79  F12  “  A80  ‘  A81  +  "  A82  * 
•• 

m  X1R  *  +  A60  “  ral  rcl  l°x  8in 

where 


A79  a 

Rbl 

-  (1 

al  + 

P  Pf  a78 

A80  ” 

P  [ 

Pfl  A72 

+  PL 

A74  3 

1 

00 

< 

M 

[Pfl  A. 

73  +  p 

1  A7  5  ^ 

A82  ” 

U  P 

1  A76 

^3  “ 

Iwl 

P1  A77 

Further 

X1R  *  A6i  +  Ag3 
•• 

when  $  and  $  have  the  same  signs,  and 
I1R  *  A6i  -  Ag3 


'  °y  coa  ’l  (D-403) 

(D-404) 

(D-405) 

(D-406) 

(D-407) 

CD-408) 

(D-409) 

(D-410) 
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W  s 


■v  V 


'  *  "  d.  ”  *  - 

iv-'! 

r.-.*. 


Expression  for  Contact  Force  F^ 


Equation  D-403  may  be  rewritten  to  obtain  an  expression  for  the  contact 


force  Fios 


F12  " 


"  I1R  *  “  ^1  ♦  “  *82  *  *80  "  *60  +  ml  rcl  l°x  8in  Y  “  °\ 


(n-411) 


Dynamics  of  Gear  and  Pinion  No.  2 


Before  the  force  and  moment  equations  for  gear  And  pinion  no.  2  can  be  writ¬ 
ten,  it  is  first  necessary  to  find  an  expression  for  the  absolute  acceleration  of 
the  gear  and  pinion  pivot  point  Oj ,  which  coincides  with  the  center  of  mass  ^  of 
this  component.  A  top  view  of  gear  and  pinion  no.  2  in  the  mechanism  plane  is 
shown  in  figure  n-9. 


Absolute  Acceleration  of  Gear  and  Pinion  Pivot  0, 


The  absolute  acceleration  of  pivot  point  O2  is  given  by: 


^C^/ground  +  *C/ ground 


(D-412) 


*C/ground  “Absolute  acceleration  of  geometric  center  C  of  mechanism 
plane,  as  given  by  equation  C-4. 


where 


Further 


A-  /  *■  "ui  x  (w  x  <R  n  )  +  J>  x  n- 


w  *■  ui  i+ui  i  +  u)  k 
x  y  z 


u^I+uyJ  +  ijjk 


(D-413) 


(D-414) 


CD-415 ) 
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Figure  D-9.  Top  view  of  gear  and  pinion  no.  2  in  mechanism  plane 


(D-416) 


(D-417) 

(D-418) 

(D-419) 

(D-420) 

(D-421) 

(D-422) 

(D-423) 

(D-424) 

CD-425) 

(D-426) 


Force  Equations  for  Gear  and  Pinion  No.  2 


A  top  view  of  gear  and  pinion  no.  2,  showing  the  contact  forces  F32  and 
Fj2  together  with  their  associated  friction  forces  is  given  in  figure  D-lOa.  A 
free-body  diagram  of  the  pivot  shaft  of  this  component  is  shown  in  figure  D-lOb. 


The  force  equation  is  again  based  on  Newton’s  law,  i.e.: 
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Gear  and  pinion  2 
rotation  is  cw 


T  *  #*  p,,F_  (Thrust  friction 

z2  z2  torque) 


mFx2u 


Z(K) 


f  Y(l) 


Plane  of  forces 

F12  *  F32  and  asso¬ 
ciated  friction  forces 


(L  L  represents  one-half  the 
distance  between  the  pivot  bearings 
in  upper  and  lower  plates.) 


^Fx„L 


Figure  D-lOb.  Gear  and  pinion  no.  2.  Normal  forces,  friction  forces,  and  thrust 
friction  torque  on  pivots. 


The  force  equation  is  again  baaed  on  Newton's  law,  i.e.i 


1  f  ■  "2  *02/«rour>d 


(D-427) 


where 

E  f  ■  Sum  of  the  pivot  forces  as  well  as  the  various  contact 
forces 

m£  *  Mass  of  gear  and  pinion  no.  2 

A.  ,  .  •  Acceleration  of  the  component  center  of  mass,  i.e.,  equa- 

02/ground  tlon  (y_423 


The  full  force  equation  is  now  obtained  with  the  help  of  figures  D-lOa 
and  D-lOb,  as  well  as  equation  B-153  of  reference  1: 


"  F23  n23  ‘  82  11  F23  nN23  +  P12  n12  "  "  81  F12  nN12 
+  Fx2u  3  -  ”  Fx2u  3  +  Fy2u  3  +  •*  Fy2u  3  +  P.2  E 

Fx2L  3  +  y  Fx2L  3  "  Fy2L  3  “  ,J  Fy2L  3 

-  nv  (0X  1  +  0y  ]  +  0Z  k)  (D-428) 

Xn  the  above,  023  and  ^N2j  are  l?lven  by  equations  D-162  and  I)-163,  respectively. 
The  unit  vectors  n^  end  n^^  were  defined  by  equations  D-308  and  D-309.  Appro¬ 
priate  substitution  and  subsequent  separation  into  x  and  y  components  furnishes: 

X-Component  of  Gear  and  Pinion  Force  Equation 

-  F23  sin  (e2  +  «2)-  =2  »f2)  CO.  C9j  +  S2)-  F12  sin  (8,  -  9,  ) 


“  *1  F12  C08  (“l  '  91  1  +  Fx2u  +  vt 


y2u 


-  F 


x2L 


-  p  F 


y2X, 


-  n.n  (0-429) 
2  x 


Y~Component  of  Gear  and  Pinion  Force  Equation 

F23  cos  (b2  +  9,,  )  -  U  s2  F23  sin  )  +  F12  cos  -  0t  ) 

+  "  51  F12  *ln  K  -  91  >  -  “  Fx2u  +  Fv2u  +  "  Fx21j  '  Fy2L  ‘  "2%  (D-430> 
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Z-Componant  of  Gear  and  Pinion  Force  Equations 


This  thrust  force  is  best  expressed  in  tilded  form,  so  that, 

f*2  '  lm2 


(D-431) 


Moment  Equations  for  Gear  and  Pinion  No.  2 


Since  gear  and  pinion  no.  2  represents  a  symmetrical  body  without  pro¬ 
ducts  of  inertia,  its  moment  equations  may  be  expressed  in  terms  of  the  projec¬ 
tile-fixed  X-Y-Z  system  by  the  appropriate  adaption  of  equation  B-13  (app  B) . 

Keeping  in  mind  that  the  angular  velocity  and  the  angular  accelera¬ 
tion  i^,  of  gear  and  pinion  no.  2  with  respect  to  the  projectile,  must  be 

expressed  in  terras  of  the  escape  wheel  angular  velocity  <*>  and  angular  accelera- 
•• 

tion  <(>, 


"  N32  ♦ 


(D-432) 


-  N32  ♦ 


(D-433) 


This  gives  the  moment  equation  the  following  form  (note  that  the  pivot  point 
and  the  center  of  mass  C2  coincide): 


5o2  •  hrf  \  +  1 


0  w  feu  +  N«_  *1  -  I  „  u>.  to  I  I 

z2  y  v  z  32  y2  y  z  1 


+  [:y2  “y  +  Ix2  “x  "k  “  lzl  “5.  H  +  N32  ♦»  J 


+  *x2  (“z  +  N32  ♦)  E 


(D-435) 


The  moment  about  point  O2  is  now  found  with  the  help  of  the  pivot 

shaft  free-body  diagram  of  figure  D-lOb  as  well  as  by  the  adaptation  of  terras  due 
to  and  Fj2»  according  to  equation  B-160  of  reference  1.  Note  that  the  thrust 
torque  (-)  p  p^  k  uses  the  tilded  form  of  Fg2  equation  D-431  in  order  to 
always  make  this  friction  moment  positive,  i.e,  oppose  the  clockwise  rotation  of 
the  component.  The  parameter  pf9  represents  the  thrust  friction  radius,  while  p 
is  the  radius  of  the  pivot  shaft.  Then  2 
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\  '  *23  'b2  E  '  “  *2  F23  *2  E  "  F12  rb2  E  +  “  *1  f12 

+  »  »M  ?l2  k  +  (lu  i  -  »j  i)  *  (rx2a  i  -  «  rx2u  J 

+  (Lu  k  -  p2  J)  x  (Fy2u  j  +  U  Fy2u  I) 

+  (“LL  £  +  p2  *)  *  (-  Fx2L  1  +  y  Fx2L  J ) 

+  (“  *  +  p2  j)  X  (-  Fy2L  I  “  U  Fy2L  I) 

The  above  becomes 


"V,  “  Lu  Fx2u  ~  Lu  Fy2u  +  y  \  Fx2L  ~  LL  Fy2lJ  * 

+  tLu  Fx2u  +  y  Lu  Fy2u  +  LL  Fx2L  +  y  LL  Fy2lJ  ^ 

+  [F23  (Rb2“M82  a2]'F12  ^rb2“W8l 

+  y  pf2  Fz2  +  y  p2  ^Fx2u  +  Fy2u  +  Fx2L  +  Fy2L^  * 

Substitution  of  equation  D-436  into  equation  D-434  yields  the 
component  expressions: 

X-Component  of  Gear  and  Pinion  Moment  Equation 

“  Lu  Fx2u  ■  Hi  Py2u  +  u  ll  Fx2l  '  S,  Fy2L 

‘  Xx2  “x  +  hi  l",  +  N32  -  V  “y  “z 

Y-Component  of  Gear  and  Pinion  Moment  Equation 

L  F  „  +  pL  F  „  +  L,  F  „  +  u  L,  F  ,, 

u  x2u  u  y2u  ix  x2L  L  y2L 


'  Xy2  “y  +  *x2  “*  “z  '  Xz2  “x  H  +  '*32  *> 


(<I  -  «i)  s 


(D-436) 


(D-437a) 
following  moment 

(l)-437b) 

(D-438) 
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Z-Component  of  Gear  and  Pinion  Monant  Equation 

F23  (Si  +  u  “2  *2)  -  F12  <rb2  -  “  ‘l  <dl  ‘  al» 

+  u  pf2  F.2  +  “  p2  (Fx2u  +  Fy2u  +  Fx2L  +  Fy2lJ 

‘  hi  +  N32 


(0-439) 


Simplification  of  Force  Equations  for  Gear  and  Pinion  No.  2 


The  X-component  of  the  force  equation,  i.e.,  equation  D-429,  is  now 
rewritten*. 


where 


”  Fx2u  "  U  Fy2u  +  Fx2L  +  y  Fy2L  “  A84  F23  +  A85  F12  +  A86  (D-440) 


A84  *  “  [8ln  (®2  +  ®2^  +  u  s2  co8  2  +  02^  (D-441) 
a85  “  ”  [sin  (Sj.  “  9i)  “  P  8i  co8  (9i  “  ei)]  (D-442) 
Aae  -  -  m2  nx  (d-443) 


fled  to: 


The  Y-component  of  the  force  expression,  i.e.,  equation  D-430  is  simpli- 


where 


W  Fx2u  "  Fy2u  '  U  Fx2L  +  Fy2L  "  A87  F23  +  A88  F12  +  A89 


A87  *  "  [w  8ln  (02  +  ®2 )  “  cos  (b2  +  92 ) 1 

a88  "  [»  ai  8ln  “  9i)  +  cos  (9i  “  9i)] 

Aq9  -  -  m2  Qy 


(D-444) 

(D-445) 

(0-446) 

(D-447) 
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Simplification  of  Moment  Equations  for  Gear  and  Pinion  No.  2 


Equation  D-437b  for  the  X-component  of  the  moment  may  be  rewritten  as: 


-  “  \  Px2u  +  K  Py2u  -  “  H  Px2L  +  LL  Py2L  -  So  +  Si  ♦  (D-U8> 


where 


SO  ■  '  [IX2  “x  +  “y  \  (S2  "  V  )1 


(D-449) 


Si  '  -  1,2  N32  “y 


(D-450) 


Equation  D-438  for  the  Y-component  of  the  moment  is  written  as  follows: 


Lu  Fx2u  W  Lu  Fy2u  “  LL  Fx2L  "  M  LL  Fy2L  “  A92  +  A93  * 


where 


^2  “  ~  lly2  “y+  \  “z  Kl  "  lz2^ 


(D-451 ) 


(D-452) 


^3  "  1z2  “x  N32 


(D-453) 


Equation  D-439  for  the  Z-component  of  the  moment  remains  as  is. 


Simultaneous  Solution  of  Pivot  Forces  of  Gear  and  Pinion  No.  2 


Equations  D-440,  D-444,  D-448,  and  D-451  are  now  solved  simultaneously 
for  the  pivot  forces.  Therefore, 


-  F  _  -  p  F  .  +  f  +  U  F  ,  B,. 

x2u  y2u  x2L  y2L  21 


uF  -  F  -uF  +F  ■  B 
x2u  y2u  x2L  y2L  22 


-  U  L  F  _  +  L  F  „  '  I!  L.  F  +  L,  F  ,  *  B 

u  x?.u  u  y2u  L  x2h  L  y2L  23 


-L  F„  -  gL  F,  -  L,  F  ,  -  ij  L,  F  ,  B,, 

u  x2u  u  y2u  T.  x2b  L  y2L  24 


where 


(D-454) 


B21  “  ^4  f23  +  A85  f12  +  ^6 


(D-455) 
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B22  "  ^7  F23  +  A88  F12  +  ^9 


(D-456) 


B23  "  A90  +  A9i  ♦  (D-457) 

B24  "  a92  +  \3  ♦  (D-458) 

To  use  the  solutions  of  equation  D-67 ,  equation  D-454  has  to  be  changed 
to  a  form  that  has  the  same  signs  as  this  set  of  expressions.  This  may  be  accom¬ 
plished  by  substituting 

d*  ■  -u  (D-459) 

(This  replaces  A^  *  ^  s,.  in  equation  D-67.)  Equation  D-454  then  becomes: 


(D-460) 


With  the  above  substitution  (i.e.,  equation  D-459)  the  coefficient 
determinant  of  equation  D-460  becomes  according  to  equation  D-75: 

D  -  [(Lu  +  1^  )(l  +  U2  )f  (D-461) 

According  to  equation  D-80,  the  determinant  Dpx2u  now  becomes  with  the 
appropriate  changes: 


F  +  w*  F  +F  -  u*  F  -  B 
x2u  r  y2u  x2L  ^  y2L  21 


U*  F  „  -  F  „  +  d*  F  nT  +  F  „  -  B.„ 

x2u  y2u  x2L  y2L  22 


d*  L  F  +  L  F  +  d*  L  F  „t  +  Lt  F  _t  -  B„, 
u  x2u  u  y2u  L  x2L  L  y2L  23 


-  L  F  +  d*  L  F  -  L  F  +  d*  1>  F  0  *  B 

u  x2u  u  y2u  L  x2L  L  y2L  24 


°F  -  (Lu  +  +  »  )l“  \  hi  +  LL  B22  ■  w  B23  -  B24  1  (D"462) 

x2u 

Now  substitute  for  the  ®2i'8  according  to  equations  D-455  to  D-458: 

°Fx2u  "  (Lu  +  h  Hl  +  >*2  )  {  -  \  [Ag4  F23  +  Ag5  Fl2  +  Ag6  ] 

+  M  \  ^7  F23  +  A8R  F12  +  *89  ^ 
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where 


where 


-  m  [a90  +  Agj  l]  -  [a92  +  A^  $]} 


(D-463) 


After  collecting  of  terms,  the  tllded  force  becomes: 


x2u 


x2u 


4  (l  +  v) 


[C 


57  +  C58 


+  C59  F23  +  C60 


F12]  (D-464) 


'57  * 

1"  LL  A86  +  w  [LL  A89  "  ^0^  “  A92l 

(D-465) 

'58  " 

|y  Agj  +  Agj | 

(D-466) 

59  “ 

lLL  A87  “  A84^ 

(D-467) 

'60  “ 

[LL  A88  “  A85^ 

(D-468) 

According  to  equation  D-90,  D„  with  appropriate  changes  becomes: 

ry2u 

DF  ,  -  lLu  +  +  t-  »  h.  B21  -  \  »22  +  B23  -  “  824]  (D-469) 

y2u 

Substitution  of  equations  D-455  to  D-458  gives: 

DF  ,  -  (L„  +  •tH1  +  “21  I  -  1“  \  [*84  F23  +  *85  F12  +  *86  J 

y2u 

LL  fA87  F23  +  A88  F12  +  AR9  J 

+  [Ago  +  Agi  $]  -  V  [A92  +  A93  $]  CD-470) 


After  appropriate  collecting  of  terms,  the  tllded  force  F  ^  becomes: 


y2u  D  H  (1  +  w2) 


[Ci;i  +  Ca?  ^  +  C63  F23  +  C64  F12 


2>  1  61  62 


(D-471) 
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4X 


C61  ■ 

1"  ll  a89  "  u  (IjL  a86  +  a92^  +  A9()l 

(D-472) 

<*2  " 

1  a9  1  “  S3I 

CD-473) 

O 

a* 

U3 

1 

ao 

< 

+ 

<£■ 

co 

< 

(D-474) 

C64  " 

K  (u  a85  +  A88  5 1 

(D-475) 

According  to  equation  D-100,  D-,  with  the  applicable  changes  becomes: 

rx2L 

DF  ■  (LU  +  h.)!1  +  (>•„  B21  -  «  \  B22  -  “  B23  '  “24)  <D'"6) 

x2L 

Substitute  equations  D-455  to  D-458: 


°Fx2L  "  (Lu  +  LlK1  +  W  )fLu  [A84  F23  +  S5  F12  +  W 
U  Lu  ^A87  F23  +  A88  F12  +  A89 ^ 

“  v  [A90  +  Si  $]  -  f.A92  +  S3  $]} 

After  collecting  of  teims,  the  tilded  force  Fx2l  becomes: 


(D-477) 


TT-  -  ,  -2,  lC65  +  °46  ♦  +  C67  F23  +  C68  P12  I  <D-478) 

lt  U  +  u  ) 


*here 


C65  !~  u  (Lu  A89  +  So^  +  Lu  A86  “  A92  I 

C66  I u  A9i  +  S3 1 
C67  “  K  (A84  *  »•  A87)l 
C68  *  lL.<  (A85  ~  M  A88)l 


(D-479) 

(D-480) 

(D-481) 

(D-482) 


According  to  equation  D-109,  the  determinant  Dp 
becomes : 


after  applicable  adaptation 


°F  "  (Lu  +  +  “  )  (  »  Lu  B21  *  Lu  B22  +  B23  “  11  B24^  <D“A83) 

y2L 

Substitution  of  equations  D-455  to  T>—458  lead  tos 

Dt>  „  ■  K  +  HID  +  "2)  (  “  Lu  t*84  r23  +  *S5  F12  +  *86 1 
y2L 


+  Lu  fA87  F23  +  A88  F12  +  A89  ^ 
+  fA90  +  ^91  ^  “  **  (A92  +  ^3 


CD-484) 


Again,  terms  are  collected  and  an  expression  for  the  tilded  force  Fyov  found. 
Therefore , 

df 

V  •  -P  •  , . 1C69  +  C70  ♦  +  C71  F23  +  C72  F! 2  >  (D-485> 

Uj,  [1  +  M  J 


where 


C69  “  K  A89  +  v  (Lu  A86  “  a92  )  + 
c70  "  lA9i  “  ^  ^3! 

C7 1  “  K  (w  A84  +  A87  )  I 

C72  "  lLu  A85  +  A88 ^ I 


(D-486) 


( D-487 ) 


(D-488) 


(D-489) 


Substitution  of  Tilded  Pivot  Forces  Into  Z-Coraponent  of  Moment  Equation 

Substitution  of  equations  D-431,  D-464,  D-471,  D-478,  and  D-485  into  the 
Z-moment  equation  D-439  is  now  required.  First,  let  the  tilde  forces  be  added: 

?x2u  +  ^y2u  +  ^x2b  +  ^y2L  "  ^4  +  ^5  ^  +  ^6  F23  +  ^7  F12  (D_490> 
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where 


Fur  Cher , 


*94 

Vs 

*96 

*97 


C57 

*  C61  +  C 

'65  * 

% 

4  (i  + 

y  J 

C58 

4-  C62  +  c 

66  + 

C70 

v  (1  + 

2  \ 

y  J 

C59 

+  C63  +  C 

67  + 

C7 1 

4  (i  + 

y2) 

Vc 

l  +  C64  + 

C684 

‘  C72 

v  (1  + 

>2 ) 

let  equation  D-431  be  expressed  as 


?z2  m  *98  '  l"2 
Equation  D-439  then  becomes: 


(D-491) 


(D-492) 


(D-493) 


(D-494) 


(D-495) 


F23  (Rb2  “  P82  F12  (tb2  U  91  ^  "  ai)] 

+  W  pf2  *98  +  w  p2  ^*94  *  *95  ^  +  *96  F23  +  *97  F12  ^ 

-  lz2  (“z  +  N32  ♦)  (0-496) 


F23  K2  '  y  S2  a2  +  P  P2  *96  1  "  F12  ^2  '  P  81  K  "  al  )  '  p  p2  *97  1 
+  V  [pf2  A98  +  p2  A94  ]  ±  p  p2  A95  £  »  Ag9  +  A10q  4>  CD-497) 


where 


99 


z2 


u> 


(0-498) 


A100  "  *z2  N32 


(D-499) 
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Now  consider  again  tha  algna  of  the  friction  moment  terms,  recalling 
that  a  reversal  in  tha  gear  train  motion  will  cause  a  change  of  the  sign  of  p  in 
the  program.  The  component  rotates  normally  in  a  clockwise  direction  and  the 
friction  moments  must  be  positive.  Also  note  that  N32  is  negative. 

The  following  friction  moments  must  be  positive  for  positive  rotation  $ 
of  the  escape  wheel.  (This  implies  a  negative  angular  velocity  for  gear  and 
pinion  no.  2.) 


p  F23  p2  A96  ^A96  is  flum  ah*°^ute  values)  (0-500) 

p  F12  P2  A97  ^A97  8Um  0t  ab80*ute  values)  (D-501) 

P  [p(2  A98  +  P2  ]  (A^  and  A98  are  absolute  values)  (D-502) 

The  sign  of  the  term  containing  $  must  be  determined  by  the  sign  of  ^ 
alone.  When  |  is  positive,  the  gear  and  pinion  no.  2  turns  negatively,  and  the 
friction  moment  must  be  positive.  Therefore,  the  term  must  have  a  positive  sign, 
and  the  absolute  values  of  p  must  be  used: 


+  I  p|  P2  A95  I  (D-503) 

Mote  that  A95  in  an  absolute  value. 

With  the  above  considerations,  equation  D-497  becomes: 

F23  Kl  -  v  82  *2  +  "  °2  A96  1  '  *12  Irb2  *  U  *1  K  '  “l  >  '  “  C2  S?  1 
+  p  [ Pj: 2  ^98  +  p2  ^94  ]  +  I  wl  ?2  ^5  ^  "  A99  +  A100  ^  (D-504) 


Finally,  the  above  is  solved  for  F23! 

A102  F12  “  A103  "  A104  +  +  A99  +  A100  ^ 

p  am  — .  »  -  .  ■  -  . - - 

2‘<  A101 

where 

A101  ■  Rb2  -  U  82  a2  +  (j  02  A96 
A102  "  rb2  "  M  S1  (dl  "  ai  )  ~  V  P2  A97 


(D-505) 


(D-506) 

(D-507) 
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A103  "  W  [p£2  A98  +  o2  AgA  ] 


(D-508) 


A104  "  I  wl  P2  A95 


(D-509) 


Dynamics  of  Combined  System  in  Coupled  Motion  (Applicable  to  Both  Conf igurations) 


As  in  reference  1,  it  is  now  necessary  to  develop  a  single  differential 
equation  for  coupled  motion.  This  is  accomplished  by  first  substituting  equation 
D-411  for  Fj2  Into  equation  D-505  for  F23*  The  latter  expression  is  then  made 
part  of  the  combined  differential  equation  for  the  escapement,  i.e.,  equation  0- 
240  or  D-278.  Therefore  substitution  of  equation  D-411  into  equation  D-505 
gives: 


or 


1  ,102 


F23  "a  lA 
A101  A79 


f  102  f  **  •  *2 

(a  1  [“  *1R  ♦  “  A81  $  “  A«2  ^  '  A! 


82 


'80 


-  A,_  +  m,  r,  (o  sin  y  -  0  cos  y)l 

60  1  cl  v  x  y  J 1 

•  •• 

-  A103  -  A104  ♦+  A99  +  A100  «j>} 


(D-510) 


23  A 


1  f  :  r  102  l  .  it  r  102  ‘81  ,  .  1 
—  Mr  T—  I,p  +  A,nn  J  ♦  A—  +  AinA  j 


A—  A8] 


101 


79 


79 


*2  f  A102  A82  ,  r  ^A80  +  A60  ^  A102  .  ,  i 

+  ♦  I - A— ]  +  l - T7 - A 103  +  V  j 


^9 


A79 


+  ~~  m.  r  .  [O  sin  y  -  0  cos  y)} 
Kyg  1  cl  v  x  1  y  '  ’  ' 


(D-511) 


Th*  above  la  now  substituted  Into  aquation  D-240 

.2 


IA51  *PR  U  ‘  ^9  +  tA51  (A)2  U  +  *PR  ^  ‘  ^9  *4®)  ♦ 


A  A 

+  A51  Ajj  U  ♦  -  2^io^9  {[Al00  A. 


A102  X1R 


79 


4> 


-  [■ 


A.„„  Aa,  .  A1ft0  Aft,  -A102  ^A80  +  A60^ 


k102  "81  .  .  if  102  "82  *2 


'79 


+  Ainz.  ]  ♦  “ 


‘104 

A 


MPWI 

/9 


^  -  [■ 


'79 


+  *103  -  *9?!  +  rf  "l  rcl  (°x  ,ln  1  -  °y  “•  01 


’79 


+  *29  *50  -  *51  (*9  +  *30 1  +  *51  ",  %  K  atn  8  "  K,  c°*  8)  (D"512> 


Finally,  the  above  becomes* 


•2 


A105  *  +  A106  *  +  A107  * 


*108  +  *109  K  aln  T  ‘  V°*  y1 


+  AU0  [\  Sln  P  -  Ky  C0S  0] 


where 


A29  A49  A100  ,  A29  A49  A102  t 
■  Ae,  I„„  U  -  Aon  I_„  -  T  +  "  ^  a  «  1 1R 

*101  101  *79 


U05  “  a51  lpr  u  rt29  ^  zs 


A29  A49  A102  A82 


A106  “  A51  ^A32  U  +  XPR  V)  “  ^9  A48  +  A1(J1  A?9 

,  „  ,  A29  A49  rA102  ^1  ^  A  > 

A107  “  A51  A31  U  A,m  '  *  A79  A104 


101 


(D-513) 


(D-514) 


CD-515) 


(D-516) 


*  The  value  of  the  signum  function  8y  decides  whether  entrance-  or  exit-coupled 
motion  is  described  by  the  differential  equation  D-513. 


3 


r,  ^  v 


A29  A49  rA102  ^0  +  A60  ) 


+  A29  A50  “  A51  *A9  +  A30  ^ 


A29  A49  A102  _ 

Al09  "wsT  1  01 


Ag-  ] 


(n-317) 


(D-518) 


A110  “  A51  mp  rcp 


CD-519) 


Contact  Forces  for  Coupled  Motion 


The  contact  force  F23  is  given  by  equation  D— 51 1 t 
Alll  *  +  Al  12  ^  +  A1 13  *  +  A114 


(D-520) 


where 


.  _  A102  T1R  ,  , 

AU1 - aZ — +  Aioo 


(D-521) 


A102  A82 


(D-522) 


.  A102  A81  ,  ,  ^ 
Ai  13  (  A?9  +  A104  5 


(A  _  +  A,.  )  A... 

v  80  60  ;  102 


"  A 103  +  A99 


CD-523) 


+  ™1  rcl  l°*  8in  Y  '  °y 


(D-524) 


The  contact  force  F^  is  found  with  the  help  of  equation  D-505: 


F.  * 
f  12 


F23  A101  +  A1 1 5  +  A104  *  ”  A100  * 


(D-525) 
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wh«r« 
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Note  that  this  Is  alto  now  indapandent  of  antranca  or  exit  condition 


Now, 

lead*  tot 

A1 17 

whare 

AU7 

a118 

aU9 

a120 


F23  (aq  n-520)  la  aubatltuted  into  the  above  expression.  This 
*+  *118  *2  +  *119  *120  (D-532) 


,  *49  11 1 

A101 

(D-533) 

.  *49  *112 

148  A101 

(D-534) 

49  A113 
aioi 

(D-535) 

(*50^) 

(D-536) 

Contact  Forces  During  Free  Motion 


..  ^  Equation  D-531  may  be  solved  for  the  free  motion  contact  force  Fj,^  once 

$  and  <(>,  for  free  motion,  are  known: 


F 

F23 


♦  +  A48 
A.  ft 


(D— 537  ) 


(The  additional  subscript  F  stands  for  free  motion.) 


Equation  D-505  was  derived  for  gear  and  pinion  set  no.  2.  It  may  now  be  modified 

to  obtain  the  free  motion  contact  force  Fj>j2»  using  free  motion  values  of  Fp23* 

"  • 

<t>,  And  <)>: 


PF12 


fF23  a101  +  A103  +  a10A  $ ~  A99  ~  A100  ♦ 
^102 


(D— 538  ) 


Finally,  with  aquation  0-526 t 


FF23  klOl  +  AU5  +  A104 


A 


F12 


100  6 


102 


(0-539) 


Impact  Equations 


The  Impact  equatlona  may  be  taken  directly  from  reference  1  if  care  la  taken 
to  adjust  all  parameters  to  the  present  notation.  Just  as  in  reference  1,  only 
the  kinematics  relative  to  the  fuse  body  counts,  since  all  other  angular  veloci¬ 
ties  are  common  to  both  pallet  and  escape  wheel. 

For  entrance  impact  the  angle  a,  must  be  used,  while  for  exit  impact  time 
the  angle  a, „  is  applicable.  Entrance  and  exit  impact  equations  are  identical 
(ref  l).  Therefore, 


♦l  frsTOT  °i2  -  S  *!2  )  *  4  ^  *1  +  er  )  Di 

X  A’ ^  +  I  1  n' ^ 

S  1 

a;  -et  (♦,»;  -  a;  ) 


TV 

“STOT  1 


Di 


(0-540) 


(D-541 ) 


where 


^STOT  •  +  lz2  4  +  4 


(D-542) 
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Until  the  time  when  actual  aeroballistic  data  can  he  incorporated  into  pro¬ 
gram  SAF.ROV,  the  following  expressions  for  the  projectile  kinematics  will  be 
used:* 


Spin  Simulation 


Assuming  a  constant  spin  velocity,  obtain  for  the  spin  kinematics: 


•• 


*  0 

(E-l) 

•  _  RPM  x  2  it 

"  npHIE  "  60 

(E-2) 

and 

<t.E  -  PHXE  -  ^  t 

(E-3) 

Precession  Simulation 

Assuming  that  the  precession  velocity  is  also  constant, 
obtained: 

the  following  is 

rP  : 

it 

o 

(E-A) 

where 

•  DPHIE 

%.  m  DPSIt  ‘  -  KP 

(E-5) 

KP  *  Kp,  is  a  divisor  to  obtain  the 
the  soin  velocity 

Kp  r.  100 

\  -  PSIF,  -  ^  t 


precession  velocity  as  a  fraction  of 

(E-6) 
(  E— 7 ) 


*  For  nomenclature  see  appendix  A. 
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Nutation  Simulation 


The  nutation  angle  is  assumed  to  vary  sinusoidally  about  some  initial  angle. 

Then 

0E  -  THET  »  THETIN  +  TVAR  stn  (k  t  )  (E-8) 

where 

THETIN  -8  degrees,  the  initial  cone  angle  (E-9) 

TVAR  "  2  degrees,  the  maximum  change  in  cone  angle  (.E-10) 

«  6  to  8,  multiplier  of  precession  angular  velocity  ijj,  (E-ll) 
to  obtain  maximum  nutation  velocity 

With  the  above 

§R  -  TVAR  *  Kn  *  cos  (Kn  ^  t ) 

0E  -  -  TVAR  *  Kn2  *  sin  (Kn  ^  t ) 

Drag  Deceleration 

"»•  it 

The  deceleration  z  ■  zk  of  the  center  of  mass,  due  to  drag  and  expressed  in 
the  projectile-fixed  system,  is  given  by 

z  =*  DDZ  «  -  386.4  *  10  (E-14) 


(E-12) 
C  E— 13) 
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APPENDIX  F 


COMPUTER  PROGRAM  SAEROV 
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